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ABSTRACT 

An experimental study was conducted to investigate the occurrence 
and growth. of macroscopic instabilities in immiscible displacement of oils 
from porous media. A second objective was to investigate any effects 
polymer additives might have on the instabilities. 

Viscous fingering was photographically studied in a 43.2-cn 
wide x 53.5-cam long x 1.3-cm thick transparent bed packed with beads of 
either 0.47 cm or 0.15 cm in diameter. Dow Corning 200 Fluids with 
viscosities or: 100,° 1000, or 12500°cs were used’ as-the o11- phase. Water, 


various concentrations of 


TQ 


lycercl, and polymer solutions were used as the 
displacing phase. Breakthrough recovery and saturation profiles were 
obtained for flow rates ranging between 0.6 to 15.9 cm/sec. 

A new analytical approach, which eliminates use of relative 
permeabilities was developed to describe the displacement behavior. In 
this approach a saturation dependent viscosity is used in Darcy's law to 
describe flow in the flooded zone. Assuming that the power-law model 
adequately describes viscous properties of the polymer solution, the 
present approach was generalized for polymerfloods. 

It was determined that the advance rate of the zero water 
saturation plane, which corresponded to the tip of the longest finger, 
was linear with time or pore volume injected. Also, the pressure eradient 
of the mixed zone was found to be constant and independent of the position 
of the interface until breakthrough occurred. Good agreement was obtained 
between predicted breakthrough recovery and the available experimental data. 

The predicted saturation distributions were very close to the 


experimental profiles for cases in which severe fingering prevailed. However, 
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disagreement occurred when displacement was more or less uniform. The 
functional form of the average viscosity defined for a mixture of two 
immiscible fluids could be responsible for this disagreement. 

The expression for the breakthrough recovery of a polymer 
flood reflected an effect of the power-law exponent, n. In particular, 
lower values of n, which corresponded to the more non-Newtonian behavior 
of the polymer solution, resulted in lower recovery efficiency. The 
compatibility of the present analysis and the experimental results leads 
to the conclusion that the displacement behavior of a polymer flood 


could be adequately described by a power-law model. 
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NOMENCLATURE 
cross sectional area, on? 
an arbitrary function of x, y, z 
an arbitrary constant 
radius of the cup, am 
portion of the cross sectional area occupied by connate 
water and/or residual oil saturation, an 
an arbitrary function of x, y, z 
an arbitrary constant 


radius of the cylinder, cm 


a 


solute concentration, g cm 2 

an arbitrary function of x, y, z 

net number of counts obtained from scanning 

tortuousity factor or geometric constant, dimensicnless 

solvent volume concentration, cn of solute/an> of solvent 

greatest lateral dimension of the bed, cm 

minimum and maximum pore opening, respectively, cm 

dynamic relative permeability, dimensionless 

mean particle diameter, cn 

effective viscosity ratio, dimensionless 

volumetric recovery efficiency at breakthrough,dimensionless 

friction loss, erg/sec 

fractional Flow rate in terms of relative permeabilites, 
dimensionless 

fractional flow rate, dimensionless 

friction factor defined by Eq'n. (B-25), dimensionless 

2 


gravity acceleration, cm sec 


viscosity level parameter, defined by Eq'n. (B-24) g cm sec 
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ap 


height of the cylinder, cm 
heterogeneity effects parameter 
absolute permeability, one 
effective permeability, cm® 

apparatus constant, cm 

conversion factor defined by Eq'n. (F-1) 


relative permeability, dimensionless 


-torsion bar constant, dyne cm per micron movement of the 


transducer 

length of the bed, cm 

natural logarithn 

viscosity ratio defined by Eq'n. (II-43), dimensionless 
power-law parameter, Agee Secu 

solute molecular weight, g-mole 

viscosity ratio aeiites Oy; Eq'n. (11-64), dimensionless 
mobility ratio, dimensionless | 

shock mobility ratio, dimensionless 

power-law parameter, dimensionless 

Deborah number defined by Eq'n. (C-2), dimensionless 
Reynolds number defined by Eq'ns. (B-27,28), dimensionless 
first partial derivative of z w.r.t. x 


-2 
pressure, dyne cm 


pressure drop, dyne cm”? 
first partial derivative of z w.r.t. y 
total volumetric flow rate, cm? sec”! 


radius of the tube, cm 
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gas constant, g om? sec? g-mole? °K 


radial distance, cm 

-1 
hydraulic radius, cm 

mean radius, cm 

cumulative oil production at breakthrough, % of oil in place 
saturation, dimensionless 

torque transducer sensitivity, micron voit + 
volumetric average water saturation defined by Eq'n. (II-29) 
dimensionless 

temperature, °C 

time, sec 

absolute temperature, °K 

torque, ergs 

movement of torsion head transducer, volts 
superficial velocity, cm sec + 

average interstitial velocity, cm sec + 
average velocity, cm Bec: 

velocity in the z-direction, cm seci 

total volume of the box, cm” 

void volume, cm” 

dimensionless x, i.e. x/L 


position in the bed, cm 


Greek Symbols 


B 


contact angle of the system oil-water-solid material 


empirical constant, dimensionless 


Lx 
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average nominal shear rate, sec + 


shear rate, sec! | 


wall shear rate, sect 


apparent viscosity of the polymer solution defined by 
Eq'n (11-59), poise 

apparent viscosity of the polymer solution based on the 
wall shear rate, poise 

relaxation time, sec 

duration of a process, sec 

wavelength, cn 

viscosity, poise 


S14 159>. =. 


density, g on? 


: ‘ : =3 
density difference betwwen water and oi], g c@m 


summation 


a 


water-oil interfacial tension, dyne am ~ 


¢ 


‘ Sib 
effective or pseudo-interfacial tension, dyne cm ~ 


aie] 
shear stress, dyne cn ~ 


wall shear stress, dyne on? 
porosity, dimensionless 


angular velocity, rad sect 


defined by Eq'n (B-17) 


Subscripts 


average 
breakthrough 


critical 
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CHAPTER I 
INTRODUCTION 
1.1 General 

The paper presented by J. G. Debanné (1) reveals a clear picture 
of the vital role of the oil and petroleum in the new world political 
and wealth order. No matter what politics dictate, the proper recovery 
and use of the petroleum is in the interest of everyone. The period of 
cheap oil is gone and it is widely recognized that global petroleum 
reserves are limited. More effort and research should be conducted to 
find better and more efficient ways of oil recovery. 

The need for more efficient production and higher prices have 
recently caused some considerable improvements in the Secondary Oil 
Recovery techniques. Enhanced recovery schemes are being looked upon 
more seriously. 

The problems and difficulties involved in different schemes are 
common. How is it possible to obtain better volumetric sweep efficiency 
which requires improvement in areal and microscopic sweep efficiencies? 
What factors are responsible for such a low oil recovery efficiency, 
which in the case of heavy oil could be even lower than 20 percent ? 
Presently the problems are known but the answer to how-to-deal-with in 
the desired direction is scarce. For example, it has been known for 
quite a long time that development of a stable oil bank ahead of the 
water zone in the case of water flooding is responsible for high recovery 
efficiency. However, in almost every field operation water enchroach- 
ment into the oil zone has been inferred. This phenomenon is called 
Fingering. The mechanism of this phenomenon is not fully understood. 


Therefore, a well-paved path to eliminate the inception or growth of the 
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2 
fingers cannot be established. No general prediction method of the oil 
recovery at or after breakthrough is presently available. Even the 
effect of very basic parameters such as viscosities of the displacing and 
displaced fluid are not fully understood. Theoretical approaches and 
experimental evidence have always been subject to considerable ambiguity 
and contradiction. The term "due to the complexity of the geometry" is 
often encountered in the literature. Is it possible to resolve this | 
complexity and obtain a clear vision of the process? This question 
Temains to be answered in the future. 

I.2 Viscous Fingering 

The ideal piston-like displacement seldom occurs when an 01il is 
displaced by a less viscous fluid. In 1950 Engelbert and Klinkenberg (2) 
visually observed the intrusion of the water phase into the oil phase in 
a finger shaped manner and coined the term ''Viscous Fingering''. There- 
after many investigators (3-32) farts a on this phenomenon from different 
points of view. It is well documented (33) that in a homogeneous and 
uniform porous medium viscous fingering occurs whenever the mobility ratio 
of invading to original fluid exceeds the value of unity. 

Viscous fingering can be studied from two Aiecone approaches. 
The first approach is to look at the finger as an individual entity of 
the system. This might also be termed a microscopic approach. In this 
approach one expects to obtain information about the inception, growth 
and the number of the fingers. Mathematically speaking this can be 
studied by applying perturbation theory to the flood front (3-12). 

In the second approach one overlooks the fingers and attention 
is focused on the average properties of the bed, say saturation. A 


continuous distribution of saturation is being considered from one end 
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3 
of the bed to the other. A stable or no-finger process corresponds to 
the piston-iike displacement. More precisely saturation is discontinuous 
from one phase to enother and is constant in each phase. The second. 
approach which was mentioned above originates from Buckley-Leverett theory 
(13) therefore, it will be called Buckley-Leverett approach. 
I.2.1 Instability Theory 

It is claimed that a critical velocity exists at which the 
displacement tends to be unstable. It seems from the literature that 
Hill (14) was the first who calculated the critical velocity at which the 
transition zone in the vertical displacement of immiscible liquids became 
unstable. The experimental results of Blackwell, et al. (15) confirmed 
the existence of the critical velocity. A thorough discussion of the 
critical velocity is given by Hawthorne (18). 

Instability and fingering can be demonstrated mathematically by 
perturbation theory. In general, any distrubance can be reoresented as a 
Fourier Series, of which only the least stable terms need be considered. 
This is applicable oniy to the very beginning of the occurence ofthe 
instablities where the displacement equations can be linearized. Chuoke, 
et al. (3) applied the above approach to thie eee displacement and 
showed how the prevailing flow potentials affect the various spectral 
components of a displacement front. Their analysis was based upon a 
paratlel plate model in which a moving interface separated two regions of 
constant, unequal mobilities. However, they generalized the results to 
the porous media by assuming an effective or pseudo-interfacial tension o*. 


They concluded that instability would occur for all rates greater than a 
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critical rate defined by UL» such that, 


Pet es oe : x ae j 
Kye: Kp fec P> - 2,}) g Cos (22°) = 0 (1-1) 


and provided the perturbation contains wave lengths greater than a critical 


wave length given by Ae» 


eZ 


c* 


= ee i ee eS 1eZ 
E ae Ve (1-2) 
Be ENT), 


In (1-2) U is average volumetric velocity, or injected volume of liquid per 
unit time per unit of total area normal to Z. It is sometimes referred to 
as superficial velocity. And Z is the axis directed from fluid 1 to fluid 2 
perpendicular to the macroscopic interface; ZZ“, the angle between the 

Z-axis and the Z<axis which is directed vertically upward. The variable 

o # is an effective interfacial tension, which in the case of displacement 
between parallel plates, equals the ordinary buik interfacial tension. 
Equation (I-2) indicates that wave lengths in a perturbation which are longer 
than ro will grow. . Linear theory also predicts that there exists a wave 


length of maximum instability given by, 
ee om ge he (1-3) 


This maximum instability wave length is expected to be the peak-to-peak 
separation of viscous fingers for the two-dimensional perturbations. A 
similar attempt applicable only to miscible fluids has been undertaken by 


Perrine (4, 5). 
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As Scheidegger (6) points out, the weakness of the Chuoke treat- 
ment (3) is inherent in the fact that in the absence of capillary 
pressure the theory predicts the faster rate of growth for smaller wave- 
lengths. Or in the extreme limit the rate of the growth becomes infinite 
for the wave length approaching zero which is unrealistic. Scheidegger 
(6, 7) included a random heterogeneity spectra for the porous medium which 
he believed was responsible in starting the fingers. He assumed that 
the phenomenon of fingering was not solely governed by the prevailing flow 
potentials, but also by the spectrum of heterogeneities in the porous 
medium. This enabled him to overcome the difficulty encountered by Chuoke, 
et al. (3). Scheidegger (7) revealed that the generation of fingers is 
independent of the displacement velocity. He concluded that since the 
growth of the finger after its inception was entirely controlled by the 
macroscopic flow potentials, therefore, the whole process of fingering 
should also be independent of the displacement velocity. 

Rachford (8) included a transition zone behind the flocd front 
and eliminated the need for introducing pseudo-interfacial tension as in 
the case of Chuoke. Predictions of Rachford's analysis did not correlate 
with.the parallel plate model. His results seem to be very sensitive to 
the shape of the chosen relative permeability and capillary pressure 
curves. From his analysis, Rachford concluded that the effect of flow 
velocity on the onset of instabilities was small and although the 
recovery efficiency decreased the instabilities did not necessarily 
become more severe with increasing viscosity ratio. In fact the viscosity 
ratio appeared to have a damping effect on instability which was in 


contrast with the parallel plate model. 
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Outmans (9, 10) Bencratized the analyses further and included 
non-linear terms in the equations describing conditions at the interface of 
the liquids. His initiative was based on the fact that the non-linear 
terms in hydrodynamics problems cannot always be neglected and are sometimes 
quite important (72). Inclusion of the non-linear terms makes the 
mathematical treatment much more rigorous and difficult to handle. By 
- solving a fourth order approximation, he showed that no particular sinusoidal 
disturbance existed which along its entire length would grow more rapidly 
than any other. This means that the wave length with maximum instability, 
Ay in equation (1-3), predicted by linear theory is at least questionable. 

Perturbation theory and its predictions, if valid, is only 
applicable to the very beginning of the occurence of the instabilities. 

The drastic simplification involved in the theoretical treatment obviously 


makes the predictions questionable. 


Wee wee suck ley Leverett Approacn 


The Buckley-Leverett theory (13) is the first systematic approach, 
in this field. The theory did not predict the growth, size or shape of 
fingers and even could not predict that fingering would occur at ail. How- 
ever, it had the advantage of being easily adapted to obtain the important 
results such as oil recovery and water - oil ratio as a function of pore 
volumes injected. According to the Buckley-Leverett theory the rate of 


advance of each saturation plane Sw in a linear flood is given by, 
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where fi. the fractional flow rate of wetting phase, is defined by, 

f = a (1-5) 
It is clear that a similar equation can be written for the non-wetting 
phase. Equation (I-4) is sometimes referred as the frontal advance 
equation. 

In applying the above theory to the qualitative determination of 
the course of a displacement they fell into the difficulty of the triple- 
valued saturation problem. In their application they considered the case 
where gravity and capillary effects were negligible. To resolve the 
difficulty of the multiple-valued saturation they stated that the real 
saturation profile was discontinuous and the point of discontinuity was 
determined from a material balance. 

In 1951 Holmgren, et al. (17) and Terwilliger, et al. (18) made 
use of the frontal advance equation to describe gas-oil gravity drainage 
performance. In their paper, Terwilliger et al. (18) developed the concept 
of the stabilized and Re eerie! The stabilized zone which 
should not be confused with stabilized flood was defined to be the lower 
region where all the saturation planes moved downstream at the same velocity, 
while in the non-stabilized zone the saturation planes continued to move 
further apart. The stabilized zone is characterized by having a higher ~ 
saturation gradient than the non-stabilized region. They showed that by 
drawing a tangent to the fractional flow curve the saturation jatasthe, up- 
stream end of the stabilized zone could be defined. Later Welge (19) 
showed that construction of the tangent to the fractional flow curve was 
indeed equivalent to the material balance applied by Buckley and Leverett 
(13) to ind the point of the discontinuity. The theoretical and experi- 
mental work of Rapoport, et al.(20) indicated that the non-stabilized 


zone was responsible for the rate dependency of the laboratory water 
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8 
flooding recovery. They concluded that it is the stabilized zone where 


the Buckley-Leverett theory is applicable. 

Most of the investigators, including Buckley and Leverett, have 
contended that the multiple-valued saturation would not have arisen had 
they been able to take capillarity effects into account. The results of 
numerical studies published in papers by Fayers, et al. (21) and 
Havanessian, et al. (22) confirms this original belief. However, the 
ambiguity of the non-capillary displacement remains to be resolved. 

Sheldon, et al. (23) studied the discontinuity of the saturation 
profile from the standpoint of shock formation using the method of 
characteristics and the concept of shocks. Their treatment is indeed 
equivalent to the Buckley-Leverett approach but more logical. Buckley 
and Leverett introduced a discontinuity in the saturation profile to 
resolve the difficulty of the physically impossible multiple-valued 
saturation while Sheldon, et al. considered the muitiple-vaiued satur- 
ation as a point where the characteristic lines intersect and form a shock 
wave. This shock wave proceeds along the bed with a definite speed called 
the shock velocity. The physical picture of their approach is a thin, 
moving region of high saturation gradient in are exact flow 
behaviour is unknown. But does this layer exist at all? 

Generally all attempts have been to rationalize the prediction of 
multi-valued saturation encountered during the course of application of 
the Buckley-Leverett theory. This situation arises mainly due to the use 
of S-shaped fractional flow curves which may be due to neglect of 
capillary forces. It is evident that the experimental relative 
permeability curves may not be a proper choice for the unsteady state 


displacement process. Cardwell (24) seems to be one of the few 
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9 
investigators in the literature who brings up this point. However, he 
does proceed with the same justification as Sheldon, et al. (23). 

Van Meurs, et al. (25) questions the validity of the laboratory 
measured relative permeability curves based on the fact that the same 
distribution of the water and oil present in practice could hardly be 
reproduced in the laboratory during the course of relative permeability 
measurement. They BronOsel an idealized flow mechanism initiated from 
their visual observation (26). They pictured the flow pattern as water 
fingers with immobile protrusions and oil pockets, assuming that the 
fraction of both immobile oil and water in a cross section of the form- 
ation is independent of time and position. Their ultimate results for 
the fractional flow function could be achieved by assuming a linear 
dependency of relative permeabilities with saturation which was later 
derived by Scheidegger, et al. (27). | 

In 1963, Koval (28) introduced the K-factor method for miscible 
displacement which was analogous to the Buckley-Leverett method. He 
viewed the fingering to be brought on by viscosity differences and also 
by channeling and longitudinal dispersion. The effects of the factors 
such as channeling and longitudinal dispersion were included via a 
single parameter, H,, called heterogeneity effects. The transition or 
mixing zone developed due to the solvent fingers was looked upon as an 
effective displacing agent. The ratio of the viscosity of thesoil-te 
that of the assumed effective displacing solvent, was defined as the 
effective viscosity ratio, E. The product of He by E was termed the 
K-factor. In his analysis, the assumption of a single valued viscosity 
characterizing the mixed zone is inherent. This average viscosity was 


calculated from a correlation based on the fourth-root mixing rule, 
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1 Pie VOr.rractions® , Vol.Fraction C 
(Bec) | (4B) % (4c) % 
They also assumed linear dependency of relative permeability with satur- 
ation. Their reason justifying the linear dependency was that they 
assumed no interaction between solvent and oil, that is even if they mix, 
they would retain their individual identities. The linearly dependent 
relative permeabilities with saturation was also used by Verma (29) who 
studied fingering behavior from a statistical point of view in an 
heterogeneous porous medium with capillary pressure. 

Dougherty (30) introduced three more parameters in acon to 
the heterogeneity factor, to account for the fives of dispersive mixing. 
He, too, assumed the relative permeability of each phase to be propcr- 
tional to its volume fraction. Peaceman (31) applied numerical calculations 
to multidimensional miscible displacement taking into account the influence 
of gravity, permeability distribution, diffusion, and fluid viscosities 
and densities. In their numerical calculations they introduced an averaged 


viscosity for the mixed zone given by, 


ae alee Mann ga rae (hea eXd) In uy (I-7) 
where, 
C 
V 
X= oe (1-8) 
C, +e 6u(4 Cy 


and. C, is the solvent volume concentration. The subscripts "o"' and ''so" 

refer to the oil and solvent respectively, and 8 is an empirical constant. 
In a-recent paper, Hagoort (32) modified the Buckley-Leverett 

displacement model by assuming a shock front in the saturation profile 


followed by a zone of gradually changing saturation. Further, he 








: = 
4 7 | 
- . 
£ 
: —-. oP 
p id hee re Ce)’, “ 
BE id) SR OS vt Se ee 
. rae 
. ‘ = ’ ! 
* ‘ ; 
7 ee ae a +s» Late manne Sic Melts alee 
Trev : 7 PS at < wee 
. * ’ Ep 
. aes - f : ta : 
fe) a 2A7) T194n 
oe id 4 Y . % * e “ ~ : ; 
. - . al ; ns > ‘ 5 ¥ 
* ° ~ 
} in! “EPS Pay ie 
. 
* 


“ i 
- . a 4 bs 


Via 
“nets . oi tale 








ip aes. rire rer % "acd Men <A tan oy 


ayer * 





tr 


(2 a. on 


14 
assumed the saturation gradient to be negligible behind the shock front 
so that he could choose a single valued saturation for this zone called 
the shock saturation. Consequently the Buckley-Leverett displacement 
model can be formulated in terms of the Muskat segregated flow model (33) 
with the difference that the displacing phase is a mixture of oil and 
water at the shock saturation. Hagoort introduced a new term called the 
shock mobility ratio, M.. The shock mobility ratio was defined as the 
mobility ratio of the fluids behind and ahead of the shock front. He 
revealed that the displacement was unstable if the shock mobility ratio 
was greater than unity provided the wave-length of the instabilities was 
smaller than the width of the experimental model. He also studied the 
effect of the capillary forces on the instabilities via an energy approach. 
He reasoned that in an unstable displacement the instabilities would give 
rise to an extra production of capillary energy and an extra energy adis- 
sipation caused by viscous forces. Both terms in general depend on the 
wave-length, \. Only those wave-lengths can occur for which the net 
energy dissipation (viscous dissipation minus the capillary energy 
production) is positive. He concluded that the predominant wave-length 
is proportional to the capillary number. It also epee £0, bea 
function of the shape of the capillary pressure curve and displacement 
characteristics such as relative permeabilities and viscosity ratio. 
Finally he investigated the validity of the assumption that there existed 
a relatively small frontal transition zone by requiring \/ ErR eal 
where Lip is the length of capillary transition zone. He found that the 
condition was always fulfilled for moderate values of mobility ratio 


such that, 


M, < o/ (1-9) 
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I.3__ Purpose and Scope 

In the former work (34) flow distribution of Newtonian and non- 
Newtonian fluids in a non-uniform porous bed was studied in hopes of de- 
fining conditions leading to formation of fingers. It was found that 
the microscopic non-homogeneity of the porous bed affected the flow dis- 
tribution significantly. Thus it was possible to postulate that the 
inception of fingers might be due to the microscopic non-homogeneity of 
the porous bed. The effect of polymer solution was to suppress vari- 
ations in the velocity with pore opening to levels below those for 
Newtonian fluids. Therefore, in the light of these findings, a dis- 
placement experiment was performed in a similar bed having an dpieiicial 
large opening or heterogeneity. This opening was more than ten times 
larger than the openings of average pores. In spite of the presence of 
such a large opening the fingering pattern did not seem to change sign- 
ificantly. The prominent fingers which had already formed before the 
front reached the large opening continued on their path. And although 
a new finger could be created at this large opening, it did not seem to 
be the prominent one and was often suppressed after a few rows. These 
visual observations suggested that displacement characteristics such as 
the mobility ratio, viscous, capillary, and gravity forces could be more 
important than the local heterogeneities in governing the fingering 
pattern. This conclusion does not weaken the fact that the local hetero- 
geneities have a significant role in the inception of the fingers. 

After this preliminary study, attention was focused on the dis- 
placement characteristics. The experiments were designed so that 
capillary and gravity forces could be kept negligible and thus attention 


could be focused on the effect of viscosity for the Newtonian and non- 
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Newtonian fluids. 

Modification of the Buckley-Leverett equation based on a power- 
law model and consideration of the feasibility of this in description of a 
polymer flood is also a part of the scope of the present study. Ties 
hoped that the effect of the polymer solution elasticity could be demon- 
strated by comparing the experimental results with predicted behaviour of 


polymer flood from the modified Buckley-Leverett model. 
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CHAPTER IT 
THEORY 

II.1 Method of Characteristics 

The method of characteristics is a powerful tool to deal with a 
system of linear partial differential equations. These are commonly 
encountered in the propagation of the plane waves. Roughly speaking, a 
characteristic is a propagation path along which a physical disturbance 
or entity is propagated. In spite of its long historical background in 
Gas Dynamics, the method of characteristics has been only recently intro- 
duced into the displacement analysis (23, 33). In 1959, Sheldon, et al. (23) 
applied this method to solve the so called Buckley-Leverett partial 
differential equation. He made use of the definition of characteristic 
paths and obtained a differentiai equation which was identical to the 
PCIE VCIeVerect equation (1-4). | 

In the present text the general treatment of the linear differ- 
ential equations given by Abbott (35) will be followed. Consider a 


linear first order partial differential equation given by, 


dbee te Oe aac 204 O Clie) 


where a, b and c are functions of x, y and z and, 


p= 3 (II-2) 
Q- (11-3) 
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The total derivative of z is given by, 


dz =" Pdx + Q ‘dy (11-4) 
Equations (II-1,4) can be combined into the following matrix form, 
a b B C 
= (11-5) 
dx dy Q dz 
It can be proven (35) that along the characteristics the following determinants 


constructed from the above system must be simultaneously identical to zero. 


a b 

= 0 Galo) 
dx dy 
a G 

= 0 (it=7 } 
dx dz 
c b 

= 0 (iT-3) 
dz dy 

It follows that, 

Xue a yen Aaa dzZ = 
Prt S (2759) 


This*is the equationyot the line called the characteristic line,jor more 


precisely, the characteristic of zero order. 
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II, 2 Derivation of Buckley-Leverett Equation 


Let us consider two-phase incompressible flow along the x-axis 
in a porous medium. It is assumed that no flow is occurring transverse 
to x. The positive direction of the x-axis is taken from inlet toward 
outlet of the bed with S defined as volume fraction of the respective 


phases, the continuity equation for this system is given by, 


3d. oom 
ae eh roa (II-10) 
8d, 85, 
x. = =? 6 A at (II-11) 





where subscripts wand o stand generally for the displacing and displaced 


phases or simply water and oil respectively. Also we have, 


= ; a4 
Sy £ So 1 | (1T-12) 
Therefore, 
Q = aie 
ape es qi =o (115 133 


‘That is for a constant injection rate q, 
C Setoh Cee Bales Const. (11-14) 


We can make the following definitions for the fractional flow rates, 
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17 
Neglecting capillary and gravity forces the fractional flow equation in 


terms of relative permeability in Darcy's law, comes out to be, 


fhe I (11-17) 


ie Se Pw 


r,w “o 


where k's are conventional relative permeabilities. 
With the definitions (II-15, 16) Equations (II-10) or (II-11) 


can be transformed into, 


as 3S 

= + an (Gy cee = (II-18) 
where, 

ae) a (II-19) 


The subscript has been omitted because similar equations are produced 
for wetting and non-wetting phases. From Equation (I1-10 and 11) it is 
understood that qe and Io? ECAR Ls and f, are functions of 

x and t. Therefore, in arriving at Equation (II-18) the following trans- 


formation is implicit, 


ipa ey dk (II-19a) 
or, 
Te eerot eaaS 2 
Fae | (II-19b) 


Further assumption that f is a function of S only will be justified later 
when the fractional flow function is derived for the non-capillary 


displacement. 
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On the other hand, 


Se =e (xs t) 


The total derivative of S is given by, 
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Equations (II-18, 21) are combined into the following matrix form, 


Q_ gnc as 
1 $A CS) - 0 
dt dx 2 ds 
Along the characteristics, 
q ft 
1 ie (S) 
= 0 
dt dx 
1 0 
= 0 
dt dS 
0 ace? 
= 0 


ds dx 
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(it =20)) 


(11-21) 


GiieZ7,) 


(II-23) 


(II-24) 


(II-25) 
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19 
From (II-23), 


= a £1(S) (11-26) 


dS = 0 br l273) 
That is, 


S = Const. (11-28) 


which also satisfies equation (11-25). Therefore S is constant along 
the characteristic line, the equation of which is given by (I!I-26). 
Equation (II-26) is identical to the Buckley-Leverett theory given by 
equation (1-4). In derivation of equation (II-26) no assumption was 
made concerning rheological properties of the fluids involved. There- 
fore equation (II-26) would be applicable to the Newtonian as well as 
non-Newtonian fluids. That is by proper choice of a fractional flow 
function, the Buckley-Leverett theory hoids for polymerflood as well as 
waterflood. 
II.3 Volumetric Average Saturation 

The volumetric average saturation, S at or before breakthrough 


is defined as the total volume of fluid injected divided by the total 


volume of the porous bed invaded by the displacing phase, such that 





eee oe (11-29) 


in which Xo is taken as the position of the zero saturation plane. 
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It was noticed from the experimental data given in Appendix (F.4) Table 
(F-7) that the volumetric average saturation pence’ constant till break- 
through occurred. It is understood that since the interface is quite 
flat before flooding is started, no mixed zone or saturation distribu- 
tion exists at the very beginning of the displacement. Therefore 
S. at t = o is undefined. However, when displacement starts fingers are 
created so that mixed zone is developed and S., achieves its fixed value. 
As it is observed in Figures (V-3) to (V-18), the data points lie on the 
eye-fit line from the very beginning of the displacement; therefore, it 
is safe to claim that the transition time for a to achieve its final 
fixed value is very short and insignificant. 


From equation (11-29) for the constant S , 


oe Cares a 
tae oA = CI1-30) 


By comparison of the equations (11-26, 30) it is concluded that, 


= (Ti=41) 





That is the volumetric average saturation is equal to the inverse of the 
derivative of the fractional flow function with respect to S, evaluated 

at Sy = §. Also, it is inherent in the definition of the average 
volumetric saturation that Sy is indeed numerically equal to the volumetric 


-sweep efficiency, E,, at or before breakthrough, 
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Therefore, 
eee | 
ER (II-32) 
£6) 
Ss: = 0 
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II.4 Fractional Flow Function 

From the foregoing analysis it is understood that a properly 
defined fractional flow function with respect to saturation is required. 
First it was decided to proceed with the more commen method of obtaining 
the functional dependency of ae with Sy That is to measure the steady 
state relative permeability curves and calculate the fractional flow 
function from equation (II-17). In this work as in others experiment 
was carried on by similtaneous flow of water and oil measurement of the 
pressure drop at a steady state condition. It can be shown that the 
relative saturation in the bed is not simply the relative flow rates 
except when ol Bw = 1 and linear effective permeabilities with satur- 
ation pertain. During the course of the experiment it was visualiy 
noticed that the flow pattern was significantly different frem that of 
the displacement process. In the case of simultaneous flow each fluid 
tended te have its own path and channeling was quite evident. In other 
words the flow of the two fluids tended to be completely segregated on 
a scale much larger than the bead diameters. It was observed that the 
portion of the total Ayn 48 the bed covered by one phase was not equal 
to the relative flow rate of that phase. For example in the case of 
flowing 50 percent water and 50 percent oil, the portion of the total 
area of the bed occupied by the water was much less than 50 percent. 
Therefore, the saturation is not simply the ratio of the flow rates of 


the two phases. This phenonenon has been noticed and reported in the 


eos 


* 
, 
P 
. 
s 
* 
3 
+ 
i 
.* 
‘ 
fF 
ce 
} 
ss 
5 
. 
" 
be 
af 








- . ‘ ~2e 
F , ra : 
an : i 
i A 
r * ’ P oa mg 
i T we ‘ yi ny te ee | 
= a er — + ge ay 
.. 
> » 
: 7: ct) + OY oe 
7 i o ~ 
¥ 4 « » * = = 
cu" i} 
: 1 wot® [eyed See eee 
‘ ‘ ? Per - a . a a . * . . : "= ‘ ae i 
; a i, « 
4 é A ‘ yy eee, 2 ae = ~~ 7% 
f 7 wre = : in 
ts ; : : 7 
: -z be ae » =o oe 
; | At ojo Temas Se 
(ie * . . - is _. i , - Te : 
2 4 > 
i 
n e “4 ie Lo ‘2 
' rot CREME RS 
~ , oo. F = ae) rae 
~ 
«5 1). .otiitusd- aay 
2 > - ie y } . 
oo + yi « 4 iva @ 
Z iho J wate qth ml aM al 





“ A - ow 
: HN 2s Be SAT ee BSS site 
5 me . z 3 a 
“as Tw) “ * 
5 4 fe é 4s aye ‘- 2 
© Wy soa Tt os ce rah = a 
“ae! a) PSY heel tie Af SET TSG “BO: : 
7 ’ 2 “y= Fy 
» veg 
ti to; wht or Dads bed 
: ‘ z aL Me ~- Sie 
i ‘ uu + 2 - ie, had 4 
1 . e m4 ~ > “ =) * gage 2 oe 
5 ‘ ’ of ¢ 


A. - 1 va cee nef f 
* a2 .= Bh: Pui ; eee 
{ ; x oh 









a Oy 6 5) ie Bee he? hor a a! Byes LS ARS. aad Ven utr 26 
A eg > Se ee zi ; 





; > ia 





Le 
literature. For example, the experimental data reported by 
Richardson et al. (36) all indicate severechanneling. In all the cases 
the oil saturation far from the outlet was much larger than the relative 
flow rates of the two phases. They measured the 011 saturation distri- 
bution by weighing the different segments of the core after each test. 
For the oil and gas relative flow rates of ae = 0.002, 0.008 and 
0.004 the measured oil saturations were reported to be approximately 
60, 64, and 63 percent respectively. 

Jennings (37) used a radioactive tracer technique to determine 
the saturation. In this technique a radioactive component is added tc 
the oil phase so that the oil saturation in the core could be monitored 
with a scintillation counter. This technique eliminates the error in 
the saturation determination due to the unknown hold up but does not 
correct for any effect that channeling might have on pressure drop. 

The phenomenon of the channeling is quite different from end 
effect. End effects are due to the capillary pressure - saturation 
relationship. End effects can be avoided by measuring the pressure far 
from the ends while the independent flow or channeling of the phases 
occurs throughout the experimental bed. | 

In the case of displacement the oil is driven by the displacing 
phase and no independent force exists in the oil phase. Also the dis- 
tribution of the displacing phase across any plane perpendicular to the 
flow direction seems to be in a random fashion. 

The observed dissimilarity of flow patterns forced the author 
to discontinue this approach and to proceed with an alternative method 


to obtain the function dependency of f with Sw Attention was focused 
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on the pressure drop measurements. It was noticed that the overall 
pressure gradient of the mixed zone was almost ‘constant and independent 
of the position of the interface. The overall pressure gradient of the 
mixed zone is defined as the total pressure drop of the mixed zone 
divided by the length of the mixed zone. This pressure gradient was of 
course quite different from that of a single phase. The constant value 
of the overall pressure gradient in the mixed zone suggests that the flow 
regime of the mixed zone is controlled by a single-valued viscosity which 
is neither the viscosity of the displacing nor displaced phase. This 
brings up the possibility that a properly defined average viscosity could 
be sufficient to characterize the flow regime of the displacement process 
when the capillary and gravity forces are negligible. The criterion to 
consider the effect of the capillary forces being negligible was the 
absence of the residual oil saturation. 

The experimental pressure gradient of the mixed zone was intro- 
duced into Darcy's law to calculate the effective viscosity of the mixed 
zone. This calculated effective viscosity from the experimental pressure 
gradient was reasonably close to the arithmetic average of the displacing 


and displaced fluid viscosities, 


der a 
paige el OE . 
He 5 ORISS 
where subscript "'e'’ stands for "effective''. In Appendix H, it is shown 


that equation (II-33) can be directly derived from the integration of the 
local behavior of the bed. The local behavior of the bed is assumed to 
_be characterized by an average viscosity defined for a mixture of two 
immiscible fluids. An Arrhenius type correlation (39), in which the 


exponent is taken as saturation of the two phases is applied, 


ame lei 7 (11-34) 
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24 
where subscript "a" stands for "ayerage''. Further the mixed zone is 
considered to be constructed of a continuous distribution of saturation 
planes, the viscosity of which are given by Equation (1E-34). That is 
every plane of saturation has a single-valued viscosity characterizing 
its flow regime. A similar correlation of (11-34) was used by Peaceman, 
et al. ($1) in their numerical calculation which was given by Equation 
(1-7). Ignoring Sailers, effects, Darcy's law is written in the 


following for, 


KA 
i oe 4T- 
Qen. : ar (11-35) 
W 
KA 
Seen enon © dE 26° 
oO 
i = EKALS icp 2 
daa= ee (II-37) 
a 
and as before, 
LC Gere CL ng caer Cl (11-38) 


where K is the absolute permeability A, and Ay IPeRtBCRINSI LU GlOss 
sectional areas available for flow of the water and oil respectively. it 
is noticed that the steady state relative permeability term usually used 
in applying Darcy's law to unsteady state displacement process has been 
avoided. Instead the relative cross sectional areas of the wetting and 
non-wetting phases have been introduced. In connection with the conven- 
tional terminology, the group KA /u,, will be called "dynamic mobility" 
and the corresponding relative permeability, which will be derived 

later, is termed "dynamic relative permeability". From equation (II-35) 


to (11-38) one obtains, 
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ho a ae 7 (11-39) 
We : | 


In the absence of connate water and residual oil saturation, the total 
cross sectional area of the porous bed, A, is equal to the sum of the 


areas available for the flow of water and oil phases, so that, 


A = AL + a (11-40) 
The pertinent analysis for the case with connate water and/or residual 
oil saturation is given in Appendix (D). Solving for A. and A, from 


equations (II-39, 40) and applying (II-34) yields, 





S 
W 
wie ME = lg (II-41) 
M-1 
Sw 
- NV 
a, = Mita (11-42) 
M-1 


where, M is the viscosity ratio of the oil to the water, 
Me = (II-43) 


From the definition of i” given by equation (II-15), it is easily derived 


W 
Ped es 


(II-44) 
fr M-1 
This justifies the formerly made assumption that with fixed fluids or 
My, ©& 45.4 timctionvof Slonly, 
Ww W 
II.5 Dynamic Relative Permeability 
To obtain the dynamic relative permeability function with respect 


to saturation that corresponds to the preceding analysis one needs to 


write Darcy's law for each phase with dynamic relative permeability 
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term as, | 
Kd . 
quan. svi = (II-45) 
W 
Kd ee A 
q ". to (11-46) 


where d_"s are called dynamic relative permeability to avoid confusion 
with comventional steady state relative permeability. Comparing 
Equations (11-45, 46) with Equations (11-35, 36) respectively and 
introducing the results of Equations (II-41, 42) the expressions for 


dynamic relative permeabilities are, 


S 
WwW 
da can Mie (11-47) 
2 Meet 
S 
W 
d co 7 MGSEM (11-48) 
T,0O Mace 4 


once are just the same as the ratio of areas given by Equations 
(II-41, 42). In contrast to the conventional relative permeabilities, 
the sum of the wetting and non-wetting dynamic relative permeabilities 
is equal to unity. Or in general, when connate water and/or residual 
oil saturation exist, this sum is equal to a constant less than unity. 
~-TI.6 Saturation Distribution 

Differentiating <- from Equation (11-44) with respect to 5, and 


substituting in Equation (II-26) yields, 


eS 
eh M “on M 
— = a — 4 (II-49) 


upon integration with initial condition of Xa 2 05 
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hs 
at OM “ onM or 
x = a pelea ESE (II-50) 
$ Mie, | 


Hence the saturation distribution is obtained by solving Sy from 


Equation (II-50), 


Li eh etl x(M-1) 4A | ‘ 
Sw in M E qt in M (ee 


In plotting Sy vs x one should consider the fact that physically Sy 
can not accept values larger than 1. Therefore Equation (II-51) 


should be modified as following, 


S = dapat fOGed Dine xe <-e 
W SO 
ian 
2 Nba X(M-1)¢A _s 
ee 1 2nM E qt in M | (11-S1a) 
for all xe A 
Sy =i eT 


-II.7 Volumetric Recovery Efficiency at Breakthrough 


Volumetric recovery efficiency at breakthrough, ER: is obtained 
either from Equation (II-32), or directly from the saturation distri- 
bution. By substituting the derivative of fy evaluated at S_ = 0 into 


Equation (11-32), ER is given by, 


eae : 
E. * itm M : (11-52) 


It is understood from Equation (II-52) that the limit of ER is 
unity when M+ 1. However, similar to the case of saturation profile, 
there exists a physical restriction which is not reflected in Equation 


(II-52). This is simply the fact that ER cannot accept values greater 
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28 
than unity from the physical point of view. This physical constraint is 


imposed on the theoretical prediction of Equation (11-52) as following, 


1*for all M< 1 (1I-52a) 


FR 


3 Me- 1 
ER = Min M for all Mi 2 


Also from a simple material balance one obtains, 


L 
L Sw,bt Se aX) + / or dx Ci 55 } 
Il 
where, 
e meedtectiam _c 
Tl = ese 24 = (M-1)46A Cit 54) 


and L is the total length of the bed, X4> the length of the bed which 
has been swept completely and cares is water saturation at breakthrough. 
Substituting Se from Equation (II-51) and performing integration the 
expression for the water saturation at breakthrough is obtained as, 

S AOeEE oii (II-55) 
which is identical to Equation (11-52). 
II.8 Modification of the Theory for the Power-law Model Fluid 


Equation (11-36) remains unchanged while die in Equations (II-34, 


35} must be replaced by apparent viscosity, n, and these become; 


“ Sy 0 (11-56) 
We ee : 
KAY dP 





Ee = - 7 de (II-57) 
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Therefore, the expression for f, would be, 


ake a 
= SS I-58 
a i - (1 ) 


where Te and n both are functions of S. Assuming the power-law model 


to approximate the polymer solution behaviour, n is given by, 


n-l 


n = (| (II-59) 


W 


where, H is the viscosity level parameter given by Equation (B-24). The 
viscosity level parameter does not have the units of viscosity and is a 
factor which accounts for the additional dependency of the superficial 
velocity, U, on K and 9 due to non-Newtonian behaviour of the fluid. 

As it is shown in Appendix B, the viscosity level parameter appears in 
the modified form of Darcy's law obtained from solving the generalized 
flow equation for a power-law model. 


From Equations (11-37, 57, 59) the expression for n becomes, 
Dak py ail 
Ursula: cae (II-60) 
a 
Solving Equations (11-56) and (II-60) for n and Hap their final 


expressions in tems of Sw are obtained, 
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and, 
Sy n (1-"w) 
n-1l n(1-S ) + as n(1-S_) ct Sw (LT-62) 
Meee Ue ) Uy 


Consequently, breakthrough recovery efficiency, Eps can be easily calculated 


from Equation (II-32), 


aan [4 > 1| | 
i kava Me mn M, ae 


UE af is analogous to the viscosity ratio, M, in the case of Newtonian 


fluid and is given by, 
Mo = —°_ (11-64) 


It is evident that for Newtonian fluids or for whicn n = 1, Equation (TI-63) 
becomes identical to the Equation (Ii-52). 

To obtain saturation profiles one needs to know the derivative of 
f, with respect to Sot From Equation (11-58), considering the fact that 


Ha and n both are fimctions of Sa? it is derived, 


_ [ate #8 He SUE ‘ 
ao (11-65) 


ds. jaa =a ed + S ona 


fag 


Integration of Equation (11-26) with initial condition xi = 0 gives, 


L at ) 
x = ar s (11-66) 
where as. is substituted from Equation (II-65) | 
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It is obvious from (II-66) that in this case it is not possible to obtain 
an explicit relation for Sy in terms of x. However, the saturation profile 
can be easily established by assuming values between 0 and 1 for Sy and 
calculating x from Equations (II-65, 66). This procedure has been used to 
obtain the saturation profiles presented in Figures (V-28 to 30) and 


(V-34 to 37). 


II.9 Discussion and Application of the Theory 


Since the first derivative of the fractional flow fimction is only 
a function of Sy? Equation (II-26) indicates that the characteristic curve 
of each saturation plane is a straight line. The assumption of the homo- 
geneity of the bed is, of course, inherent in the foregoing statement. 
Therefore, a plot of x, the position of each saturation plane, v.s. 
qt/oA must yield a set of straight lines the slope of which is equal to the 
first derivative of the fractional flow function evaluated at that specific 
saturation, Depending on the fimctional form of f£(S) three different 
Situations cam be distinguished; 

a) The slopes of the lines are the same which means that 


£'(S) is a constant or say, 
£1(S) = eS (II-67) 
Therefore, 


E{S)ae = mea Seeeeab | (II-68) 
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where a" and b' are constants. 
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FIGURE TI-[ ScHeMATIC REPRESENTATION OF THE 
CHARACTERISTIC LINES FOR THE 
Case oF STABLE DISPLACEMENT. 


The physical meaning of this case is that all the planes of different 


saturation move with the same speed. That is the initial saturation 


distribution profile does not change with time. This corresponds with 


a completely stable displacement. 


b) The slope of the lines increases with decreasing S, which 


means that the second derivative of ''f'' is negative for 


all values of S, 


d’£(S) ce 
as? 


(II-69) 
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From the mathematical point of view, condition (11-69) corresponds with a 


curve of "f* v.s. "S' being concave downward. 
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FIGURE II-2 9 ScuemaTic PEPRESENTATION OF THE 
CHARACTERISTIC LINES FOR THE 
Case oF UNSTABLE DISPLACEMENT, 


The physical meaning of this case is that planes of low saturation move 
faster than planes of higher saturation. That is the length of the mixed 
zone is increasing. This corresponds with an unstable displacement. The 
degree of the instability or the severeness of the fingering can be studied 
from the relative departure of the planes from each other. 

(c) The slope of the lines increases with S, that is the 


* second derivative of f is positive for all values of S or 


oe 0 (II-70) 
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This is a sufficient and necessary condition for the fractional flow curve 


to be concave upward. 
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FIGURE 11-3  ScHEMATIC REPRESENTATION OF 
SHOCK FORMATION, 
In this case, contrary to the former one the planes of higher saturation 
move faster than the planes of lower saturation. At a certain time they 
will intersect each other causing the physically impossible multiple 
saturation case. This is equivalent to the moment when certain of the 


characteristics intersect and a shock or discontinuity in the saturation 


profile is formed. 


The lower portion of the experimental S-type fractional flow curves 
is concave upward which causes the difficulty of shock formation. The 
experimental S-shaped curves are based on steady state relative permeability 
measurements. As thoroughly discussed in section 1.2.2 the validity of 
these measurements in representing the real situation in displacement processes 
is at best questionable. This subject will be later referred to in 


this section. 
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The plane which can be easily studied experimentally is the plane 
of zero saturation. This plane corresponds with the tip of the longest 
finger. From Equation (II-26) and the foregoing discussion, it is understood 


that a plot of the position of this plane, Xo» V-S. qt/o¢A must be a straight 


line. The slope of this line is, 
SG | 
ENGST Eee (iio 1) 
Ow 
} 
(Se 
W 


This would hold for water flood, polymer flood, as well as miscible floods. 
The experimental work published by Benham, et al. (38) verified the above 
linearity prediction for a miscible displacement. They studied the effect 
of distance displaced upon viscous finger length or frontal distortion, D, 
for both open and packed models. They took three or four frames from the 
16-mm movie film for each run which showed the front for three or four 
different average positions in the model. By tracing the outline of the 
front, they determined the mean frontal position from the integrated area of 
the displacing fluid. Indéedyethey converted the traced front to a 
rectangular shape corresponding to the planar frontal shape. Next, they 
divided the model width into 20 equal-width slices, or elements, and found 
the mean frontal position for each element in the same way as before. The 
distance between the mean frontal position of each element and the mean frontal 
position of the whole width was taken as frontal distortion. The absolute 
values of the frontal distortion in each of the 20 elements were arranged 
from maximum to minimum and a number was assigned to each value of frontal 
distortion. This number gave the percentage of the total front having 


distortions equal to or less than the value of that frontal distortion. 
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A plot of frontal distortion as a function of mean frontal, position 
revealed a linear behaviour. Therefore, a linear dependence of frontal 
distortion D, upon the mean frontal Beene ien or distance displaced was 
concluded. The line of "99.99% of front having distortion less than D 
corresponds to the plane of zero saturation in the present study. 

The experimental data indicating the linear progress of the zero 
saturation plane will be presented in Chapter V. As will be observed in 
Figures (V-3) to (V-18) the linearity holds for Newtonian as well as 
non-Newtonian displacing fluids. 

The advancing rate of other saturation planes cannot be deter- 
mined as easily as the zero saturation plane. Therefore some theoreti- 
cal analysis, based on the assumption that the local behaviour of the 
bed could be characterized by an average viscosity for a mixture of two 
immiscible fluids, was undertaken. This made it possible to derive the 
advancing rate of any saturation plane. The advancing rate of the 
different saturation planes was ultimately linked to the breakthrough 
recovery which can be easily checked by ae experimental data. 

The dynamic relative permeability functions derived from the 
present analytical treatment are given by Equations (Ii-47, 48). It 
should be noticed that the viscosity ratio dependence of the dynamic 
relative permeability functions, 1s stronger than the dependence shown 
by published experimental relative permeabilities. In the present 
analysis dynamic relative permeability is nothing more than the relative 
flow areas occupied by each of the phases. It has already been mentioned 
in section (II-4) that in the usual experiment for determinations of 
relative permeability by simultaneous, steady-state flow of the two 


fluids, strong channeling was observed when the viscosity ratio was 
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37 
significantly different from unity. Noting that when the flow area of 
the wetting phase decreases because of channeling it is possible that 
the compensating effects of area change in the non-wetting phase may 
mask the dependence of relative permeability on viscosity ratio. A 
further consequence would be that the usual relative permeability is 
not very sensitive to the details of the flow pattern. This can be 
demonstrated by the following example. Consider a porous bed in which 
we are conducting the usual experiment to measure relative permeability 
by measuring pressure drop at steady state for simultaneous flow of oil 
and water. Suppose that the flow ratio is unity and half of the bed is 
occupied by the water, so that there is no channeling. The total dis- 
Sipation loss in the bed, Ey? is the sum of losses in the oil and water 


phases and is given by, 


ee eae a ate iaa72) 


Be = 10 Loa Gli 735) 
or, 
EY ae LL ew (t1-74) 


Now, supposing that instead of the ideal case, there is significant 
channeling so the water saturation is not equal to its fractional flow 
rate but less than that. For example, assume that water occupies half 


of the former area. In this case, 


Pea ke eee (II-75) 


where prime indicates the new situation. 
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However by proportion of areas and viscosities 


ae, eee (I-76) 
ts aie? a Al) : 
- VO ett) nO rf 3 EW eae, 


From Equations (I-75 to 77), 


meee 20 e 
E bats tata eS (II-78) 
Or, 
E 
V - 33 
Ea Sars (11-79) 


Since total friction dissipation is proportional to measured pressure 


drop of the bed, 


p= (II-79a) 


which is also just the ratio in relative permeability that would be 
eoysvietiiles This is, although the water saturation has decreased by 
half, the change in relative permeability would be approximately twenty 
percent. It is a straight forward algebraic exercise to show that this 
observation holds to the limit of severe channeling or zero flow area 
of the water. This simple calculation reveals that the experimental 
relative permeabilities which are based on pressure drop measurement are 
expected to be somewhat insensitive toward the relative area occupied by 
the channeling fluid. Since the severity of channeling is in turn 
related to the viscosity ratio, the example suggests a similarly weak 
dependence on viscosity ratio. Thus the qualitative difference cited 
above between the dynamic and conventional relative permeabilities may 


be a consequence not only of definition but also because of difference in 
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39 
the experiments which are performed in their determination. 

For polymer flood additional uncertainty arises due to the fact 
that the same pair of relative permeability curves as in the case of 
water flood have been applied (42-45, 56). Based on the conclusion 
derived from the present analysis, both displacing and displaced phase 
dynamic relative permeabilities change if either phase is changed. 
Therefore, by substituting the polymer solution for water as the dis- 
placing phase the former dynamic relative permeability curves do not 
represent the new situation. 

The advantage of the present analysis lies in the fact that the 
analysis can be easily modified to the polymer flood assuming that the 
power-law model describes the behaviour of the fluid. Due to the shear 
dependency of the polymer solution its apparent viscosity changes with 
saturation. The apparent viscosity at each plane of saturation is given 
by Equation (II-61) which clearly shows the error involved in assuming 
a single apparent viscosity applicable all over the bec. This is 
particularly important in the case that terre distribution covers a 
significant portion of the bed. The latter situation is encountered when 
viscous fingering is severe. 

Presb tinocell recovery for a linear polymer flood is given by 
Equation (11-63). It is understood from Equation (II-63) that the recovery 
efficiency depends not only on the viscosity level of the polymer solution 
but also upon the power-law parameter n. Larger n produces larger recovery. 
Therefore, a Newtonian fluid, n = 1, at an equivalent viscosity ratio 
defined by Equation (11-64), is expected to yield higher recovery than a 
polymer solution. It is well known that the polymer with higher molecular 


weight produces smaller n and larger elasticity effect. Therefore, by 
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40 
choosing two different molecular weights of polymers and adjusting 
their concentration to cbtain the same approximate viscosity level, one 
could check the effect of the polymer solution elasticity on the dis- 
placement recovery. If the elastic behaviour of the polymer solution 
has any significant effect on the displacement process, the experi- 
mental recovery efficiency of the higher molecular weight polymer 
solution would deviate from the prediction of Equation (II-63) which is 


based on purely viscous behaviour of the polymer solution. 
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CHAPTER III 


SCALING LAWS 


aa ee General Scaling Laws 

Scaling laws can be obtained either from dimensional analysis 
or from a more mathematical approach called inspectional analysis. 
Dimensional analysis and physical similarity have been used in other fields © 
for a long time and the pertinent theory is adequately given in the 
literature (46 - 49 ). In 1941, Leverett, et al. (50) applied dimensional 
analysis to the two phase flow displacement and studied the similarity 
criteria for radial and also linear displacement. Following Leverett's 
work, Engelberts et al. (2) studied the scaling laws of linear flow systems 
in more detail. They believed in five important scaling factors, two of 
which were considered to guarantee the geometrical similarity such as 
length and inclination of the bed and the rest which are given below were 


considered to preserve the physical similarity, 


Ko AogK Ao g LvK 
y 2 U Hw o Cos « 


In these groups, Ao is the density difference between water and oil and «, 
the contact angle of the oil-water-solid system. Their experimental studies 
showed that breakthrough recoveries depended on most of the factorsnexcept 
Ap gL vK/o Cos =. They concluded ee capillary forces were not expected to 
have substantial effects on water drive processes under field conditions. 

A continuation and extension of Engelberts et al. 's experimental 


work was later undertaken by Croes, et al (51). Rapoport and Leas (20) 





42 


reported the results of scaled experiments on waterflooding and brought 

up a discussion of the stabilized zone which was mentioned in section (I-2.2). 
They revealed that a linear flood can be characterized by means of the 
scaling factor U L ae Later, Rapoport (52) obtained a more general form 

of scaling coefficient by including the water-oil interfacial tension, o, 

and the contact angle, «, His treatment was.based on the more theoretical 
inspectional analysis. He summarized the conditions required for proper 


modelling of water-oil flow processes as: 


1. Geometrical similarity such as boundaries, well distribution, 
wel penetration, etc. 

(gs Initial distribution as well as the succession and distribution 
of operations (water injections and all oil withdrawals) must 
be the same for model and prototype. 

3. Relative permeability functions and the oil-water viscosity 
ratio must be the same for model and prototype. 

4. Capillary pressure functions applying to model and prototype 
must be related to each other either by direct proportion- 
ality or by a general linear transformation. 

5. The ratio of gravitational gradient to capillary pressure 
gradient must be the same for model and prototype. 

6. The ratio of capillary pressure gradient to flowing pressure 


gradient must be the same for model and prototype. 


In 1956, Geertsma, et al. (53) combined the two methods of 
dimensional and inspectional analysis and applied the scaling criteria to 


three types of displacement namely cold water drive, hot water drive and 





43 
solvent injection. Their similarity groups to be considered for cold-water 
drive were reported to be, 


HO KApg : 0 Cos « VKb ; he Ue 
UW i. U uy, U L u 


In most of the practical models the influence of inertial forces is not 
Significant, therefore, the Reynolds' group can be left out of consideration. 
The latter scaling factors were applied by Van Meurs (26) to design his 
transparent three-dimensional model for studying the mechanism of flow 
processes in oil reservoirs. 

So far the relative permeability and capillary pressure relations 
have been eareiatred to be the same functions of saturation in the model and 
its prototype. Perkins, et al. (54) redefined the relative permeabilities 
and saturation in a way to permit the model to have different relative 
permeability and capillary pressure relations. Their new dimensionless 


saturation and relative permeability terms were, 


Si AE: OPS SOG SCN chy Taare 
k a 

ee : W,Yro 
k. Xo 

ee O,CW 


where, Sro is the residual oil saturation, See the connate water saturation, 


Lae the specific permeability to water at the residual 0il saturation and 
> 


K, ee is the specific permeability to oil at the connate-water saturation. 
> 
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44 
The dimensionless parameters produced from the equation of continuity 


combined with Darcy's law are given below, 


Ls be Ki ,ro%o ete: oy ore 
Ly by Ky ww U ub K 


Ky, 10 oe By 

au ew a 
where Lj> L, and L. are some characteristic lengths along the three 
rectangular Cartesian coordinates. 

The principal difference between the scaling criteria outlined 
by Perkins, et al. (54) and those formerly mentioned lies in the modified 
definitions of dimensionless saturations and relative permeabilities. It 
is understood that in this treatment mobility ratio appears to be a 
scaling criterion instead of viscosity ratio. 

Nielsen, et al. (55) numerically studied the performance of a 
hypothetical gas reservoir subject to water drive and its scaled model to 
find the effect of unscalable variables such as exact boundary conditions 
or capillary and relative permeability curves. They revealed that the 
performance of the reservoir and model was very ene to the shape of 
the relative permeability and capillary curves. However, the situation 


did not appear to be that serious for the boundary conditions. 


III.2 Scaling Laws for Unstable Immiscible Displacement 

In the case of studying the fingering phenomenon, some other 
scaling factors seem to be inevitable. During their theoretical and 
experimental investigation of unstable immiscible displacement, Chuoke, 


et al. (3) revealed that a properly scaled model for the study of the viscous 
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45 
fingering phenomenon not only should satisfy the conventional similarity 


groups but also its greatest lateral dimension d, must be greater than a 
critical wave length, roe They argued that if he exceeds the greatest 
lateral dimension, d, of the model the displacement would be stable and 
lead to higher recoveries than in the reservoir. Therefore, they 


included the following additional scaling criteria, 


(; ) 
_ Oy Cn dae Se ree ee Bh re eed gueee Si ee 
: Ne at 
G2 ee eae) d 
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O W 

This in general, need not be exactly the same in both model ana prototype. 
If A /d is either large or small compared to unity in both model and pro- 
totype, either no fingers or a large number of fingers will appear in 
both and the relative effects on behaviour will be the same. Thus only 
14g A Jd is approximately unity would exact scaling be required. The new 
scaling criterion proposed by Chuoke, et al. has been recently supported 
by Hagoort (32) via an energy approach in studying instability of water 
drives in water-wet reservoirs. 

Chuoke, et al's analysis was based on a parallei plate model 
and was challenged by Rachford (8) on the basis that the performance of 
such floods might not properly represent waterfloods in strongly water- 
wet porous systems. Rachford introduced a transition zone behind the 
flood and carried on a first order analysis of the perturbed displacement 
equations numerically. He contradicts the scaling requirement of Chuoke, 
et al. (3) and stated that no additional scaling requirements appeared 


necessary to insure proper treatment of the instability. He believed that 
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46 
since im any rigorously scaled model the same dimensionless differential 
equation describing the prototype, is applicable, the analysis should 
give identical stability predictions for model and prototype. 

From the higher order approximation of the instability theory 
undertaken by Outmans (10) some different criteria for similarity were 
deduced. He distinguished two different cases depending on whether the 
gravity and interfacial tension were negligible. In the case of the 
gravity and interfacial tension being negligible, the model and proto- 
type net only must be kinematically similar but also the initial fie: 
turbance should be identical. This criterion cannot be met unless the 
permeability fluctuations starting the fingers have been scaled down. 

He, therefore, concluded that this type of model study was not a reliable 
method for evaluating sweep efficiency in the reservoir. However, in the 
case when the gravity and interfacial tension were ace Rieaat a complete 
similarity required equality of a dimensionless group somewhat similar to 
the Weber number in hydrodynamics (72). 
III.3 Application of the Scaling Laws 

For the purpose of the present study, in which gravity and 
capillary effects are assumed to be negligible, it is concluded that the 
only criterion to be considered from the general scaling laws is that the 
flow regime mist be in laminar region. That is, Darcy's law must describe 
the flew rate-pressure drop behavior. The scaling laws for fingering 
phenomenon, sree to the kinematic similarity, have not been well 
established. However, the fact that the critical wave length must be 
smaller than the width of the bed seems to be easy and safe to consider. 
It is emphasized that this scaling criterion has been predicted from the 


linear theory (3) as well as an energy approach (32). Therefore, a 
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stronger ground is furnished to believe its validity. 

The width of the packed bed employed in the present investi- 
gation was one foot. A multiple number of fingers even at low flow 
rates were observed in most cases. Therefore, it is safe to claim that 
the ratio of the bed width to the critical wave length was large enough 


to ensure instabilities were not influenced by an inadequate bed width. 
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CHAPTER IV 


EXPERIMENTAL EQUIPMENT AND PROCEDURE 


IV.1 Model Porous Media 

Four different beds were used during the present investigation. 
The first two were composed of regularly spaced matrices of 6 mm diameter 
glass rods ina triangular arrangement. Figure (IV-1) illustrates a unit 
cell of triangular pitch. The last two beds were packed with different 
sized beads. The multi-cylinder matrices produced two-dimensional flow 
while the packed beads provided three dimensional flow at the microscopic 
yevel: 

Displacement tests were performed in the second muliti-cylinder 
matrix (Figure IV-3) and the motion pictures obtained were used to show the 
existence of a linear advance rate of the frontal plane or plane of zero 
water saturation. Results were identical to those presented in Figures 
(V-3) to (V-18). However, no pressure drop measurements were performed 
on this kind of geometry. 

Since the packed beads are presumably a better simulation of flow 
in real porous media except for qualitative description and completeness, 


only the data points of the packed beads will be analysed and presented here. 


IV.1.1 Multi-cylinder Matrices 

The first multi-cylinder matrix was exactly the same one that was 
used in a former study (34). The cylinders were contained in a rectangular 
channel of inside cross section 5.2 cm x 6.5 cm - Figure (IV-2). The 
model contained 40 rows of cylinders with each row alternatively 9 or 10 


cylinders in width. This arrangement left a large space near the wall in 
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FIGURE IV-1 TRIANGULAR-PITCH UNIT 
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50 
each row with nine cylinders. To eliminate by-pass flow at the walls and 


to effect uniform resistance of each row, the 3 mm diameter cylinders were 
placed next to the wall in each of these rows. The minimum and maximum 
pore openings, shown as d, and d, in Figure (IV-1), were approximately 0.5 
and 1 mm , respectively. The volumetric porosity of the bed was 0.36. 

In this work the middle rod of row No. 16 was removed to create an artificial 
large opening at the center. This large Opening was produced to study the 
effect of a large heterogeneity on incipient fingers. 

- The second multi-cylinder matrix was based on the same design as 
the first one but of different dimensions. In this case the cylinders were 
contained in a rectangular channel of inside cross section 1.1 cm x 13.1 cm 
as shown in Figure (IV-3). The shorter depth was chosen to eliminate 
fingering at different depths of the bed. The top and the side wall facing 
the light source were made of 1/2" clear glass plate. To insure that the 
tops of the rods were touching the top glass plate, a 1/2 cm_ sheet of 
spongy material was placed beneath the glass rods. The bottom plate was 
made of 1/4" Aluminium plate. A one-inch thick aluminium plate with 
necessary holes was placed on the top of the spongy material and served to 
hold the cylinders in place. 

To minimize inlet and outlet disturbances, two flow distributors 
25 cm apart were placed before and after the rows of cylinders. The first 
flow distributor was located 1.3 cm. above the inlet and the same distance 
was kept between the second distributor and the outlet. The model contained 
41 rows of cylinders, each row alternatively having 19 or 20 cylinders across. 
The minimum and maximum pore openings were as before, 0.5 and 1 m hes@e ous 
ively. The volumetric porosity, 4, of the bed was 0.30. As mentioned 


before, to eliminate the by-pass channels, smaller cylinders of 3 mm diameter 
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were inserted at the large openings next to the walls. 
IV.1.2 Packed Glass Beads Model 

As shown by Figure (IV-4)the model was made up of two trans- 
parent plates of 3/4 in thick UPVC (Unplasticized Polyvinyl Chloride) 
with dimensions 18 inches wide by 24 inches long. The plates were 
spaced 1/2 inch apart. The internal dimensions of the packed region 
were. 43.2 — cn wide x 55.5 -*cm ions x [25'- cm thick, To obtain a 
uniform injection front, the injection head was composed of three 
inlets. Im addition a flow distributor was placed 1 inch downstream of 
the injection head. This distributor was made up of 1/16 - in thick 
brass plate with holes distributed as shown in Figure (IV-5). A second 
distributor consisting of two brass plates (IV-6) holding a 30 - mesh 
stainless steel screen was placed 1% inches downstream from the first. 
The beads were resting on this screen. To eliminate fluidization and to 
maintain uniform flow distribution to the end of the bed a similar 
rectangular 30 - mesh screen was placed at the top of the packed bed 
just before the outlet head. The outlet head was composed of three out- 
lets to further assist in maintaining wmiforn flow. The distance 
between the outlets and the screen was 1/2 inch. 

To eliminate bulging due to pressure within the bed the back 
and front plates were supported by six 2 - inch wide steel bars. 
These bars, positioned laterally pressed two vertical bars against the 
surface of the sure! and front UPVC plates - refer to Figure (IV-4). 

Oil was pumped into the bed through five inlets located on the 
surface of the front UPVC plate and between the first and the second 


distributor. 
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FIGURE IV-4 PACKED BEADS MODEL 


a3 





54 





FIGURE |IV-5 FLOW DISTRIBUTOR NO1 
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FIGURE |V-6 FLOW DISTRIBUTOR AND 
SCREEN HOLDER 
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To ensure uniform packing of the beads the following procedure 
was used to fill the bed. The bed was moeturoned vertically with the 
outlet head unattached. Three vibrators were connected to the body of the 
bed. Then beads were poured from the top while the vibrators were on. 
Approximately after every inch of beads had been added a T-shaped bar was 
used to tap the beads from the top. After the bed was filled, the outlet 
head was attached and high flow rate water from the domestic tap was run 
from the top of the bed for 24 hours while the vibrators were on. The 
outlet head was then opened and additional beads were added to completely 
fill the bed. 

As mentioned earlier two sizes of beads were studied. The larger 
bead size was 0.47 + 0.05 cm diameter and the smaller was 0.15 + 0.04 cm 
diameter. The smaller size beads were screened through 14 and 9 mesh 
sieves. The volumetric porosities, ¢, of the beds were obtained by direct 
measurement of the amount of fluid necessary to fill the void volumes. 

To do this the Instron was set at a constant low flow rate, q, then the time 
t, required to fill the packed region of the bed was measured by a stop 


watch. Therefore, 


c 
= 
where, V is the total internal volume of the box. The measurement was 
performed several times and at different flow rates. The average volumetric 


porosity obtained for the beds were calculated to be 0.42 and 0.39 for the 


large and small beads respectively with reproducibility se Ue 
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IV.2 Experimental Fluids 

Dow Corning 200 Fluid was used as the oil phase for all of the 
experiments. Dow Corning 200 Fluid is the trade mark of a clear silicone 
liquid available in viscosities from 1 to 100,000 centistokes. Only 
three ranges of viscosity were used for the present study namely, 100, 
1000, and 12,500 cs. Their specific gravities were 0.96, 0.971 and 
0.975 respectively. Oil was usually separated from the water phase by 
decantation and recycled. 

Water, different concentrations of glycerol, and polymer 
solutions were used as the displacing phase. The displacing fluids were 
coloured with a small amount of red food colouring to visualize the 
fingering pattern and enable photographs to be obtained. The displaced 
fluid was colourless. 

The polymer used was a partially hydrolized polyacrylamide 
product of Dow Chemical Co., trade name Separan AP-273. Its solution in 
water shows effects of elasticity. Domestic tap water was used for 
preparation of polymer solutions and care was taken to avoid mechanical 
degredation of the polymer molecules by stirring the solution very gently. 

To compare with Separan AP-273 several additional runs were made 
using a less elastic Pusher 500 solution. Pusher 500 is also a partially 
hydrolized polyacrylamide with smaller molecular weight than Separan 
AP-273. Rheological properties of the experimental fluids are reported 
in Appendix A. 

IV.3 Pumping 

A constant-rate Zenith Gear pump driven by a variable speed trans- 

mission was used for initially filling the bed with oil phase. Filling 


rate was very slow - at a rate of 0.7 cm? sec to eliminate the 
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possibility of trapping air bubbles in the porous bed. The oil phase 
was introduced into the bed from the five taps installed in the trans- 
parent plate just below the packed region and above the injection head. 
The displacing phase was introduced into the bed from the injection head 
by using a single-acting, 6 - inch I.D. X20 - inch hydraulic cylinder 
powered by an Instron model TTEM tester. In displacement tests this 


2 Se thiee. vial Ve Gee 


allowed flows from 0.159 cm 
The hydraulic cylinder was evacuated each time before filling 
with the experimental fluid. This eliminated the air pocket in the 
cylinder so that an instantaneous steady state flow from the Instron 
could be obtained. 
IV.4 Motion Picture Photography 
A slow motion Beauliau movie camera with 25 mm focal length lens 
was used to record the progress of each displacement run on 16 - mm 
coloured movie film. Usually the minimum speed of the camera which was 
2.1 frames/sec was used. In the case of low flow rates the camera was 
stopped periodically for a known length of time as measured by a stop watch. 
The camera was placed at a distance approximately 2 meters from the verti- 
cally positioned bed. 
A light source was located behind the bed. The light source 
was composed of three 18 - inch 32 watts fluorescent lamps type BL 118. 
The lamps were positioned in a 66 - cm x 48 - cm box made up of 22 guage 
galvanized sheets. To obtain as much light as possible a 46 - cm x 61 - am 
mirror was placed behind the lamps. 


IV.5 Film Scanning 





A Wooster Automatic Recording Microdensitometer, Mark II, was 


used for scanning the films to obtain saturation profiles. This 
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58 
instrument was coupled to a Hewlett-Packard eet ‘chart Recorder with an 
Integrator which allowed a continuous measurement of the light intensity. 
The Microdensitometer allows determination of concentration of the coloured 
material by measuring the intensity of light transmitted to a photocell. 

The necessary optical adjustments were performed according to the 
instruction pamphlet supplied with the equipment. These included the 
adjustment of the microscopes and prisms so that the light beam passing 
through the slit plates fell onto the same region of the sensitive surface 
of the photocell. The glass film-carrying plate was set parallel to the 
instrument rails and the desired wedge having suitable gradient was chosen 
and inserted into the wedge carriage. The proper choice of a wedge is the 
wedge which has a minimum density just less than the minimum density of the 
film under test and a miximum density just greater than that of the filn. 

To increase the amount of light falling on the film and wedge, 
lenses of 50 mm focal length were inserted under the film and wedge. The 
lens holders had a pinhole stop to reduce the effects of double scatter or 
to increase resolution. 

To allow scanning in any direction, the glass film - carrying plate 
was cermorted on two carriages which provided two perpendicular movements 
of the film. The position of the carriages could be read from a millimeter 
scale. 

The present experimental work was not originally performed for the 
purpose of the scanning. Consequently, one of the major difficulties was 
the fact that the density of the colour of different batch solutions was 
not the same. Therefore, measurement of the diffused light or saturation 


had to be performed on a comparison basis. 
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Two frames of each displacement run were chosen: one before flood- 
ing or when the bed was full of clean oil, and the other at breakthrough. 
The desired frames were converted into 35 mm negative black-and-white 
slides for the purpose of scanning. The positions of the bars on the 
frame were used to locate positions within the bed. Scanning along 
several lines perpendicular to the flow direction was performed at the 
same position on each frame. The amount of coloured material encountered 
along a line was registered as a number of counts by the recorder. The 
number of counts after flooding were subtracted from the counts before 
flooding to obtain the net number of counts due to flooding. The necessary 
calculations to obtain the saturation profile are given in Appendix F.6. 
The speed of scanning was 9.6 mm/min and the strip chart was run at a 
speed of 1 in/min. 
IV.6 Pressure Drop Measurement. 

Three pressure taps were located on each side wall of the bed. 
The pressure taps were 24.5 cm apart. Pressure drop measurement was per- 
formed by using Validyne DP i5 differential pressure transducers, manu- 
factured by Validyne Engineering Corporation of Northridge, California. 
Diaphragm plates rated at 1, 5, and 20 psi were used depending upon the 
range of pressure drop encountered in the bed. A Model 7100 B Hewlett- 
Packard Strip Chart Recorder with two channels was used to continuously 
record the pressure signal of the transducers. The transducers were 
always positioned below the lowest pressure tap so that a constant static 
head was applied to both sides of the transducers. 

The calibration of transducers and the pertinent data points are 


presented in Appendix E. 
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IV.7 Procedure for Displacement Tests 

The experimental bed was positioned vertically in all cases. 
Prior to a displacement test the bed was always cleaned and dried as 
described later. Therefore, all the displacement runs were at zero 
connate water saturation. The displacing phase was first injected from 
the injection head to a level just below the oil taps. Then the oil taps 
were opened and the oil was pumped in at a very small rate to eliminate 
trapping of air in the bed. In the case of the more viscous oil the 
filling period was as long as three or four hours. After the bed was 
filled the pressure taps were opened and the transducers zeroed. The 
Instron was set at a very small flow rate such as 0.16 cm? sec + of-dis- 
placing fluid to bring the displacing fluid-oil interface to a level 
just below the screen. The Instron was stopped and reset at the desired 
flow rate. Camera and lighting system were set and the experiment was 
then started. 

The displacement test was usually stopped a few minutes after 
breakthrough. The recovered 011 was collected in a cylinder and was 
allowed to separate from the water phase over night. The recovered oil 
was reused in subsequent experiments. After each run the bed was washed 
with water, varsol, alcohol, and again with varsol. After washing was 
completed, air was blown through the bed for several hours until the 
beads were dry. 

To shes aie quality of the pictures taken, and also to illustrate 
the fingering configuration, two series of four frames are presented in 
Figures (IV-7) and (IV-8). The flow rate of Figure (IV-7) is 1.59 cm? sect 


3 = : 
sec +, These pictures are for water 


and that of Figure (IV-8) is 6.35 cm 
displacing 100 cs Dow Comming. Therefore, the oil-water viscosity ratio 


was 96. 
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Photograph A Photograph B 
t=O eyo  OMOlS) 





Photograph C Photograph D 
qt/pAL=0-11 Breakthrough 


FIGURE IV-7 FINGERING IN A LINEAR 


WATERFLOOD,FLOW RATE - 
1.59cmsec (Roll No.11,Run No3) 
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Photograph A Photograph B 
t=O qt/¢AL=0.03 





Photograph C Photograph D 
qt/gAL=0.11 Breakthrough 


FIGURE NM-8 FINGERING IN A LINEAR 
WATERFLOOD,FLOW RATE : 
ee eee (Roll No.11,Run No.2) 


° 
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Photograph A of each Figure was taken at t = 0 in which the interface 
between oil and displacing fluid is completely flat. Photograph B which 
was taken at 0.03 pore volume injected reveals that the fingering pattem 
has already become established. Photograph C is included to present the 
growth of the fingers and show the advance rate of zero saturation plane 
which is discussed in sections VI.2. Finally Photgraph D was taken at 
breakthrough. 

It can be seen from comparing Figures (IV-7) and (IV-8) that the 
width of fingers is reduced at the higher flow rate while the number of 
fingers is increased significantly which helps to achieve higher recovery. 


This point will be considered again in section VI.2 
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CHAPTER V 


RESULTS 

V.1 Pressure Drop 

The measured pressure drop of a single phase Newtonian fluid of 
known viscosity was used to determine the absolute permeabilities of the 
beds. The data points of these pressure drop measurements are reported in 
Appendix F.3. The arithmetic average of the different runs in Table (F-3) 
and (F-4) has been chosen to approximate the permeabilities of the beds 
which were: 


1.97 x 10°" cn’ for the large beads 


D (Os UO! Ge Fass dis aueat SKE 


K 


K 


Figure (V-1) shows the friction factor vs. Reynolds number for 
Newtonian and non-Newtonian solutions. The equations used in this plot are 
equations - (B-25, 27,. 28). The Separan data points of this plot are recorded 
in Tables (F-5,6) and those of the Newtonian fluids in Tables (F-3,4). 

The pressure drop of the mixed zone was determined by subtracting 
the calculated pressure drop of the unflooded portion from the total 
pressure drop measurement. The procedure and the pertinent data points are 
given in Appendix F.4, Table (F-7). Friction factor - Reynolds number 
behaviour for the mixed zone is shown in Figure (V-2) and the correspond- 
ing data points are given in Table (F-8). Reynolds mmbers of the mixed 
zone are. based on the effective viscosity defined by Equation (11-35). ihe 
equations of the friction factor and Reynolds number used in Figure (V-2) 


are the same as those of single phase flow with modification of the viscosity. 


V.2 Frontal Plane Advance 


The position of the zero saturation plane defined in section II-9 is 
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Legend For Figure V-1 


Large Beads 
oO Tap Water 
4a 100 cs Dow Corning 
a 12500 cs Dow Corning 
Vv 0.05% Separan AP-273 
< 0.1% Separan AP-273 
Small Beads 
e Tap Water 
A 100 cs Dow Corning 
= 0.05% Separan AP-273 
Vv 0.1% Separan AP-273 
& 


0.2% Separan AP-273 


Equation (B-26) 


66 


Airietiong bacon yt): Dimensionless 
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7 6 a3 ZA 2 2 
1K. 0 10 10 10 1O me sD iO 
Reynolds Number, Noe Dimensionless 


FIGURE V-1 FRICTION FACTOR VS. REYNOLDS 
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Legend For Figure V-2 


Large Beads 
fe) Water Displacing 100 cs Dow Corning 
= 0.05% Sep. Displacing 100 cs Dow Corning 
4 0.1% Sep. Displacing 100 cs Dow Corning 
0.2% Sep. Displacing 100 cs Dow Coming 


0.2% Sep. Displacing 12500 cs Dow Corning . 


Small Beads 
be Water Displacing 100 cs Dow Corning 
te Glycerol Solutions Displacing 100 cs. Dow Corning 


a 0.05% Sep. Displacing 100 cs Dow Coming 
oo 0.1% Sep. Displacing 100 cs Dow Corning 
+ 0.2% Sep. Displacing 100 cs Dow Corning 
v 0.1% Sep. Displacing 1000 cs Dow Corning 


¢ 0.2% Pusher Displacing 1000 cs Dow Corning 





Equation (B-26) 
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2) 
10 


Friction Factor, fy , Dimensionless 





5x10 10 © 10> 10. 102 2x10 
Reynolds Number,N_ Dimensionless 


FIGURE V-2 FRICTION FACTOR VS. REYNOLDS 
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plotted vs. qt/A¢ in Figures (V-3) to (V-18). The pertinent data points 
and the specification of each run are given in Appendices F.1, 2. 

V.3 Volumetric Recovery Efficiency at Breakthrough 

Volumetric recovery efficiency at breakthrough is plotted vs. 
viscosity ratio for the Newtonian flood in Figure (V-19). The theor- 
etical curve predicted by Equation (II-52) and some experimental data 
points available in the literature are also included in the Figure. 
Figure (V-20) presents the volumetric recovery efficiency at breakthrough 
for the polymerfloods. Also, in the same Figure, the theoretical curves 
Micaulcted by Equation (11-635) for ny =.0. 45, 0758, .0.67 and: 1.0 ‘are given. 
The theoretical curve of n = 1.0, which corresponds to the Newtonian 
flood is identical to that of Figure (V-19). Power law parameters 
pee 0-45, 0.58, and.(.67 correspond’ to 0.2.00" Itand 0.05 percent 
Separan AP-273 solutions respectively. The pertinent data for Figures 
(V- 192620) arespresented’in Tables (F-9 to 13%0 
V.4 Saturation Profile 

Figures (V-21) to (V-37) are plots of water saturation vs. 
Posi cicieex, 1 the bed. | the pertinent data dre presented in dabies 
(F-14, 15). The procedure to obtain the saturation profile from the 
data points of the film scanning is discussed in Appendix F.6, Figures 
(V-21).to f¥-27) are-for waterflood, —(V-31) to (V-35) "tor elyeeror 


solution flood and the remainder are for polymerfloods. 
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| = 20.6354 cm3/sec 
Roll No.8,Run No.5 
=e G<)) Ssle'S ai ee 
Roll No.8,Run No.1 
eae Be hacmansec 
Roll No.8,Run No. 2 
__n— q=6.354cem>/sec 
Roll No 8,RUn No3 
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13 14 
qt/@A,cm 


FIGURE V-3 WATER DISPLACING 100cs DOW 
CORNING , Large Beads 
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; ae Begena: 


ai Se ICE eielee cm/sec 
Koll No.11,Run No.1 
auseee GOL EH 1 Penne 
Roll No.10,Run No.9 
—o— q=0.6354 cm/sec 
Roll No.10,Run No. 8 
_-g-- Q21-5885cm7sec 
RO INO Oeuf Om, 
pe q=3.177cm/sec 


Roll No 10, Run NoiO 


——— © = 6/5 54 Ciieicce 
Roll No.11, Run No.2 


OFM) 24.3 4a pO Oo PO heel 3 
CULO Avierm 


FIGURE V-4 WATER DISPLACING 100 cs DOW 
CORNING, Small Beads 
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Roll No.11, Run No.3 


Roll No10,Run No. 7 





CR ea tae he) AE ay lap or Lae eneey OL ae oe 
qt/¢A ,cm 


FIGURE V-5 WATER DISPLACING 100cs DOW 
CORNING ,4=1-5885 cm >/sec 
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emcee © Che OlVCErO| 
q=0.6354cm~/sec 
RONG la hun NO. Ss 


eOepr Glycerol 


C= ocSScim Sisec 
ROMBNOWCE RURSIING. 


B22 GDaGlyeono! 
q=3177cm “/sec 
Roll No12,Run No4 
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FIGURE} V2Z6 GUYCEROE SOLUTIONS’ DISPLACING 
100cs DOW CORNING, Small Beads 
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—--1— q=1.5885cm?/sec 
VS) Roll No.12,Run No.6 
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FIGURE W27 20cpeG YG EROLSSOLUTION 


DISPLACING 100csS DOW CORNING, 
Small Beads 
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45 
40 
35 
30 

E 25 
20 


1D Eegend: 
° Roll No. 12, Run No.2 
10 o ~©6r RoI! No. 12,Run No. 3 





O 5 
Ooi se aT ae ae) ite, Oh eee Beeb aS 
qt/é¢A,cm 
FIGURE V-8 17cp GLYCEROL SOLUTION 
DISPLACING 100cs DOW CORNING, 
Ga 329/ ACiie 1 SeG 
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45 
4O 
35 


30 


Z ey, vy Legend: 
x Vi [= 918885 cm sec 
Ts Large Beads 
vi Roll Noi@Run No.3 
- pci ea q-1.5885cm2/sec 
Small Beads 
Roll No11, Run No.6 
—s— q=15.885cm 3ysec 
10 Ve Large Beads 
ROIMNGd Onn UnshN@ns 
ie Peer ee i) = 885cm sec 
5 Y/ Small Beads 
Roll No. 11, Run No 7 





o¥ | 
QO. 2 46 8 102, 3416) ie eo 24 26 23 
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FIGURE V-9 0.05% SEPARAN AP-273 
DISPLACING 100cs DOW CORNING 
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20 Legend: 
4 Roll No10,Run No.1 
15 ae RONG LO Run No 4 
o §6©Roll No10,Run No.6 
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FIGURE \V-10 0.05% SEPARAN AP-273 
€ DISPLACING.100cs DOW CORNING, 
q=6.354cm~/sec 
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45 
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20 Legend: 
q=6.354cm/sec 
Large Beads 

15 Roll No.9, Run No. 7 


q=6-354cm sec 
Small Beads 


A Roll No.11, Run Nog 
| q=15.885cmsec 
Ss Large Beads 


Roll No.9, Run No.8 





45 OA OS Ome A Ge Cen er a Oe 
qt/¢A,cm 


FIGURE V-11 0.1% SEPARAN AP-273 DISPLACING 
100cs DOW CORNING 
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Begenc: 3 
—o— q=1.5885cm-~/sec 
Roll No.9,Run No.1 


re q=3.177cm/sec 
Roll No.9,Run No.2 


------- q26-354cm 4/sec 
Roll No.9,Run No.3 


——+— q=15.88 5cm-/sec 


Roll No.9,Run No.4 





O 5 10 We 20 oD) S@) a) 
qt/@A,cm 


FIGURE \V-12 0.2% SEPARAN AP-273 DISPLACING 
100cs DOW CORNING,Large Beads 


~ 
‘ 
- 
7 Ss 
¥" 
? 
* 
Tay 
j ~~ \@ 
* 





a | it ee 








Vein Low lo} 





ial r = ‘ a “ent van Cts sem 1 Amaia 


1 hes =u! y ie 


ann =) 
ie 


ae ea en 


oe 





80 


Small Beads 
Roll No.11,Run NoiO 


Large Beads 
Roll No.9,Run No.2 





O 5 ce) 5 20 BS 30 cS) 
qt/¢A,cm 


FIGURE V-13 0.2% SEPARAN AP-273 DISPLACING 
100cs DOW CORNING,G=3177cm /sec 
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Legend: 


See - 1 G-O.31 77cm Isec 
Roll No.9,Run No.6 


——o— q=0.6354 cm/sec 
Roll No.9,Run Nod 





O 1 } 2) 4 aD 6 i 
qt/¢A,cm 


FIGURE V-14 0.2% SEPARAN AP=273 DISPLACING 
12500cs DOW CORNING,Large Beads 
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FIGURE V-15 SEPARAN AP-273 SOLUTIONS 


DISPLACING 100cs DOW CORNING 
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----2----Q.05%. Sep AP-273 
Roll No.11, Run No 4 
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FIGURE V-16 SEPARAN AP-273 SOLUTIONS 
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O > 10 


AS 20 OD Sh) 35 
qt/¢A,cm 


FIGURE V-17 100cs DOW CORNING DISPLACING 
WATER, q=3-4 cm>/sec ,Large Beads 
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G=© Noo 85cm/sec 
ROH NOiSs RUN NG4 


O2vo Pusher DOO 
=) 03.04 1Chn Sec 
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0.1%. Sep. AP-273 

q-0.15885cm~/sec 

ROlMPNOMssRUneNess 

O1%. Sep. AP=273 

q=0-6354cm/sec 





V Roll No13,Run No.4 
Cee 5 Ss / Se Ton eae 
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FIGURE V-18 POLYMER SOLUTIONS DISPLACING 
1000cs DOW CORNING,Small Beads 
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Legend For Figure V-19 


° Data Points from Everett, et al. (Table F-10) 

i Data Points from Engelberts, et al. (Table F-11) 

a Data Points from Jones-Parra, et al. (Table F-12) 
re) Data Points from Croes, et al. (Table F-13) 

r Present Work: Water Displacing 100 cs Dow Corning 
a Present Work: Glycerol Solutions Displacing 


100 cs Dow Corning 
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Legend For Figure V-20 


Large Beads 
as 0.05% Sep. Displacing 100 cs Dow Corning 
a 0113 “Sep Displacina sl 00sesum pow Corning 
Oo 0.2% Sep. Displacing 100 cs Dow Corning 


Vv 0.2% Sep. Displacing 12500 cs Dow Comming 


Small Beads 
© ~=0.05% Sep. Displacing 100 cs Dow Coming 
Ms 0.1% Sep. Displacing 100 cs Dow Coming 
A 0.1% Sep. Displacing 1000 cs Dow Corning 
= 0.2% Sep. Displacing 100 cs Dow Corning 


* 0.2% Pusher Displacing 1000 cs Dow Corning 


——~ Fquation II-63 
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CHAPTER VI 


DISCUSSION OF THE RESULTS 


Wied Pressure Drop 


It was mentioned in Section (V.1) that the arithmetic average of 
the absolute permeability values in Tables (F-3, 4) was chosen to approximate 
the suite Nine of the beds. In Figure (V1-1) ap/uL.is plotted against 
superficial velocity, U, to show the quality of the data and also the 
applicability of the Darcy's law. The slopes of the lines shown in this 
Figure are the inverse of the corresponding average absolute permeabilities. 

The absolute permeability predicted by Equation (B-21) with the 
commonly used value (see Appendix B) C. = 25/12 = 2.08 was very close to the 


measured value in the case of the small beads, while a different geometrical 


constant is necessary to fit the data for the large beads. This discrepancy 
is probably due to a wall effect. The ratio cf bead depth to diameter of 
beads was approximately three in the case of the large beads. The minimum 


ratio of column diameter to bead size fora packed colum free from wall effects 
is reported (58) to be eight; thus some wall effect is expected with the 
large beads. 

The geometric constant calculated for the large beads based on the 
absolute permeability measurements was o = 3,45, This value was used in 
calculating the Reynolds number from Equations (B-27,.28) for the taree beads. 

The linearity of the friction factor v.s. Reynolds number for 
Newtonian fluids is evident from Figure (V-1) which also shows the applicability 


of Darcy's law. 
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For the three different polymer solutions, upward deviation of the 
data points from the straight line is observed. - This upward deviation 
could be attributed to the effect of elasticity, adsorption and plugging 
as has been discussed in the previous work (34, 59). 

Deviation from linearity can also be demonstrated as in an earlier 
work (68), by plotting f Npe against the Deborah number. This is per- 
formed in Figure (VI.2). The Deborah number is defined and calculated 
for the polymer solutions in Appendix C. The upward deviation at high 
flow rates is emphasized in this Figure. Different concentrations of the 
polymer solution tend to show independent trends of upward deviation. 

That is 0.05% Separan solution shows a significant effect of elasticity 

at a Deborah number of ie while the 0.2% Separan solution does not seem 
to show any elasticity effect even at a Deborah number of unity. The 

early occurence of the upward deviation is not due to adsorption or plugg- 
ing effects. Because, if this was the case, the more concentrated solutions 
would have appeared to show higher upward deviation at an even smaller 
Deborah number. It is observed from Figure (VI-2) that this is not the case 

Figure (V-2) is similar to Figure (V-1) but is for the displacement 
process. The Reynolds number and friction factor were calculated from the 
same equations as for single phase flow, Equation (B-27, 28) and (B-25), 
except for the modification of viscosity. The viscosity used was the 
effective viscosity defined by Equation (II-33). In the case of polymer- 
floods, Ly in Equation (11-33) was replaced by Nyt The arithmetic average 
of the pressure gradients with astrisk in Table (F-7) was used in calcu- 
lating the friction factor of the mixed zone. Figure (V-2) shows that 


a Single valued viscosity correlates the data points satisfactorily. 
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Although the data points of the polymer solutions scatter above the 
Newtonian line, there is no consistent upward deviation as was observed 
for the single phase flow polymer solution data. As a matter of fact, 
the range of the flow rates studied for the displacement tests is not 
sufficient to make. any conclusive statement on this matter. Figure (V-2) 
is an indication that the chosen effective viscosity gives a reasonable 
reproduction of the data points. 
VI.2 Frontal Plane Advance 

The linear advance of the frontal plane, or plane of zero water 
saturation, is predicted by Equation (II-26). The experimental data of 
Figures (V-3) to (V-18) support this prediction. The linearity is observed 
in waterflood, glycerolflood, as well as polymerflood. Figures (V-5, 8 
and 10) show reproducibility of the data points when different runs with 
the same conditions were conducted. 

The slope of the eye-fit lines is equal to the first derivative of 
the fraction flow function evaluated at Sy = 0. It is also understood 
from Equation (11-32) that the inverse of this slope gives the breakthrough 
recovery efficiency. Therefore higher slope indicates less recovery or 
more severefingering. The data points of different flow rates tend to 
behave similarly in a reasonable range of flow rates. The effect of flow 
rate is being shown in Figures (V-3) and (V-4) for the large and small 
beads respectively. It is observed in Figure (V-4) that the data points 
of the different flow rates ranging from q = 0.32 to 3.2 cc/sec could be 
presented by a single line while the highest flow rate q = 6.35 cc/sec 
and smallest q = 0.16 cc/sec, tend to have smaller slope which means 
better efficiency. More stable displacement at very low flow rate is not 


surprising because in the present experiments there exists some small 
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effect of capillary and gravity forces which become important at very 
low flow rate. This tends to stabilize the displacement process. How- 
ever the better recovery at high flow rate was not expected. Reraetine 
to the motion pictures of the different cases, it was noticed that in the 
case of high flow rate many small fingers were created which helped to 
sweep the area more efficiently. A similar behaviour was observed for 
large beads. This behaviour can be observed in the sample photographs 
presented-in Figures (1V-7) and (iV-8). 

dhe effect of flowirate on the slope oS the; frontel.advance lines 
was more pronounced for glycerol solution floods. This is shown in 
Figure (V-7) which reveals that the smaller flow rates have less slope 
corresponding to higher efficiency. The effect of flow rate. on recovery 
is presented in Figure (VI-3) for water and glycerolfloods. 

The effect of viscosity ratio becomes evident by comparing the 
results of the different concentration glycerol solutions and also the 
results for waterflood. Higher viscosity ratio defined by Equation (II-43) 
produces less stable displacement which is in thes same directionsas pre- 
dicted by the mobility ratio considerations (41). This will be further 
explored in the next section. 

From the present data points some effect of the bead size is notice- 
able. The slopes of the eye-fit lines in Figure V-3 are smaller than 
those of the corresponding lines in Figure V-4 which means smaller 
recovery for the smaller beads. This can be verified with the visual 
observation. While fingering bok oh were Similar in the axial direct- 
ion, in the case of the small beads the fingers did not cover the depth 


of the bed. This effectively caused smaller concentration of displacing 
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fluid than if the fingers were the thickness of the bed. For the large 
beads fingers were nearly the thickness of the bed causing a highe 
concentration of displacing fluid. This can be attributed to the effect 
of the front and back walls on finger configuration which forced the 
bed to behave more like a narrow open channel. 

In their study of viscous fingering in miscible displacement, 
Benham, et al. (38) showed that the displacement in an open channel is 
more stable than in a packed bed. Therefore, the behaviour of the bed 
with larger beads is in the expected direction. Indeed, as already 
discussed the pressure drop measurements for the same bed showed some 
wall effect indicating that the bed did not fully represent a packed bed 
free from wall effects. The effect of the bead size is much less 
_ pronounced for polymerflood than waterflood. 

It is understood that the term qt/o¢A is the position of the 
average interface. That is the position of the hypothetical interface 
of a completely stable displacement. Therefore, for a complete piston- 
like displacement, the zero saturation plane would coincide with the 
position of the average interface. Or in terms of the displacement 
figures, a plot of X, VS- qt/oA is a line with slope unity which is the 
limiting case. The limiting case of slope unity was shown experimentally 
by displacing water with oil. The data points of this experiment are 
shown in Figure {V-17). Due to the highly favourable viscosity ratio a 


complete piston-like displacement was expected which is not surprising. 
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It is evident from Figures (V-9 to 16) that all the polymer runs 
have smaller slopes than the waterflood which indicates better efficiency. 
Higher concentration of polymer produces better efficiency. This is 
shown in Figure V-16. However, the increase in recovery from 0.1% 
Separan to 0.2% was not as large as that from 0.05% to 0.1% Separan. 
This suggests a limiting effectiveness of increasing the polymer concen- 
tration. 
V1.3. Volumetric Recovery Efficiency at Breakthrough 

The slope of the frontal plane advance lines, as Aienesecdigs the 
former section, was inversely related to the volumetric recovery efficiency 
at or before breakthrough. To further examine the effect of viscosity 
ratio on the volumetric efficiency, the latter was plotted against 
viscosity ratio in Figures(V-19) and (V-20) for Newtonian and non- 
“Newtonian floods respectively. The theoretically predicted curves of 
Equations (II-52, 63) are also presented. For the waterflood some 
available literature data points have also been included to show the con- 
sistency of the theoretical prediction Rseeieeed in Chepterdl1l. ~Consider- 
ing the fact that the experimental work of the other investigators has 
been performed with quite different and diversified situations such as 
the size and type of the bed, beads, etc., the data points fit the 
theoretical curve in a satisfactory manner. The effect of the flow rate, 
which has not been reflected in the analysis seems to be inevitable. 
This effect has already been discussed in the preceding section. 

The data points of the polymerfloods in Figure (V-20) are without 
exception below the theoretical curve of Newtonian displacing fluids. 
This means that at the same viscosity ratio, as defined by Equation (II-64), 


the recovery of the polymer solution is less than Newtonian fluids. 
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This is in contrast to conclusions that one might arrive at by reviewing 
the abundant literature material on this subject (42, 43, 45, 57). 
Improving the mobility ratio of waterflood via increasing the viscosity 
of the displacing fluid has shown increased recovery efficiency. Also, 
it is well documentated (42, 43, 60-62) that the addition of polymer to 
water may decrease its mobility many times more than would be expected 
from the solution viscosity alone. However, this by itself does not 
justify the claim that the recovery efficiency would also increase as 
many times as the mobility reduction. The earlier literature attributes 
mobility reduction of the polymer solution to viscosity as well as the 
entrapment or adsorption of the polymer molecules, and possibly the 
effect of polymer solution elasticity. Only that portion of mobility 
reduction due to the viscosity can be safely considered for recovery 
increase. Adsorption, pore blockage and any elasticity effect may or 
may not change recovery efficiency. No evidence or direct REE aa 
has been made to show that the recovery efficiency of the polymerflood 
will increase as many times as its mobility has been reduced. Different 
arguments may be necessary for heterogeneous porous media. As shcwn in 
former work (34), due to the effect of elasticity of the polymer solution, 
the effect of heterogeneity seems to be suppressed by adding polymers. 

Lee, et al. (71) in studying a Five-Spot Hele-Shaw model, also 
reported a lower areal sweep at breakthrough with polymer solutions than 
with Newtonian fluids of comparable viscosity. Their analysis only allows 
for an effect of polymer through a variation of viscosity with shear rate 


due to the radial flow. Im linear flow, therefore, they would predict no 
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effect. They did not include the additional dependency of the polymer 


solution apparent viscosity on the saturation distribution as discussed in 


Section II-9. 


VI.4 Saturation Profile 

Saturation profiles presented in Figures (V-21) to (V-37) are the 
profiles obtained at breakthrough. Due to the physical difficulties 
encountered in obtaining the correct data of water saturation, some correction 
had to be applied. The procedure and its justification is discussed in 
Appendix (G.3). 

From the experimental measurements of the saturation profiles 
presented in Figures (V-21) to (V-37) two different types of saturation 
distribution are observed. In the case of waterfloods a sharp decrease in. 
water saturaticn is noticed at the beginning of the bed, then the change in 
Saturation becomes more gradual - Refer to Figures (V-21) to (V-27). A 
sharp decrease in saturation corresponds to high slope in frontal plane 
advance lines which is an indication of severe fingering. The sharp 
decrease in saturation is consistent with the theoretical prediction.  Ihis 
behavior is not restricted to waterflood but occurs whenever the fingering 
1s severe such as Figures (V-28) to (V-30). 

In the case of more uniform sweeping such as Figures (V-31) to 
(V-37), a gradual decrease followed by a sharp decrease of saturation is 
noticed from the experimental data points. There is no apparent major 
difference between Newtonian and non-Newtonian displacing fluid. In this 
group the predicted saturation profiles are qualitatively different from 


the experimental data points. . This difference can not be justified from 
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the error analysis which is presented in Appendix (G.3) Therefore, the 
experimental data points indicate a different regime of saturation distri- 
bution compared to the prediction of the present theoretical analysis. 

The shape of the experimental profiles is similar to the Buckley-Leverett 
profile with frontal transition zone. The assumption of having a mixed 
zone being characterized by a mixture of two immiscible fluids may not hold 
for a uniform displacement. 

As mentioned earlier in section I.2.2, Hagoort (32) showed that 
the analysis based on shock front or a relatively small frontal transition 
zone was not justified for Mp 2 57 where Mp is the conventional mobility 
ratio. Taking the viscosity ratio M, of the present analysis equivalent 
to Mp» it is understood that the displacement tests of Figures (V-21 to 27) 
lie in the group that shock front analysis is not justified, while the 
glycerol solution floods, i.e. Figures (V-31 to 33) remain in the domain of 
the shock front analysis. 

It is concluded that in the case of severe fingering the present 
analysis predicts the saturation distribution very accurately while the 
experimental data points of the more uniform swept displacement tests 


resemble the Buckley-Leverett profile with frontal transition zone. 


VI.5 Effect of Polymer Solution Elasticity on Displacement Process 

The theoretical analysis of Chapter II is based on purely viscous 
behaviour of the displacing fluid. To independently explore the possible 
effect of the polymer solution elasticity on the displacement process, the 
following approach was undertaken. Two polymers, one more elastic than 


the other were chosen. These polymers were Separan AP-273 and Pusher 500. 
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The former one having a molecular weight in the range of 10 million is 
more elastic than the latter with molecular weight approximately 2-3 
million. (Approximate values quoted by Dow Chemical - it is not known 
how they were estimated. Independent normal force measurements support 
the relative magnitudes). 

The concentrations of the two polymer solutions were adjusted 
to obtain a similar viscosity level at shear rates encountered in the 
displacement tests. Their apparent viscosity and shear stress - shear rate 
behaviour is compared in Table A-3 and Figure A-5. The relaxation time 
defined by Bueches' theory (63a) is chosen as a criterion to compare the 


elasticity effect of the two polymer solutions so that, 


where 6 represents the relaxation time. 
Writing Equation (VI-1) for both polymer solutions and dividing 


the respective equations, one obtains, 








u 
= Es a is C 2 | (VI-2) 
Pp Sle 2p 


where the subscripts "'s'' and 'p'" stand for Separan and Pusher respectively. 
In obtaining Equation (VI-2) the viscosity of the solvent has been neglected 
compared to the viscosity of the polymer solution at zero shear rate. 

Relative viscosity, Up defined by the ratio of the polymer 
solution viscosity at zero shear rate to the viscosity of the solvent for 
large values of molecular weight is given (63b) by, 
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By substituting (VI-4) into (VI-2) and simplifying, 


a C 4 7 4.4 

Suk S is F 
tt ie M (VI-5) 
p p : 15D 


In the present investigation the concentration of the Separan solution was 





chosen to be 0.1 percent and that of pusher 0.2 percent. nerefore, by 
taking Mies! gs = 3 one obtains, 
9 4 4.4 
Ss 1 
on = (+) (3) = 7.86 
Pp 


That is, the relaxation time of 0.1% Separan solution is approximately eight 
times larger than the 0.2% Pusher 500 solution while the viscosities are 
approximately equal. If the elasticity of the polymer solution has any 
Significant effect on the displacement behavior it should be evident from 
comparison of the results of these two displacements. If the displacement 
behavior, such as the advancing rate of the zero saturation plane, break- 
through recovery, and saturation profile are similar, it should be safe to 
conclude that the polymer solution elastic behavior has no significant effect 
on the behavior of the cusp taceaene process. 

Figure (V-18) gives a poneenison of the advancing rate of the zero 
saturation plane at two different flow rates. The slopes.of the two eye-fit 


lines are very nearly identical. The small difference in the slope is 
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attributed to the small difference in viscosity level. 

The breakthrough recovery efficiency evaluated from the slope of 
the eye-fit lines of Figure (V-18), also, by using a stop watch to measure 
the time of the injection period till breakthrough is compared in Table (F-9) 


which is reproduced below: 


Table VI-1 Comparison of the Separan 


and Pusher Displacement Recoveries 





E % Of Recoverable Oil | 


Pusher 500 
0.15885 30 28 5S 31 
0.6354 18 18 ily 7 


Obviously the breakthrough recovery efficiencies of Separan are 
Closely in agreement with those of pusher. 

Figure (V-30)compares the breakthrough saturation profiles of the 
twortluids ata flow rate of aq ==) 9026554 cc/sec. The theoretical curve 
predicted by Equation (II-66) withn = 0.53 for 0.1% Separan AP-273 is also 
included in this Figure. The experimental data points of the two different 
displacements are surprisingly close to each other which leads to the 
conclusion that elasticity has little if any effect on the displacement. This 
is in agreement with the earlier conclusion concerning the effect of polymer 


as shown in Figure (V-20). 
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CHAPTER VII 


CONCLUSIONS AND RECOMMENDATIONS 


ViI.J Conclusions 


A new analytical approach has been presented to describe displace- 


ment behavior in porous media. This analysis differs in a major way from 


conventional methods in that relative permeability is eliminated. The 


present analysis and the experimental investigation afford the following 


conclusions. 


V. 


The pressure drop of the mixed zone can be adequately predicted by 
assuming an effective viscosity which is obtained to be equal to an 


arithmetic average of the displacing and displaced fluid viscosities. 


The progress of the zero saturation plane was found to be linear 
With time. This linearity holds for waterflood as well as © 


polymerflood. 


The experimental relative permeability data obtained from steady 
State conditions does not reflect the phenomenon of channeling 
faithfully. This deviation tends to be worse at larger viscosity 


ratios. 


Estimation of the breakthrough recovery from Equations (II-52, 63) 
for water and polymer flood respectively is satisfactory compared to 


predictions based on presently used methods. 


Saturation distributions predicted by Equations (II-51, 66) for 
water and polymer flood respectively are very close to the experi- 


mental profiles in the case of severe fingering. However, the 
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experimentally measured saturation distributions, for the cases of 
viscosity ratios less than about ten, more closely resemble the 


Buckley-Leverett profile with a frontal transition zone. 


Over a range of flow rates recovery efficiency goes through a minimm. 
At low flow rates no growing fingers occur, while at high flow rates 
many fingers are created which help sweep the cross section area of 


the bed more uniformly. 


The displacement behavior of a polymer flood can be adequately described 
by power-law model with power-law indices measured from the viscometric 
behavior of the polymer fluid. The effect of the polymer solution 
elasticity on the displacement behavior was independently checked by 


doing tests with two different molecular weight polymer solutions. These 


showed almost identical behavior. This, of course,. does not overrule 


what might be the controlling effect of polymer solution on inception 
of the fingers due to the local heterogeneity or, in general, on the 


fingering behavior in a heterogeneous porous media. 


The recovery efficiency of the polymer flood is very much higher than 
for a conventional water flood. However, in spite of the general 
expectations, it is less than the case of Newtonian displacing fluid 
with comparable viscosity. This behavior is reflected in Equation 
(II-63) via the power-law exponent, n. A smaller n, which corres- 
ponds to the more non-Newtonian behavior of the polymer solution, 


results in smaller recovery efficiency. 
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Initiation of an individual finger does not guarantee its growth even 
at an unfavorable viscosity ratio. From the photographs both 
acceleration and deceleration of a single finger along its path in the 


bed was noticed. 


VIL.2 Recommendations 


1. 


The present theoretical analysis can be adapted to the polymer slug 
flood to predict the minimum slug size. Experimental slug flood 
should be undertaken to justify the prediction. Modification of the 
present analysis to the polymer slug flood is simply achieved by 
considering a non-Newtonian slug followed and preceded by Newtonian 
fluids. Slug breakdown occurs when the plane of zero water satur- 


ation meets the mixed zone of polymer and oii. 


Inclusion of the gravity and capillary forces in the present analysis 


could give an interesting insight of the gravity and capillary effect 


if the resultant equations are amenable to the analytical treatment. 


The present approach should be adapted to the radial flow case and 
ultimately to pattern flooding, especially five-spot patterns, and 


comparisons made with the available literature data. 


An experimental study of the different kind of polymers with different 
degrees of elastic effect performed on a more heterogeneous porous bed 
might establish useful information about the effect of the polymer 


solution elasticity on displacement behavior. 


This study did not reveal a large effect of bead size, however, a more 


systematic approach to study the effect of bead size, if any, on 
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fingering phenomenon would be useful. Studies by reducing bead size 
should be performed cautiously so as not to include the effect of 
Capillary forces. Alternatively a proper treatment to consider the 


capillary effects must be undertaken. 
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APPENDIX A 131 
RHEOLOGICAL PROPERTIES OF THE 


EXPERIMENTAL SOLUTIONS 


A WEISSENBERG R-18 Rheogoniometer with cup and cylinder arrangement 
was used to measure viscosity of the glycerol and polymer solutions. 
Consider two vertical coaxial cylinders, the outer one of which is rotating 
with an angular velocity 2 by application of a torque To: This torque 
must balance the torque exerted by the fluid on the inner face of cup so 


that, 


ft 2 : 
i = 2rb, h TR (A-1) 


where, b, represents the radius of the cup, h, the height of the cylinder in 
contact with the fluid, and Tp is wall shear stress. 
| The torque, Th is determined by measuring the torsion in a T-bar 
suspending the cylinder, so that, 
qT, a 8 


C -7? 
mks UG (A-2) 


where, Sj represents the torque transducer sensitivity which was calibrated 
to be equal to 1.008 micron/volt, At, movement of torsion head transducer in 
volts, and K.; torsion bar constant which was equal to 106.6 dyne cm per 
micron movement of the transducer. The dimensions of the cup and cylinders 
were as follows: 


radius of the cup, ES hes WA Sein 


4) 


radius of the cylinder, by =) 2n5°cm 


height of the cylinder, h = 5 a@m 
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Therefore, equations (A-1, 2) become, 


| 
i] 


196.35 Tp ergs (A-4) 


ri 
I 


107.455 At. eres (A-5) 


For Newtonian fluids, it can be shown (64, 65a) that the torque 


is related to the viscosity of the fluid by, 
T= Kun (A-6) 


where Ky is a constant for any given apparatus and equal to, 


4nha,“b, : 
K. = po = 2286.49 cm (A-7) 
: 


oqeen 7 


From equations (A-5, 6) uw for Newtonian fluid is calculated by, 


Les 9 02 047 aa poise (A-8) 


shear rate of Newtonian fluid in Couette flow is given (64) by, 
2 
2a, by Q 


y= — (A-9) 
Cr 


where eee is the mean radius, 
L 


ay he 
& ik us ci 
Therefore, 
yie=210157 @, aseeat (A-11) 


The shear rate of a power-law fluid in Couette flow is given (63C) by, 


a9) 
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(A-12) 
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The shear stress-shear rate relationship of the power-law model 
is given by equation (B-9). Therefore Dy plotting log 1 se vlog @,. the 
power-law exponent, n, is obtained from the Slope. The Y-intercept of 


such a plot at 2 = 1 is equal to, 


n 
Y-interncept =. m ae (A-15) 
1 n 
n iC 
1 
from which m can be determined. Or, by knowing n, y is calculated from 


equation (A-12). A plot of log t v.s. log y is constructed with the slope 
equal to n and Y-intercept at 7 = 1 equal to m. 
The apparent viscosity of the polymer solutions is calculated 


from the following equation, 
ee ay (A-14) 


Plots of shear stress v.s. shear rate are given in Figures (A-1), 
(A-3) and (A-5). Figure (A-1) is for the Glycerol solutions, while the 
viscometric behaviour of polymer solutions is presented in Figures (A-3) 
and (A-5). Figures (A-2) and (A-4) are plots of t v.s. 2 used to obtain 
n for the polymers. All of the data points of the viscometric measurements 
are included in Tables (A-1 to 3). 

It should be mentioned that most of the polymer floods were Ger 
formed on Separan solutions with rheological properties given by Figure A-3. 
A new batch of 0.1% Separan solution was prepared for the last few displace- 
Mente Cestseto compare with 0.2% pusher 500. The rheological properties of 
this new solution with those of pusher 3900 are presented in Figures (A-4, 5) 


and Table (A-3). 
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VISCOMETRIC MEASUREMENTS OF CLYCEROL 


SOLUTIONS 





Goa Glycerol Solution, T° OC 
1 2.84 3.14 1.05 sy O17 
2 9.0 9.96 235 1288 aly 
3 28.4 31.44 10.6 5.80 Oui 
4 56.8 62.87 21.8 1103 Oly 
5 9.0 9.96 3.4 1.86 Ooig 
71% Glycerol Solution, T= 20°C 
1 284 Sill 1.6 0.88 0725 
2 OO: 9.96 pal 2.79 0.25 
3 28.4 31.44 16.0 8.76 0.25 
4 56.8 62.87 Bos 17.68 0.26 
5 2.84 oud 1.6 0.88 0.25 
71.5% Glycerol Solution, T = 20.5°C 
1 2.84 3.14 165 0.9 0.26 
2 9.0 9.96 5.30 2.9 0.26 
3 28.4 31.44 16.6 9.08 0.26 
4 56.8 62.87 3507 18.44 O27 
5 2.84 3.14 1.65 0.9 0.26 
72.5% Glycerol Solution, T = 20°C 
1 2.84 3.14 1.8 0.99 0.28 
Z 9.0 9.96 518 518 0.29 
3 28.4 31.44 17.6 9.63 0.28 
4 56.8 62.87 36.8 20.14 0.29 
5 2.84 3.14 ee 0.99 0.28 
77% Glycerol Solution, T = 20°C 
1 2.84 3.14 25 1.78 0.51 
Z 9.0 9.96 OES as 0.52 
3 28.4 31.44 32.9 18.0 Ous7 
4 90 99.62 105 57.46 0.52 
5 2.84 3.14 3.25 1.78 0.51 
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VISCOMETRIC MEASUREMENTS OF SEPARAN 


APE= 9275 SOLUTIONS 
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APPENDIX B 


MODIFIED FRICTION FACTOR 


AND REYNOLDS NUMBER 


The capillary tube model, has been successfully applied to porous 
media to predict the pressure drop or friction factor-Reynolds number 
relation of purely viscous flow (40,66). The model is based on cylindrical 
poiseuille flow using hydraulic radius of the medium instead of radius of 
the cylinder. Consider the flow of a fluid through a circular tube in the 


Z-direction as shown in Figure (B-1). 





FIGURE B-1 = Flow IN CYLINDRICAL TUBE 
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where q is volumetric flow rate, u Velocity Ine@thesZ-direction, 1, radial 


Z 
distance and R, radius of the tube. Integrating equation (B-1) by parts and 


assuming no slip at the wall yields, 


q x du 
— — = if 2 _Z dr (Be 2) 
T 6 

dr 


Bee Use aia (B-3) 
where t represents thee and Ap is pressure drop along the length L of the 
tube. From equation (B-3) the wall shear stress, Tp» is obtained, 

= Ap : 
Tp oT R (B-4) 


Perrseevident from-equations .(B-3, 4) that, 
dt (B-5) 


substituting equation (B-5) into (B-2) yields, 


oR} 
eee ih ie (B-6) 
s Oo 


where T is defined as an average nominal shear rate, i.e. 


u 
— 
R (B=7) 
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du 
and u_ is average velocity, 7 = - rca represents true shear rate. For the 
m ie 
assumed steady flow y is a function only of t. For Newtonian fluids, 
e a - Tv a 
m F (B-8) 
and for the power-law model, 
Te emely | 3) 


where |¥| is the absolute value of Y , m and n are power-law parameters. 


To modify equation (B-6) for a porous media, it is sufficient to 


introduce U/®¢ in place of ek 1.€. 


S 
rl 


4 (B-10) 


Also R should be replaced by 2Ris where Ry is hydraulic radius and L by L! 
the length of the tortuous path of the fluid. Taking C, as the tortuousity 


factor or geometric constant, 
1 eee (hat) 


in which C, Should be obtained experimentally. Analysis of a great deal 
of data in the literature has led to the value 25/12 for packed colums 
(40, 65d). 


The modified form of equations (B-4, 6) are respectively, 
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where the subscript » indicates that the equation has been modified for 
the porous bed. 
Following Bird, et al (65C), the hydraulic radius may be expressed 


in terms of the porosity and the mean particle diameter, Ui 





D, o 
consequently, 
"$ 
2 ICES, Od alla ie ieee (15) 
> D 2 es 
p ¢ “¢ oO 
where, 
= w<AP B-1 
bia g 0G. me 
c 
$ 
w = = (B-17) 


For a Newtonian fluid by using equation (B-8), equation (B-15) becomes, 
a = i (B-18) 


Solvaneyfor ipeandssubstituting (B=15 216,717). gives, 


2a 
© aDp Som ee ap (B-19) 
: 72C, (1-4)? LU 
From Darcy's law, 
Te 3 (B-20) 
Y U L 
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Therefore, 


D > 
K = + ___, (a2) 
72C, (1-4) 


Applying Equation (B-9) to the Equation (B-15) gives, 





n 
7, = ( ae =) (B-22) 


K AP 
i. = aoe (B-23) 


where H denotes a viscosity level parameter defined be 


1-n n 
2 


H = iT (72C,, Ko) ( (B-24) 
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Equation (B-24) is identical to that presented by Christopher, et al. (40) 
Choosing C.. =) 25/12 for thesspherical beads. “It is evident that H = 1 when 
n = 1 as for a Newtonian fluid. Equation (B-23) is the modified form of 
Darcy's law for a power-law model fluid. 

et us use Ergun’s'derinition (607) on the £riction factor: 
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The Reynolds number is arbitrarily defined so that, 





= 1 x 
fe = > (B-26) 


After straight forward algebraic manipulation of the preceding Equations 


et 2-n 


ee 7 7 C-aH ee 


For a Newtonian fluid n = 1; therefore, 
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APPENDIX C 


DEBORAH NUMBER 


The dimensionless group Deborah number represents a ratio of time 


scales of the material and the flow process. It may be defined (68) as, 
Iv) 
= f 3 
oa 6 fen 
pr 


where 8 is the relaxation time of the fluid, and 6 thes duratvon of 2 


pr’ 
process which represents either the duration of a given deformation state 
or, equivalently, the reciprocal of the rate at which the deformation 
process changes. Marshall and Metzner (68) applied the convected Maxwell 
model to approximate the behavior of the viscoelastic fluid and showed that 
the appropriate measure of the process time in a packed bed of spheres was 
the reciprocal of the deformation or stretch rate in the flow direction. 


They showed that in this special case of a porous medium the Deborah number 


then had the following forn, 


—_ u _ 
Neb = e295 Be a (G7) 
p 
where, 
ape a 
u Ke (G25) 


is an average velocity in the pores. 
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The relaxation time, Be for the present study was estimated from Bueches' 
theory, Equation (V1-1). Viscosity at zero shear rate of 0.2% Separan 
AP-273 solution was taken equal to 42 poise - Refer to Fig. A-3 of reference 
(34). The limiting viscosity at zero shear rate of 0.1% and 0.05% Separan 


was estimated by using Equation VI-3 such that, 


u G . 
Loyd bi 

U a ( G (C-4) 
2,2 2 


where subscripts 1 and 2 indicate the different concentration of the 


Separan solutions. Therefore, 
He (0.13 Sep.) =e ole poise (C=53 
= A C- 
Pe (0.05% Sep.) foecy a5) 


Consequently, 6, was estimated to be, 


®F (0.2% Sep.) = 3.126 sec (C-7) 
°F (0.1% Sep.) = 0.19 sec (C-8) 
°F (0.053 Sep.) = 8.34 mili. sec (C-9) 


Calculated Deborah numbers for the three different concentrations of polymer 


solutions are presented in Tables (F-5, 6). 
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APPENDIX D 


BREAKTHROUGH RECOVERY WITH 
CONNATE WATER AND/OR RESIDUAL 


OIL SATURATION 


The procedure for calculating the breakthrough recovery with 
non-zero connate water and/or residual oil saturation is exactly the same 
as presented in Chapter II. The only modification needed to be taken into 
account is the fact that the total cross sectional area of the porous bed 
available for flow is now reduced because of the connate water or residual 
oil. Assuming that the distribution of the connate water saturation and/ 


or residual oil saturation is uniform along the bed, one can write, 
Aw a OA # oo K +e) A C=} 
where A. indicates that portion of the cross sectional area occupied by 


connate water and/or residual oil. Solving A. and ne from Equations (11-39) 


and (D-1) and applying Equation (II-34) yields, 


= oP ne LS Pepe iV 
A, ea ON re Gee 
Sw MA 
eye ; 
RS re OnEsS, (D>) 


Therefore the new.expression for fo will be, 
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Consequently, breakthrough recovery efficiency is obtained by, 
A 
FS Mer roel Z 
BR OM an (: 43.) Ory. 


This shows that the effect of the connate water and/or residual oil 
Saturation is simply a parallel downward shift of the breakthrough recovery 


Curve. 
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APPENDIX E 


PRESSURE TRANSDUCER CALIBRATION 


The calibration of transducers was accomplished through the use 
of the scheme presented in Figure (E-l). The calibration vessel was a 
cylindrical bomb which could be connected to either a mercury or a water 
manometer depending upon the desired calibration range. Valve Vv, connects 
the bomb to the positive side of transducer while the negative side was 
open to the atmosphere. Valve V, connects the bomb through a regulator 
Oran aix pressure supply. Prior to the calibration, valve Vy was closed, 
the others were opened, and the transducer zeroed. Calibration was per- 
formed by pressurizing the bomb to various levels and reading the manometer 
and corresponding signal on the recorder. 

Table (E-1) presents the calibration data of all the transducers 
used during the course of the experiments and Figure (E-2) is to obtain the 
pertinent equations of the transducers. The working equations derived 


from Figure (E-2) are as follows: 


Trans. No. 2221 with 5 psi Diaphragm: 


1B = WE 


in HO 

Trans. No. 2010 with 5 psi Diaphragm: 
Pe Hone 

Z 

Trans. No. 2012 with 1 psi Diaphraom: 
Bad H,0 =. 0 sO0L60R 

Trans. Nos 22255wath 5 °psi Diaphracm: 
a de Us 005) 

Trans. No. 2221 with 20 psi Diaphragm: 


ae Hg =) 0 01 5.8 


epeniget) vay 2 ia $105 of sehigeT 2 . 
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where, p represents pressure and R, reading on the recorder at attenuation 


1 in inches. 
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TABLE -E=1. 


PRESSURE TRANSDUCER CALIBRATION DATA 





Run Manometer Recorder 
No. Ing-On H0 Ing atx 





Trans. No. 2221 with 5 psi Diaphragm 


2 ZaS 488 
2 6.5 12035 
5 aL 225 
A LG 348 
5 30.8 610 
6 SF. 85 750 
Z, 50.6 1006 
8 67.9 1364 
9 5.6 114.4 
10 ey 125 
i 16.55 5505 
12 26.12 530 
Trans. No. 2010 with 5 psi Diaphragm 
i) Bel 136 
2 LO) BS 
3 HOSS 404 
4 2012 638 
Trans. No. 2012 with 1 psi Diaphragm 
i eee 363 
Z a. O5 604 
5 T5205 940 
4 TAA 1450 
Trans. Now 2223 with 5 psa Diaphracm (ing of He) 
1 L285 493 
2 ek 413 
5 225 741 
4 AP LZ 1100 
| 
Trans. No. 2221 with 20 psi Diaphragm (In of Hg) 
ai 2045 165 
Z ae. 247 
3 S00 384 
4 7305 Sa, 
s oe2 626 
6 DL S6 782 
7 Lies 1208 
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CALIBRATION 


4 
. 
7 
. 
. : r se iad 
’ 
‘e 
— 
i 
; . 
r af ay, 
+ <= * : ‘ "ei 1 
. ; 4 ao = 
' ¥, : 
. Nis j ¢ 
* 
a 
o~ % 
e a> cous 
‘ 
* =I § 3 
. a 3 
- 1 
«A 
Tg , 
a 
~ 
~— 
s 
i 
f 
> 
— 
4 
© 
‘ 
; ’ 
ve be 
ix 4 \ 
. i atl ® 
L 
. 





OITA’ ne ai Be ‘ \ . de 





USy 
APPENDIX F 


EXPERIMENTAL DATA 


rel Specification of the Experimental Runs 


The specifications of all the experimental runs reported in this 
Study are given in this section. The experimental rms are labelled by 
two numbers, namely Roll No. and Run No. Roll number reflects the consecu- 
tive number recorded on the original rolls of film on which the pertinent 
runs were recorded. These are the same numbers recorded in the note book 
of the raw data during the course of the experimental work. Therefore, 
reference to the raw data and also unreported data of the glass rods bed 


is directly viable without any confusion. 
TABLE F-] 
Specification of Experimental Runs 


1.  Rums being performed with the large beads, 


Bead Diameter = 4.7 mm 
Label Specification 
Roll No. 8 Water Displacing 100 cs Dow Corning 
Run No. 1 q = 1.5885 cm/sec 
Roll No. 8 Water Displacing 100 cs Dow Corning 
Run No. 2 Cues ai ese 
Roll No. 8 Water Displacing 100 c¢s Dow Corning 
Run No. 3 de ==30, 554 cm/sec 
Roll No. 8 Water Displacing 100 cs Dow Corning 
Run No. 4 q = 6.554 cm/sec 
Roll No. 8 Water Displacing 100 cs Dow Corning 


Run No. 5 q = 0.6354 er feee 
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Label 


Roll 
Run 


Roll 


No. 8 
No. .6 


No. 9 
Nos 1 


No. 


INOS. 2 


No. 9 
No. 3 


No. 
No. 


No. 9 


No. 


No. 9 
No. 6 


No. 9 
No. od 


No. 9 
No. 8 


No. 
No. 


No. 
No. 


No. 
No. 


No. 
No. 


q 


0. 
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Specification 


100 cs Dow Corning Displacing Water . 


3.4 ciibieed 


Sep. Displacing 
1.5885 cm/sec 


Sep. Displacing 


weil ci yeee 


Sep. Displacing 
6.354 cn /sec 


Sep. Displacing 
jee) Git Bee 


Sep. Displacing 
0.6354 cm/sec 


Sep. Displacing 
Ae cm” /sec 


Sep. Displacing 
6.354 Gi /sec 


Sep. Displacing 
15.885 cn”/sec 


Sep. Displacing 


62354 ae flaws 


Sep. Displacing 


Tesessten eee 


Sep. Displacing 


6.354 en reee 


sep... Displacing 


15.885 ar /ace 


100 


100 


100 


100 - 


12500 


12500 


100 


100 


100 


100 


cS 


ES 


eS 


GS 
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Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


cs Dow Corning 


ES 


eS 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 
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Roll No. 10 


Specification 


0.05% Sep. Displacing 100 .cs Dow Corning 


q= 0,554 cm/sec 


Water Displacing 100 cs 
q = 1.5885 cm/sec 


Water Displacing 100 cs 


CG ieeee 0554 cm” /sec 


Water Displacing 100 cs 
q = 0.3177 an’/sec 


Water Displacing 100 cs 


Gi) a 52277) cm/sec 


Water Displacing 100 cs 
rz 
Q = 0215885ecn /sec 


Water Displacing 100 cs 
q = 6.354 ce 


Water Displacing 100 cs 
q = 11.5885 cn” /sec 


0.05% Sep. Displacing 100 


q ‘= ss clagecmafsee 


0.05% Sep. Displacing 100 


q = 6.354 cm/sec 


0.05% Sep. Displacing 100 


q = 1.5885 Sie see 


Runs being performed with the small beads, 
Bead Diameter = 1. 5 m 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


Dow Corning 


cs Dow Corning 


cs Dow Corning 


cs Dow Corning 
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Label Specification 
Roll No. 11 0.05% Sep. Displacing 100 cs Dow Corning 
Run No. 7 Gea 152885 cen face, 
Roll No. 11 0.1% Sep. Displacing 100 cs Dow Corning 
Run No. 8 q = 6.354 cm/sec 
Roll No. 11 0.1% Sep. Displacing 100 cs Dow Corning 

v4 

Run No. 9 ah Re ann eee 
Roll No. 11 0.2% Sep. Displacing 100 cs Dow Corning 
Run No. 10 q = 3.177 am/sec 
Roll No. 12 26 cp Glycerol Displacing 100 cs Dow Corning 
Run No. 1 q = 1.5885 cm/sec 
Roll No. 12 17 cp Glycerol Displacing 100 cs Dow Corning 
Run No. 2 Ge =e sells cm/sec 
Roll No. 12 17 ‘cp Glycerol Displacing 100 cs Dow Coming 
Run No. 3 q = 3.177 cm/sec 
Roll No. 12 j 52 cp Glycerol Displacing 100 .cs Dow Corning 
Run No. 4 q =) 354177 cm/sec 
Roll No. 12 25 cp Glycerol Displacing 100 .cs Dow Corning 
Run No. 5 Ge = ie Sal/7 cm”/sec 
Roll No. 12 25 -cp. Glycerol Displacing 100 cs Dow Corning 
Run No. 6 Gael 5585 cm/sec 
Roll No. 12 “25 Cp Glycerol Displacing 100. cs -Dow Coming 
Run No. 7 q = 6.354 cm/sec 
Roll No. 12 28.5 cp Glycerol Displacing 100 cs Dow Corning 


Run No. 8 q =. 0.6354 cm/sec. 
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Specification 


0.2% Pusher 500 Displacing 1000 c¢s Dow Corning 
0.15885 a ieee 


rQ 
iH] 


0.2% Pusher 500 Displacing 1000 cs Dow Corning 
0.6354 piinieee 


Q 
" 


0.1% Sep. Displacing 1000 cs Dow Corning 
WEE Gace 


0.1% Sep. Displacing 1000 cs Dow Corning 
7a 
q = 0.6354 cm/sec 
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Pee Frontal Plane Advance 


The procedure undertaken to obtain the data points of Figures 
(V-3) to (V-18) is described in the following. The position, Xo» of the 
frontal plane or the plane of zero saturation was determined by projecting 
the movie pictures frame by frame onto a white wall. The distance of the 
projector from the wall was set so that a magnification of 1:1 with 
respect to the actual dimensions of the bed was obtained. The mumber of 
the frames was counted to calculate the time elapsed, t, by using the 
previously determined camera calibration factor. The time of breakthrough 
was checked by a stop watch which determined the time period between when 
the displacing fluid passed the inlet screen and breakthrough. The distance 
from the screen to the tip of the longest finger, Xo» Was measured from the 
projections bya precision steel ruler with lmm divisions. In the case of 
low flow rates to reduce film consumption, the camera was stoped after taking 
a few pictures then restarted again. The time lapse between every two sets 
of pictures was measured by a stop watch and added to the time calculated 


from counting the number of the frames. 
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TABLE F-2 


POSITION OF THE ZERO SATURATION PLANE 


X, t at/oA 
cm sec cm 
Roll=No. 8, Run No. 1 
S00 9352 0.66 
Soh 19.05 5 1. 
8.0 28257 1.97 
JPRS 38.10 2205 
1S yee 47.62 3weS 
Were Sia. 3.94. 
Loo 7 66.67 4.60 
7a ely) 76.19 ite FAS 
245 85778 5.91 
Zos5: 95.24 6.57 
552.0 104.76 [=25 
SoD 114.29 7.88 
39.0 123281 8.54 
43.0 T3557 95 9720 
47.0 142.86 9.85 
5025 152258 LOST 
Ho.) 155574 10.74 
Roll No. 8, Run No. 2 
4.5 Vela 0.99 
8.0 14,29 Voy 
EDs, Cle AS 2.96 
18.0 23257 5.94 
EAD 55 4.93 
LS. | 42.86 Do Suh 
28.0 50.00 6.90 
ew) See 7.88 
36.4 64,29 8.87 
44,3 71.43 9.95 
48.0 184.57 10.84 
Seles, 91043 P2500 
Roll=Noz 785) Rum «Nos 3 
Del 4.76 deseo 
Lass Sis 4 2205 
1Be5 14.29 3.94 
18.0 19.05 Dis LD 
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qt/ oA 


8A Liod 


S-4 ay 
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fi 






X t qt/oA 

cm sec cm 
Roll eNow &; Run No. 6 

4.0 25.61 5250 

125 VheOZ AOS 
ie [MES 10.54 
TANS SAL T4206 
17 DLEOSO5 Lha37 
20.3 142.86 Zae 0S 
PAP ES 190.48 SEN 
5502 258210 Rhoral Ls 
43.0 Yas toe 8 ALS 
5128 COE ae) aoe 2. 
55.40 346.19 5b. 0 
RekleNo, 9. ¢Rum=No. 1 

Zr8 25 ok 1.64 
a20 Miao 3.28 
6.0 TEAS 4.93 
920 95.24 OF5a 
9.5 112.86 Fel 
Lies 297.96 Z0nao 
VAS ES 529.86 LL eS 
Soo 429.06 29259 
S80 456.20 31.46 
42.0 550560 36.59 
AS. 571.08 59559 
5200 652.78 43.64 
Doeo 652485 44.95 
RebeeNO. 95 skuneNoO. Z 

4.5 25785 3.28 

8.0 46.62 G25 4 
1245 He, 45 9 85 
LGs5 94.76 13507 
32a) 194.96 26.89 
5050 219. 25 50.24 
5859 228. teh Sloss 
5b 0) 294.57 40.63 
Be) 505..52 AQ. ee 
RO JING we 9D). ser NO. ES 

Za) 4.76 Dees 


TABLE F-2 (Continued) 
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Xo ic qt/oA 
cm SEC. cm 
Roti Now 95° Run Nov 3 
AS. 5 338 53 22.99 
49.5 55¢ 7h: ZA%.03 
5a 380 10 26227 

RO1 Novae Ohm Rum 'Nocwe: 
Cee 4.76 bere | 
6.9 Oey 2 OS 
a. 14.29 3.94 
16.1 OOS Dao 
SES Zoe On 6n.57. 
fags Zon Oo 7.88 
Zoe i ASS) 9.2 
Zoe EL 38.09 10352 
550 42.86 11382 
She MeOe Ieee: 
Sced 5758 Tae 45 
42.8 oN esa 1 S275 
AT.2 61.90 17208 
Se 66.67 Heese 9) 
Doe 70.00 TORS 
Roll No. 103. Run No. 3 
Ane dh Wve teh L.64 
AS A Shon dein) 
14.9 aod, Ok 
17.6 LEST ales) 9.45 
24.6 (AES STS 14.86 
TAS 239724 16.50 
3185 247 eon L7. 09 
5720 500 1a L007 
Ae 6 52492 TEA TAM, 
43.0 5a oe 7m 24505 
44.0 5 Sieae 24567 
49.8 404.08 Lie OW 
5525 449,79 351802 
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12.9 LOS 0.66 
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TABLE F-2 (Continued) 


x t qt/oA x t qt/oA 


ie) 
Oo 


RolioNo. LOS > Run No. 4 Ro lis Now 05 Run) Not! 6 
54 Te XA L..97 ah a 38.09 LOaS 
8 11.90 528 29.4 42.86 I sS2 

ha ToOn67. 4.60 55 oa. 47.62 oes 
pS Alt AS 5.91 36.8 See 14.45 
7 5 26.19 Ls 22 AG SS Gi A iS S26 
20.8 30.95 8.54 47.0 61.90 IW 208 
aa Boel 9.85 BaD HOAs 19.44 

2655 40.48 Lely 

BO05 45.24 12.48 Rodi) Nowell Rin Nowa, 

oes ate Oe 5.0 jaa tase 

40.5 54.70 Sr alal. 
9.0 14.29 1.06 

4255 59.52 16,42 
eS Fil Wig ES ler) 

AS 30 64.29 LATS 
16.0 Zoe 51 Dot 

50.0 69.05 19.05 
53.5 74.76 20.62 19.5 SiS il ale 2.64 
; VS ea | 42.86 oral 
= De, 5020 3.69 

Roll No. 10; Run ‘No: 5 98.5 57.14 4.22 

4.5 2.358 Pood S726 64.29 Les 
10.6 4.76 S78 S47 (aR) Des 
PASS Ps LA 4.93 59 S35 Se Fook 
735 9.52 Gas A4.0 aes al G255 
zea) 11.90 Sk A7.5 92.86 6.86 
a24 50) 14.29 9.85 SEO TOGO 7.39 
ZBL 16.67 11.49 BS 55 LOS 7 Horone 
Sfaul 19.05 13.14 
Steps LW. 435 14.78 Rom Noe 05 Run No. 
oe eo | Bub one 23.81 0.70 
AG £5 26.19 18.06 

Low ATNO2 1.41 

50.8 2S yh 19.70 9 
63.5 30.48 71.02 ESO 7245 FEL 
: Nie) 88.57 ZAGZ 
Dovel LOS E47 4.89 

ROWE No. eu) Run No. <6 16.0 172.14 5.09 
4.0 AO 5d: 28.5 194.04 Sip hes: 
6.9 9.52 IOS 39.5 25 Oe 7.43 

LOS! 14.29 3.94 Ses Lio Selo 
14.0 19.05 Sas SR 294.2 Su720 
LO56 25.8 6. 5k 
EOL ZS, Sel 7.88 
2520 eB 9.2 
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Xo t qt/oA 
cm. See cm 
ROLIENow? 10%" Rum’ No. 9 
oe, Day onl 0.35 
4.8 47.02 0.70 
6.0 74.29 110 
Bes) 141.19 2.09 
14.4 165.00 ape i 
Z4*0 ZAG C20 3.64 
Te 7 2097.05 3.98 
30.0 344.65 5.09 
oe U 509.41 5.46 
oy ieee 452.01 6.68 
AZ 474.87 Pea? 
A7.5 Doe OF ical lay 
5055 Sy MOS 8250 
5270 O50) 5 9.65 
DLO 672228 oS 
oe 752.98 USS Ss: 
ROLL No. 105, Run No. 20 
DU) 4.76 0270 
Sars Siapey4 hese): 
14.0 14.29 Force 
E750 19°05 2.81 
Vi vues ZOO. See 
24.0 Zouot AS eee 
28.0 ieee) oS 
Son 38.09 5.65 
500 42.86 OF 55 
A2.5 AO 7.04 
A3.5 Si 8 the 
48.5 Sioa tt 8.44 
50.0 61.90 Shee kes 
sere) 66.67 9.85 
ROLL Noe. lie shuns Nowe 
526 96.67 OSL 
4.6 160297 1.19 
oa5 LSZ. OH. A rihas 
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Xo t qt/oA 
cm sec cm. 
Roll Ne. dl, ~ Rum sNo 5 
a5 ae Gs ae ge 
8.0 O..52 thea 
IDS 14.29 AneZ 
16.6 19.05 B05 
20.5 Sent 7204 
BS 28)..51 8.44 
Z29R 0 oi) ees) 9.85 
Bye 5 Soe 0 0 
SonU 42.86 L207 
ET A Oe 14.07 
47.0 Bo MS) 15.48 
49.5 Ses 16.89 
See) 61.90 TSU 
ROL Now Ly Rum No.6 
Sua ESM 6 
8.0 47.62 5252 
10 G2n, 38 Aol 
19.0 HL Gy aes: 9.74 
2785 156554 157 
34.0 223.04 16.48 
ieee 246.38 187,21 
MS) 286.28 Zils 
AS. 7 503. 42 22a? 
5.5 Sie ie 24.07 
Roi eNo. ly Run Nos 7 
GR Diao) 1.70 
AAS) 4.76 5 oe 
15.50 Tele BeZs 
Ae) SG) On 57 7.04 
2525 1290 8.80 
2O.5 14.29 LOY 56 
50n:0 Los 67 eZ 
Ane 0 {Oe 05 14.07 
48.0 Lie AS 15563 
50.0 Zon Ol 1 7e59 
bo Sie) ie ae 18365 
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Peo Pressure Drop Measurement of Single Phase Flow 


These measurements were performed to obtain the permeability of 
the bed and also to establish the friction factor - Reynolds number relation- 
ship. Pressure drop was measured either for the total length of the bed or 
for the top or bottom half section - The small portion of the bed before 
the first and after the last pressure tap is of course excluded. Close 
agreement between the absolute permeabilities calculated for each section 
was an indication of uniform packing. The absolute permeability, K, was 
calculated from Darcy's law, Equation (B-20). Friction factor and 
Reynolds number were calculated from Equations (B-25, 27) respectively. 

Single phase flow data are recorded in Tables (F-3) to (F-6). 
Tables (F-3, 4) are for Newtonian fluids while Tables (F-5, 6) present the 
data points of the polymer solutions. Deborah number recorded in Tables 
(F-5, 6) was calculated from Equation (C-2). Experimental fluids involved 


in these measurements were discussed in Section (IV.2). 
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F.4___ Pressure Drop Measurement of the Displacement Tests. 

It was mentioned in Section (IV-6) that the pressure drop history 
of each displacement test was recorded on a Strip Chart recorder. The 
pressure drop recorded on the chart, AP», Was the total pressure drop which 
consisted of pressure drop of the flooded zone, APE plus pressure drop of 
the o1lizone, AP The value of AP, for the single phase region could be 
calculated by knowing the length of the bed above the plane of zero saturation. 
This was obtained as described in the following paragraph. 

The time when the displacing fluid passed the inlet screen was 
marked on the pressure drop chart. Therefore, by knowing the speed of the 
chart, the corresponding time, t, of eachypoint on the chart could be 
obtained. This value of t was used to calculate qt/¢A. From the eye-fit 
Fane of Fipures ‘(V-3) to -(V-18)»the value of X, or length of ficoded zone 
was determined. Therefore the length of the unswept portion of the bed could 
be easily calculated by difference and subsequently used to determine AP: 

Table (F-7) includes the pressure drop measurements of the flooded 
zone and the overall pressure gradient of the flooded area, o ; 

Volumetric average water saturation, Ce which has eee been 
defined in Section II.3 was calculated at each position in the bed. This is 
“included in column four of the Table (E=7)...- Correction of the pressure drop 
reading due to the small difference in specific gavity of the displacing and 
displaced fluid has also been considered. Occurrence of the breakthrough is 
marked by B. T. which appeared in Columns 3 and 4 of Table (F-7). 

Table F-8 presents the friction factor-Reynolds number relationship 
of the flooded zone. Friction factor and Reynolds number were calculated 


from Equations (B-25, 28) respectively. The only modification involved was 
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replacement of the viscosity in the Reynolds number by effective 
viscosity defined by Equation (II-33). For the polymerflood, Equation 
(II-33) with apparent viscosity of the polymer solution instead of Me 
was used to obtain the average viscosity for the Reynolds number 
calculation. Estimation of the apparent viscosity of the polymer 
solution was based on the shear rate in the bed. The reason behind this 
choice is discussed in Appendix HH. 

Tomcaiculate. the.iriction.factor ofthe flooded zone the arith- 
metic average of the overall pressure gradients with astrisk in Table (F-7) 
were used. The initial data points were ignored to minimize the error 
contribution due to the initial segment of the bed which is discussed in 
Appendix G.1. 

For the sake of comparison, in addition to the experimental values 
of the overall pressure gradient of the flooded zone its predicted values 
based on the effective viscosity and Darcy's law are also reported in 


Table (F-8). 
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TABLE F-7 


OVERALL PRESSURE GRADIENT OF THE FLOODED 
ZONE 


* The arithmetic average of astrisked overall gradients is being used 


to calculate friction factor of the mixed zone. 





a a RS 





t qt/oA Xo S, AP er APE, (dP/dx) ,, 
sec cm am. cm of 110 cm of H,0 cm of H0/ 
cm 
Roll No. 8, Run No. 1 
0.0 0.0 0.0 0.0 ee - = 
30 Ze 07 One) On2t 6.4 0203 0.004 
60 4.14 18.8 0.22 5.8 0.16 0.011% 
90 Gar 28 0.22 BG 1€12 U2045% 
120 8.28 Dat Oe22 50) 05 0.046% 
156 10.34 46.7 Ove2 4.4 eo A 0.044% 
180 jeagul 5525 Oe 5 Ae Ze2 0.047% 
210 14.5 ale Bale 3.8 18 OO 57 
240 Te. 5 2 = eA Thott 0.029 
300 2057 = - Sed ee 0.024 
360 24.8 - ~ one 1.0 G202 
420 29 = “= 2.9 0.9 0.018 
480 55.1 - = Mots 0.8 0.016 
Roll No. 8, RunvNow 2 
0.0 0.0 0.0 0.0 1625 = 3 
‘SS Ae 5 19 0.24 Se Tae Onr5* 
48 6.62 28 0.24 ae Ge 85) OL 
78 10.76 A Ss 0224 1083 Ons G2i5% 
108 14.9 Bole Belz oil 6. 74 Glas 
138 19.03 = = ee ie Oe ae 
168 (eS = - 6.4 4,44 0.09 
198 2.351 - - Sac 3.84 0.08 
228 38.45 2 = 528 3.84 0.08 
258 55659 = - ge 3.34 0.07 
; Roll No. 8, Rum No. 3 
0.0 0.0 O20 0.0 52.95. ae = = 
15 4.14 14 0.30 30815 SOS 0.47% 
30 8.28 292 0.28 25.95 11.74 OQ ASt 
45 1741 AZeS 0.29 C554 5 On 0.42% 
60 16.55 Bob Bau 2055 18.59 Ueo0s 
90 24.83 = - 7 Loss0 os 0.356 
120 SanlU = - L515 PSL Wes 
150 41.38 = - 8.8 6.84 0.14 
180 49.66 = = 1s) 5.59 Heal 
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TABLE F-7 (Continued) 


t qt/oA i Sy AP AP (dP/dx}, 
cm cm. of H,0 cm of H,0 Cneor H,0/cm 
_ sia ae a eae ae ee eee tee 


Roll No. 8, Run No. 4 


0.0 0.0 0.0 0.0 Lesa te) 7 7 

iS 4.14 14 0. 30 26.42 3.45 Oras 

30 8.28 292 0.28 Fame Seine 0.34% 

45 Type a AZeS G29 olathe! WSS 0. 38* 

60 1lo.55 But Bol. 20509 Te. 71 On 30" 

90 24.83 = = 19.42 17.46 C3507 
120 ao! LO = IG./5Z 14.56 0.30 
ie 37 24 = = Is Bie) Gig 0.20 
150 41.38 = = Sa55 Gao O.13 
180 49.66 z - 7.04 5.08 0.10 
210 i 2 = igo 4.73 OG 
240 66.21 = = 6.02 4.06 0.08 
270 74.48 - ‘ SOS Boo 0.08 


Rol Noy*s,8 Rum No.5 


0.0 0.0 0.0 0.0 2.94 = : 

Vin 0358 ae O23 B50 0.43 0.86 

a 1.41 De! ORZ5 era) UE SS) 0.09 

81 PETERS 930 Q225 Op tks) Omri 0.05 
im 3.06 P25 0.24 Bas 0.44 0.04% 
141 S09 a7 OZ 5 515 U6 0.48 Qs04% 
yp 4.72 dig, 0325 510 0.48 O705* 
201 5.54 22a 0325 Becks 0.41 2 
216 5.96 24.2 U25 ws95 0.42 0.02% 
Z51 O35) 20.0 O25 2.84 OC35 0.014 
246 6.79 PoE 0°25 2385 Cer 0.01% 
261 7 a20 29 20 OF25 2.80 Orav Odes 
Bo 8.03 5235 Qa25 TA ei) Os05 Ds QaG7% 
SA 8.86 36.0 UAGAS 2.84 0.56 0, 016% 
381 10351 Uae O25 Zo k 0.76 0. 019: 
44] A 49.0 325 2a 1 Osyt a0 TS 
501 135.82 B.T Bal 2.42 0.46 0.009 
561 15.48 : = 2.45 0.49 0.010 
621 Laks = = Za50 O59 O5.012 
681 13579 . = 2011. 0.08 0.002 
741 20.44 = _ BMS 0.36 0.007 
801 ZLabO = = RACY, 0.02 0.000 
861 Gogo 7 = Deep Z 0.16 0.003 
921 25.41 = = 2.34 O56 0.008 
981 Zi = = Tie hs, 0.22 0.004 
1041 ZO. z, % = ACIS Oley, 0.003 
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TABLE F-7 (Continued) 


cm of H0 


Roll No. 9, Run No. 1 


USO 1.2o 
0.91 UES) 
0.86 OO 
0.85 ae 08 
0.84 F208 
0.84 14.88 
0.84 16.2 
0.85 18.64 
0.84 50 
0.84 24.36 
0.84 201.96 
0.84 29350 
0.84 G2)./50 
0.84 Sree) 
0.84 50276 
Bel 56.26 
36.26 


Roll Nos 9, Run No. 2 


0.0 14.35 
0.80 18.6 
Sica) 22605 
0.85 onal 
0.84 20.50 
0.84 S1s5 
0.84 Rory’ 
0.84 57 38 
0.84 Al.45 
0.84 44.60 
0.84 46.60 
Bal, 46.25 
= 46.84 


AP (dP/dx) » 
an of H,0 cmvot H,0/cm 





0.02 0.0004 
0.02 0.0004 
G02 0.0004 
0.0 0.0 
0.0 0.0 
AS 1.24 
4.36 Oeo5 
O02 0.93 
8.07 Ooi 
LOES6 0.86* 
Vaee2 0.79% 
16.48 One 
20.00 0.75" 
Zoe LS OR ee 
20427 0.2" 
Z0m9 WR Slog: 
55208 Ooi 3° 
34.80 Ones 
34.80 0.70 
34.30 0.70 
D0: 1.58 
9.64 lye) 
13.94 ie ag 
18.4 105" 
22.84 1.025 
Pi thsitss Ls01* 
51.67 Os oy 
50855 0296" 
40.90 Oso7* 
44.13 0, 94% 
44.29 0.90 
44.88 0.92 
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TABLE F-7 (Continued) 
t qt/ oA x ci AP Ap (dP/dx) y 
cm W a F F 


cm cm. of H0 cneOt H.0 cm. of H,0/cm 


OAS SS ee SNe Ra gee ee em eee ee ee eta 


Roll No. 9, Run No. 3 


0.0 0.0 O06 0.0 Reveal 5 * 

24 6.62 Sao 0.78 3740) 12.68 1.95 
54 14.90 18.8 0.79 43.80 24.60 1. 46= 
84 ye ee yj LOS 0.79 Syeone 0) 59 CIS 14a 
114 5145 59.6 0.79 58.60 SoA Se 
138 38.07 48.1 0.79 64.60 61.09 Lass 
156 43.04 55.65 0.80 64.80 62.84 Ve2os 
E74 48.00 pag tee Bet 63.40 61.44 125 
204 56.28 = = 65:40 Gia Hi 9(6) 

Rodd sNo..*9, = RuneNos 4 

0.0 0.0 0.0 0.0 Hsee! = 

9 Gaal: Sis ORE SoZ 22.84 oad. 
15 LOZ55 bees O..72 OA 37554 5.05% 
raat 14.48 ne 2 O72 99.7 51.58 2.054 
ZT 18.62 Tl Sete) On7Z 109.2 69.01 Loge 
a Zend 8 Sine Oe dias ee Sin Oi! Payee 
39 26.90 Sioa Wey 2 iors 94.45 Zou 
4S eee 43.0 OF72 19-10 LO S765 2.55" 
51 55. 47 48.7 O-iZ (LALSe P16.79 oA 
57 BOS a Belts B.T WN DUS 74 Deon 
63 43.45 = = V5 27 RO ED Bese! 2.28 
69 47.59 * = 114.0 112504 2.29 
75 5 La/S = 3 EEA 1.64 Lote 
105 iowa d - = aS tar ie 2.28 

RoliyNo] 9. Run No. 5 

0.0 0.0 0.0 0.0 LIS ace - ee 

30 0.83 Tad Ost 294.76 we) 5 8.65 
Ag oe iy ORL 28/210 64.23 6.02 
78 ZS 19.0 Orda BMG ES OO me Se oe 
108 2.98 26.5 Ou 275.89 135.30 5.54% 
138 Cate! Soe Oxd 1. 269.07 169357 Ook 
168 A 63 Ade Oech Zal <0 200.09 Sq ltt be 
198 5.46 48.7 Oak 248,12 Deyo 1S) A. 99% 
228 O22 9e. Dee Bul 25031 254.16 4,78 
258 dak? = = 251.50 PARSE SY 4.68 
318 Si. = = 2178 209.380 4.28 
378 10.43 = = 1ST 04 189.08 3.86 
438 12.08 : = 179.86 177.90 5.65 
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TABLE. F-7 (Continued) 


x S AP AP (dP/dx) 


cm cm. of H,0 cm. of H,0 cm of Hm 


r a 166.32 164.36 OB, 21a) 
‘ > 148.78 146.82 6400 
: = 1297524 LP lee 8 2760 
4 = 120.06 JEL Seal Boe. 
; (DTN JRL JOZ Za28 
ss zi AO oy 100.56 2.05 
# 5 O76 S104 1.84 
: "4 88.16 86.2 16 
. . 718.98 17.02 Masi 
5 2 66.62 64.66 ee 
. F 57.44 55.48 Leal 
Roll No. 9, Run No. 6 
0.0 0.0 Le sO - = 
5.0 Ons 76.07 4.01 AYO 
6.0 Onl 72.94 Pay 2.88 
120 Oni4 166.16 26.08 PRA 
18.0 0.14 160.64 41.74 LeOl* 
24.0 0.14 eae! 50:10 Desa 
5052 0.14 ee a 2 Ou 2.505 
SOs. 5 0.14 13%. 94 84.63 ae 
A2.5 0.14 eZ O;. 79 OA aT) Dear 
48.5 Os 4 ie Oo 108.79 ie 
Bole le ie 110.84 108.88 Dcnlad 
= = 100.48 98.52 ON 
= z 94.30 SVAe oul 1. 88 
= - 82.48 S0Na2 1. 64 
= - 14.595 Laat 1.49 
= = 68.48 66552 le Se 
= > 65.40 63.44 LaZ9 
= = O72 O5..16 SS 
= = 50,20 54. 30 ee sh 
= = Seeks a0 a9 IeuZ 
= = 48.58 46.62 (Oe) 
= ; AS. 94 41.98 0. 86 
Rol No. 9’ Run No. 7 
O20 0.0 ZoLoo = - 
LOe5 0255 2O.29 4.58 Oa54* 
1825 0.54 20249 8.68 Qe55° 
26.0 0.54 ZOnoo 7.0 455" 
Lee Oras 26.69 TOO Oeo1 
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f qt/oA x S AP. AP. (dP/dx) r 
Sec cm an W cm of H,0 cm of H0 C.on H, 04m 
81 TREN 40.8 Onoo 25.96 19.03 0, 49* 
96 26.48 48.6 OF 54 LD) 22.19 0.46" 
Hee 30.62 Bal. iSegl ie. 24.19 21S HOD ey = 
126 34...76 ~ - 25 Lt AOD 0.44 
141 38.90 - - Ae Clas) 0.44 
17i 47.17 = - 22.89 20.293 0.43 
201 55.45 = - 22.89 20.93 0.43 
Roll No. 9, Run No. 8 
0. 0.0 0.0 0.0 65.05 4 = 
12 S28 A 56 0.55 66,57 19 503 jae 
24 16.55 30.2 0.55 65, 81 Sr Ou 131s 
36 24.83 4585 O55 64.92 55.88 1g23* 
48 55. 20 Bal Bale 57 8L 55.05 Las 
66 BS ee 2 - : AG. 71 44.75 O21 
84 oS - - 45.40 43.44 0.89 
102 70.35 = = a8 eee 0.89 
Roll Wo 105 Rum No. 1 
ae 0.0 (bei) 0.0 7a 8 = : 
24 6.62 1G 0.38 23.0 5.48 Cis 
oy 14.90 39.6 Omo6 Z050 12073 O. 54% 
84 aaa cialis Bs le re 0 13704 0.2 
114 5E.45 = = 13.75 LAeeg eet 
144 58 ae = = 13.0 11.04 G25 
162 44.69 = 2 as 10.54 OR ce 
RO] No. 10, Run No. 3 
On 0.0 0.0 0.0 One = = 
a6 Lapggass 4.0 Oey On ORY, 0.29 
63 454 ESO Qeesy! 6.6 ORSS Oee es: 
93 6. 41 jhe! 0.56 Ce Nhe) 0.15% 
R25 8.48 [Sak R56 6.8 diapewe Oapilae 
183 2 bZ ii 0256 6.9 POSS Op, ee 
243 16% 6 30.0 0.56 ao atl Oe 
303 20 00%: ye 0.56 7.06 eS) (ne 
OS ZS AAS) 0:56 ele Mase Diepu lige 
423 Zoe Bare 0.56 6.9 4.94 Og Os 
483 555 SL Bale Bole GEG 4.64 0.09 
543 5h AS = = 6.0 4.04 0.08 
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t qat/oA x S APp APY. (dP/dx) F 
sec cm cm cm of H,0 cm. of H,0 cm of H,, 04m 
Roll Now 10, Rime Now 4 
0.0 0.0 0 O20 24.9 = = 
15 4.14 HES) 0.38 24 of5 4.06 0.45% 
30 8.28 LLiwh Wea 24.0 S705 Oaa2* 
45 Lag4y 556k On 37 DS as 132k6 0.42% 
60 14.55 44.0 0.38 Faas Lao Odi 
1S 20.69 Bot Berk Z20inS 18.54 Gis SS 
90 24.83 - - 19.0 17.04 O25 
105 23.97 = - G5 14.54 0.30 
120 53.10 - = US aes: LSe79 0.28 
155 5724 = = 15.5 13454 Ou28 
150 41.38 - - YS.25 13.29 OZ? 
180 49.66 - = TAs 12.79 Ou26 
204 56.28 - = 14.9 12.94 W.26 
ROLL No. 210) kim No. 5 
0.0 0.0 0.0 0.0 Sui) = = 
pee one a 14.7 os 55.0 9.30 Wea 
15 hOsS5 Lhe 0.38 AS.0 14.49 OVSe* 
20S ee AD. 7 0.38 50...0 15.035 0.40% 
30 20.69 Bele Ba lee 25.0 25.04 0)..47 
i 2 86 = = Zia 25.04 0. 4:7 
2.5 56.21 - ~ 25.0 253 04 0.47 
Rol Now 10" Run No. 6 
0.0 0.0 0.0 One 24.7 - = 
5 A, la 10.8 0.38 2A 4,09 0.46% 
30 8.28 Dea Oey Teas 8.17 0.40% 
45 D2 ad: 5300 0.38 ORS 12.49 0.39% 
60 lass 44.0 0.38 OS 16.04 Oy 3582 
75 20.69 Bit Bebe Oe 5 Liebe 02 56% 
90 24.83 = - fe 6) 15.04 05. Sal 
105 28.97 - - Se eo 13.49 0.28 
20 5540 - - 1D. 0 13.04 0. 279. 
155 owl eal! = - a0 13.04 Oye 
150 41.38 = - 15 12 4 0.26 
165 AST 52 - - Lae5 IZ 54, 0.26 
195 53.79 = - 14e25 12.29 Oe25 
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TABLE F-7 (Continued) 





t qt/oA x = Ape AP, (dP/dx) 
sec cm. cn M cm.of H,0 cm of H,0 am. of H, 04m 
Roll No. 10, Run No. 7 

0.0 0.0 0.0 ee) 52.48 7 7 

Zt tea oa0) Dike 49.98 Ut sed 0.59 
36 2206 Loe5 On 16 AS Oaugll 0.60% 
66 4.88 30.0 0.16 A222 A ORS 0.028 
96 7.09 Ae 0.16 36.28 26413 026052 
126 Ole Sh Bae B.T 29.17 CAN 0.56" 
156 e535 = - 24.98 25.02 0.47 
216 15296 = = 20a Les0Z 0.38 
276 20.39 : - 18375 16.79 0.34 
336 ZARBS - = 16.98 1 02 Oe sk 
456 5501/0 = 2 15.46 13.50 0.28 
576 42256 : - 14.18 V2 a22 Oe2a 
696 5g 4s : 2 ban 4 7, lip be sy 0223 
816 60.30 = = T2495 1037 Olez 
936 69.16 - = 12.09 CaaS 0.21 

Roll Ne, 10," Run’ No. s 

0.0 0.0 0.0 0.0 22400 : 

30 0.89 Sore, Un iG LZbad ® NS 0.44 
60 ie 7h 1025 Cie ly 21255 Lvehe On. 55s 
90 2.66 16.0 Onli 20 par? A Oe Sas 
120 5a 55 2125 Op ley 19.66 5 yo OF 272 
150 4,43 ZOn0 Orsi le? 18255 6.62 0. 23" 
180 be o2 oe Ore 17.44 Og 0.25" 
210 Oe ZA. oy Oaks i Cree be 8.00 Once 
240 7.09 Aan Oey 8.5 9.78 Ose eie 
270 7.98 A Tae Chea 14.84 iLL el QaZoe 
300 8.87 Bam OeL7 14g 3 Lee Ce 25m 
330 9375 Bee Dek MSS 2S) Val a On 25 
360 10.64 - : 12.95 10.99 0.22 
420 Peas = = 1S) 10268 OG A 
480 14.19 = : 11.98 10.02 0520 
540 15.96 = : 11.86 9.90 0.20 
660 19.51 Z = 1/6 9.80 0.20 
720 Zin 28 0.18 
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t qt/oA x 5 AP. AP. (dP/dx),, 
sec cm cm y cm of H,0 cm of H,0 Cmso£ H, Ofm. 
Roll No. 1057) Rum No. 20 
0.0 0.0 0.0 0.0 2). 67 - : 
18 2.66 16.0 aly: TOS a Za 34) 1,615 
53 4.88 29.0 Quel? 91.39 50n5, yg 
48 7.09 as 7 QaL7 Val 22 VAR (oes! 1.35" 
63 O51 Bal. Bet. 64.44 62.48 aoe 
78 11.55 = 48.78 46.82 0.96 
93 roy74 = 38.20 36.24 0.74 
108 15.96 2 BAST. Seek 0.66 
138 20.39 = 29 £5 Oe: 0.57 
198 29.26 - 25.04 25.58 0.48 
258 58415 = ZLeh9 20.83 0.43 
318 47.00 = 20950 Te. 34 0.37 
378 55.86 - 18.62 16.66 0.34 
438 64.73 - 6655 159 0.30 
498 £5.60 - Reece £5.95 0.28 
618 O1 255 - 14.00 12.04 O25 
738 109.07 - aS TOS Wy Qa2Z 
Rodd No. SEIS Rum No.2 
0.0 0.0 ee) 0.0 250.02 = = 
9.0 2.66 125 0.20 252.25 Aa 1d 3.84* 
to.5 4.88 24.5 oe Ae) 228.80 67.96 So. 02% 
24.0 FAS) See) 0.20 ZO eS 104.79 S013" 
515 O a5) 46.5 0.20 156613 126.:65 Z2.85* 
39.0 ao Bela B.T. 112,09 TDG; dS Dee 
46.5 Lows4 - 84.81 82.85 169 
54.0 15.96 - 66.92 64.96 ree 
61.5 18.18 - 60.53 5S eos 1220 
70.5 22a - 54214 Sais 1.06 
106.5 31.48 - Ae AS OLA 0.92 
PSOE) 40.35 - Aa 27 AZ SL 0.86 
166.5 49221 = 41.92 39.96 Blake? 
Rodi Noe diy Run No. 4 
0.0 OO 0.0 0.0 126.43 - = 
30 4.43 1128 0258 VS a4 S 21.90 Lee 
60 8.87 Zaew (es 8 a43 Ae Use hes 
90 13450 5525 Oa57 IS a2 7 71.89 2.15% 
120 1S AT 2 (0738 106.34 94.73 Za 10* 
150 UIA G ee lh pele 100.93 98.97 2502 
180 26.60 = 96.74 94.78 193 
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Roll No. 11, . Run No. 
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cm of H,0 cmsof H 
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pee! 
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0 cnsor H, O£m 
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£ qt/oA x > AP APr. (dP/dx) r 
sec cm cm Y cm of H,0 cm. of H,0 (wil (ose H,04n 
ROLIDNoO. Li aRun No.8 
0.0 0.0 0.0 0.0 224229 E = 
iS 4.43 Sil Ge55 Z2he92 Biles 25 SZ 
30 8.87 T6=3 0. 54 PAS ee SESS Bes 
45 13230 LOS 0.55 245.01 121.90 Wie 
60 VES Cans OF 55 247.09 161.19 Seo 
90 26.60 vey) G55 2S Oe 1 BENS OM Ba IU Oky 
120 Sete Bal Belg 247.09 TIS AMES 5.00 
150 44.34 = 247.09 ESS Aa ES 5.00 
180 5 5e20 = 242.94 240.98 4.92 
20 64207 = 241.91 239.95 4.90 
Rois Wo. Lis. Run No. 9 
0.0 0.0 0.0 0.0 TT 2p5 2 - = 
30 4,43 6.4 0.69 U1 Opa WCRI 5.0 
60 Saa7 jecee 0.68 D2Am Sz 37.00 eee 
90 5.2050) Oy, 0.68 128. 46 56.00 3. £o* 
120 View 26.5 Oelew! 50054 fe, 40 5. One 
150 ZB 35.0 0.67 133.04 91.14 2.94% 
180 26.60 39.8 O67 Sieh 110.09 ZOU 
210 S104 46.5 0.67 140.90 Zee o0 2.89* 
240 55.47 ae 0.67 140.90 138.94 204 
270 39.90 Danie Balke 144.00 142.04 2.90 
330 AS a7 - 145.04 143.08 2292 
390 57.64 = T4504 WAS a8 2.92 
450 66.50 = 145.04 143.08 2.92 
BLO Tonk = 145.04 143.08 2a Oe 
540 79.80 = 145.04 143.08 Ze 02 
RoaNor il Run wNoe 10 
BO. 0.0 0.0 0.0 el 39 = = 
15 Pee Z.6 0285 16455 6.50 10.85 
4S 6.65 Sal OV82 144.00 46.24 ANS 
75 11.08 UBS O82 170.94 S525 (ee ee 
105 15 52 LOWG 0279 195.80 12S ae 1203* 
Iie 19.95 24.6 0.81 UES Sie: Sa As CoO 
165 24 58 50raL 0.81 Zou 181.88 Oa * 
195 Z Om ee 55.46 0.81 247.60 Za: 25 eke 
225 i ae. 41.0 0.81 LOS=52 244.03 On Oe 
255 37.69 46.5 0.81 290.08 LIS 7 6425" 
285 A ee led Sule 0.81 515.059 Seale Geoo 
315 46.55 Bele Bele 326.34 BW Si3} 6.62 
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TABLE F-7 (Continued) 





ANE 
x 4 AP. APs. (dP/dx),. 
cm cm of H,0 cm of H,0 cm of H, 0fm 
: = jo4 41 SOLAS 6.78 

Si as = 340.84 338.88 6.92 

- = 344.99 345.03 7200 

: : 347.06 345.10 Fad 

= = S504 la7 SAS oak ak 


RodT Ne.el2,. Runs Nose 


0.0 0.0 AQ. 73 = = 
ee Oe GZ Ag. 20 See 0.68 
15.0 0.59 48.69 OF Le OOF 
2260 0.60 Ay 1A ier al O66* 
29.5 0.60 AG S10 17.95 0265" 
5/22 0.60 OS aay Bane a 066% 
44.6 0.60 44,2 217.99 0.66% 
San) 0.60 aS. (i 2.50) 0.66% 
Be be pesles 42.99 SloiZ 0.65 
= = 41.44 50.17 0.62 
= = 41.44 Seal! 0, 62 
r = A0.92 2905 Oil. 
= = AQ.15 28.88 0.59 
“ = Se oi) Zo 10) O37 

Rolie Nos 12, Rum No. a2 
0.0 0.0 OS 742 = = 

2 0.40 8, 42 LOLS Leo Ol 
Lewd 0539 aaa S205 ee 
SANZ 0.39 OO? 46.27 Laas 
45.8 O59 80.29 60.65 1.38% 
Bale Del, 74507 65.78 sO 

= = 68.64 Don oo iio 
= = 65227 54.98 de: ina 

Roll Noi:el2% “Ron No.s3 
(mae) 0.0 SY = = 

a ee ee 0.40 Oil el 5 1420 eos 
Zh 0.39 5.24 lieu) pes 
34.2 0.39 88.06 he te) i.Aoe 
45.8 0.39 62.10 2252 L423 
Deke Bed FOe41. 66812 os 

= = GIGS 60.94 1.24 

‘ = 67.34 Sy) FOS Le 

= = 64.75 54.46 J 

= - 62.16 splay £06 

= = 60.87 50.58 Pe0S 
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t qt/oA x = APS. Ae (dP/ dx). 
Sec cn. cn i cm of H,0 cm.of H,0 am of H,0ém 
Roll Now 12, Run No. 4 
0.0 0.0 0.0 0.0 98.94 = = 
30 4.43 Z0 065 99.46 OF Sy Tei, 
60 Sew 14.0 0.63 100.49 UTES de Og 
90 15250 ZA 40 O65 101.01 55,08 besten 
120 Dial 28.0 OSS 1G 1653 48.59 oe 
150 Zeal Boek OF65 Layo O59 NG sior 
180 26.60 A2 2 E65 103.08 fon co Eeoae 
210 5L.04 49.2 O03 104.38 88.95 1.83" 
240 55 AT Bele Bere 106.97 G4. 72 Ue ees 
ROld Now UZe Run INO 65 
0.0 0.0 O20 a0 OF 5c = = 
30 ATS OFZ 0.48 OFS) LESS eo 5 
60 8.87 TSS 0.48 Te 247,50 48% 
90 15650 TAGS 0.48 88.58 - 50645 ees 
120 RES Shak 0.48 85.47 50S 1. 43* 
150 22n)7 46.6 0.48 81. 84 63.29 ae 
180 26.60 Baal Bale 78.74 67.96 1S 
20 See = = FOALS 652571 1. 3S 
240 Soaay zs Z (BOS 60.70 eae 
270 39.90 = - 67.86 Sees ES 
RO leMING na: lee UI NOseeO 
0.0 0.0 0.0 0.0 48.07 = = 
30 O22 Sisal 0.60 48.07 1250 0.76 
60 4.43 ex 0.59 47.66 S18 0.69 
90 6405 ey 0.59 AOS 5.97 ORNS 
120 8.87 1 Sek 0259 46.41 Se Olli 0.65" 
150 11.08 18.6 0.60 45.58 LOSL4 aeons 
180 P5250 LIBA® 0.59 44.76 eo 05005 
ZV Po e52 ZOuS 0.59 44.55 12 OFon. 
240 Dials 50a 0.59 43.93 LV G25 Onol. 
270 19.95 54a 0.59 A Sak 19.87 0207. 
300 AERA! 58 al 0.58 42.89 22229 moa 
330 ZAG 5S ALS 0.58 42.68 24.89 O5635* 
390 VAsPa SYA 49.4 0.58 41.96 29.96 OR ebe 
450 Bn AS) Bal Bale Ale65 SUL Sih 0.63% 
510 37.69 = o 40.61 29.83 Cow! 
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t qt/oA ne S AP, AP, (dP/dx)- 
sec cm cm. Mi cm of H,0 cm of H,0 a of H., Om. 
Roll No. 12, Run No. 7 
0.0 0.0 0.0 0.0 195.80 FE 
je 4.43 PIES 0.39 189.59 28.90 ia Ee 
30 a On 2560 0.39 180.26 65.75 re Oe 
AS 1350 5 SS 0.39 171.98 98.90 a 0as 
60 LES 45.8 0.39 162. 65 TS7e2 5 5.023 
TS Ce Bet. Bile L505 74 139.96 2.86" 
90 26.600 z =: 142.97 132,19 pak) 
105 31.04 = = Toke 7 120279 Wh 
120 Spiiys = = DOVES 116. 65 2008 
135 39.90 4 = Wee 25 SEDI Lapa 8)/ VRP lal 
150 44.34 = ) 2 TVS = 0 LOVE SZ IS) 
165 AS. 77 = = jb LOA 7S 2.14 
ROLL’ No: 123 Run No.3 
0.0 0.0 0.0 0.0 NO ey = = 
60 Tce Pee 0.74 18.86 = = 
120 3555 46 Dh al 19537 0.43 OnLy 
180 Date jee O75 19.68 JEt6 O23 
240 7.09 9.4 0.75 19.89 1e75 OLZ4 
530 O75 15,0 0275 19.99 2.44 On2e 
420 ee a Gea O75 20. 10 Syed Os 225 
510 Per Zee hee 70520 Syeere) ORE: 
600 1 Webes: 25 On S Ze ol Aledo Die 
690 20,59 2H Ona... 20.41 5.24 One 
810 23.94 Syiees 0.75 ZO O25 Ae ee 
930 27.49 Eyes Ory 2. Ons ieloye Deas 
1050 51.04 Ages! 0.74 21.96 9.18 OO. 25" 
1170 34.58 46.6 0.74 ZO Shdlecowh 0.25" 
1290 Soya io) Si bas One DEW ESS NEPA ea 
1410 41.68 Pete Dele 25/00 sean TA! Ci Bi 
1530 AS e722 = a 2A 55 10523 OFZ 
ROL INOW iS. Run eNO 
0.0 0.0 Ose 0.0 AB See | I = = 
60 0.44 ie 0.40 AWS = = 
120 0.89 2.6 0.34 LVRS) = = 
210 PASS AS 0834 A783 5 50 0.62 
300 Lee 6.6 0.34 AG. S31. oe Se O55 
390 Vipetes Bad O20 45.69 Sook 0.58 
480 byes 10.6 Uni55 45.69 eye 0.66% 
600 4.43 Res Oe 44.76 eas OF66* 
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TABLE F-7 ( Continued) 
a—06c6—0ennnnnan»@ma{={=9w9NS88NSSSSSS nS 
t qt/oA x , AP, APE (dP/dx )e 
cm of H,0 an. of H,0 cm of H,, 04m 


g 
g 





720 Sie owas O35 eras ale 9.49 G2 07" 
840 On4l cho Pat OFSS 44.24 Zoe OseZs 
1020 7.54 Zon On35 ANE A Aye dee! Ooi2* 
1200 8.87 7d es Ona 41.96 18.09 Ona 
1380 LOSZ0 SL a6 Uy 40.92 20.96 Oe als 
1560 DBixt55 oy | oz SRE 25.55 Coos 
1740 12.86 40.7 O6 82 38.54 26.68 0.70. 
1920 14.19 44.2 0242 So TAS ees) Oe70= 
2100 te toe AT OSs 505 51ess4 0.69% 
2280 16.85 ays Osa, 34.91 oe 0.68% 
2520 18.62 Bel: ely. Sel) SI06 005 
2760 20.39 = = 51708 297 Oi 0.60 
3060 22.61 : : ZOCa2 Dips Oe au 
3360 24.83 = . 28. 359 VAS SIs One 
3660 27.04 ANS) 2629 ORS 
3960 29.26 = = ZO 00 Zonelo OS 
4260 61.48 = - 25.90 Beles Ord) 
Rol No. 13, Run No. 2 
0.0 0.0 0.0 0.0 196". 62 5 = 
30 0.89 Bae Ory, 191.66 15.09 Mes 
60 ieee 10.6 Oy Lepaa2 30.10 5, 508 
90 2.06 losG Oe lay 183.89 Ay. 51 Saolee 
120 SNS 20.6 ea by NAS AWAS 60.41 Sos 
150 4.43 bees) Oleg Wy eS TSR UL Seba 
180 SSRs on ae Oey N68%.35 90.44 Balers 
210 On Z1 36.4 Ooely Voz. 13 1042 30 5,055 
240 Duo 41.7 ale 15.40 3.05 2eOIE 
270 7.98 47.0 ORLY, 149718 Ue. 2 raphe 
300 8.87 Sea Oa, Tepes 140.98 Za8en 
360 10.64 [By Bal. L30.54 TASH 2a05 
420 ak = = 118.10 NiG.05 25S 
480 14.19 z = 108.78 LOY 2.3L Pays) 
540 1.00 z Z Se ii 98750 2a0d 
660 LOT SL = = Oe) O0s/3 1262 
780 ZS = = 82.88 81.41 Lg66 
900 26.60 = Coe ge, 1.58 
1020 SU ribs z = Je LA 73.64 1.50 
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t qt/oA x & AP, AP. (dP/dx), 
sec cm. cn iy cm of H,0 an. .of H,0 an of H,04n 
RO LT Now gaye RudNG. 3 

0.0 0.0 0.0 0.0 AOS. - - 

60 0.44 es 0.29 49.73 : - 
120 0.89 a 0.30 49.73 0.98 0.98 
180 ip. oi 4.1 On 2 49.73 2 07 0.99 
240 SET Site 52 AO 7 ) 0.85 
360 2.66 8.7 0.31 48.17 5.04 Oss 
480 B85 1S (gg 47.66 7.28 0.77% 
600 4.43 14.4 jae 46.62 9.10 0.73% 
720 Senay 174 Q. 37 45.89 Nis 0.74% 
900 6.65 207 Pel 44.55 14.22 0.72 

1080 7.98 26.3 0.30 42.99 17.19 0.71% 
1320 9.75 2.0 0.30 41.44 21:26 0.71% 
1560 11.53 38.0 0.30 38.95 24.68 0.69% 
1860 13.74 Me 0.30 37.09 29.90 0.69% 
2160 15.96 52.8 0.30 33.98 aoe 0.66% 
2460 18.18 By ae 30.25 28.78 0.59 
2760 20.39 - - 29.01 27.54 0.56 
3060 2236 - : 27. O17 26.50 0.54 
3360 24.83 - : 24 25.67 6.52 
Roll No. 13, Run No. 4 
C. 0 0.0 0.0 0.0 170.94 = 3 

30 0.89 4.9 0.18 183.89 22.98 7.92 

60 eae 9.9 0.18 18289 A227 pst) 

90 2.66 14.9 0.18 180.01 53.68 4.16% 
120 5.55 19.5 0.18 16. 12 6570 Be oe 
180 5.32 29.3 0.18 164.47 87.95 3.22% 
240 Wess) 5902 0.18 149.70 AO 4d. ole 
300 8.87 49.2 0.18 134.68 126.98 2.69% 
420 12.41 Baie Bale 110.85 109.38 2.23 
540 15.96 - ss 94.02 02855 1.89 
660 19.51 - 5 84.18 82074 1.69 
720 23.05 - ; Tet 76.23 1.56 
900 26.60 - : OS 69.76 1.42 

1020 301, 15 = 3 67.86 66.39 1.35 
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TABLE F-8 
FRICTION FACTOR AND REYNOLDS NUMBER. OF 
THE FLOODED ZONE 
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Run 
q , We pee NRe fy 
3 cm. oO H, /cm. 
No. No. an’/sec. poise — 5 zi 
Pred. EXD. x10 x10 
Water Displacing 100 cs Dow Coming,Large Beads 
8 1 1.5885 0.22 0.49 0.07 0.038 To Qu27 
8 2 Sale 0.24 0.49 is 0215 38% Sif 0.26 
8 eS 6.354 Oa29 0.49 0.30 0.43 77.14 0-19 
8 4 6.354 Q.2 0.49 0.30 a5 77.14 O15 
8 5 0.6354 0.25 0.49 0205 Ore ey 0.88 
Water Displacing 100 CS Dow Corning, Snall Beads 
10 7 1.5885 OnG 0.49 0.60 0.60 Lies, 1.04 
10 8 0.6354° 0.17 0.49 0.24 QZ? 3.88 2093 
10 10 Set WHY: Oe 0.49 gh78) 13359 19.33 0.60 
11 7. 6.354 0.20 0.49 fae GrUr 38.76 O55 
Glycerol Solution Displacing 160 CS Dow Corning, Small Reads 
12 1 1.5885 0.60 0.61 0.75 0-07 deeds ic 
12 Z See 0239 Qua/ 1.39 1.49 16.66 0.65 
jx See By Oes9 UgeWs 39 1.46 16.66 0.63 
12 4 Boner 0.63 0.74 ist 1.88 1d ee) 0.82 
12 ei, 0.48 0.61 1.49 1.43 ieee ul) 0.62 
ie 6 PeSSS5 5 0E5o 0.61 On75 0.62 7.18 1.08 
12 7 6.354 0.39 0.61 2.98 5201 31.14 ORS 
Lee 8 0.6354 0.75 0.62 0.30 Dees 3.06 2.00 
0.05% Sep. Displacing 100 cs Dow Corning, Large Beads 
10 a 6.354 0.38 0.58 O55 0555 5a CAS 
10 3 1.5885 0.56 0.64 0.10 Oe12 14.76 0.84 
10 4 6.354 0.38 0.58 0.55 0.43 65-417 D719 
10 See .585 0.38 0.56 0.84 god 168.74 0.04 
10 6 6.354 0.38 0.58 Oi35 0.41 Goel) 0.18 
0.05% Sep. Displacing 100 CS Dow Coming, Small Beads 
11 Ariens 77 0.38 0.56 ey Riis 16.96 0.94 
11 5 6.354 0.34 Osos 2.69 See 34.55 0.45 
11 6 1.5885 0.47 0.58 Oat lS 8.19 1.96 
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TABLE F-8 (Continued) 
q “ He (aP/dx)_ Noe 
é Fr 
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0.6354 0.18 hs 2.03 Sed USD 
Pusher 500 Displacing 1000 cs Dow Corning, Small Beads 
Oels 6859 0255 See 0.70 0.69 8. 30x10" 
0.6554 3 017 ee de 2.64 yes Q@.35 


So Oo re 
° . e ° 


yabd «od #250 






ius a p20 -eaned 0 é 
| | EG BAEOM eee 
We: prefer: | oe etl. Oe y Teil a< 6 

’ ; Le bee. VG Bie ; dd « 


f .0ppeesy oy 
| hie b4248 rn 
9 7 6 — 





© 





200 
F.5 Volumetric Recovery Efficiency at Breakthrough 
Volumetric recovery efficiency at breakthrough was calculated 
either from the slope of the eye-fit lines in Figures (V-3) to (V-18) or 
directly from the time of breakthrough measured by a stop watch. Close 
agreement between the two methods is evident from Table (F-9). The data 


points of Figures (V-19, 20) are based on the first method. 
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In Figure (V-19) in addition to the experimental data of the 
present study the data points of other investigators available in the 
literature are also included. These data points and the pertinent 
experimental conditions involved are briefly discussed below. 

ies Everett, et al. (69) studied the effect of viscosity ratio on 
displacement recovery for miscible as well as immiscible displacement. 
The waterflood experiments were performed on a linear system made of 
cemented Silica sand grains having diameters between 65 and 270 mesh. 
The permeability to air of the various linear portions was reported to 


be 2.77, 2.44, 2.54, and 2.99 darcies. The average porosity of the 
5 


Syocem was 2/.2 percent with a total pore volume ot 980°am. The over- 
all colum was 115.8 cm long and 6.35 cm inside diameter. The oils 
used were standard grade lubricating oils. By varying the type oil and 


the temperature during displacement the desired ratios of oil viscosity 
to water viscosity were obtained. The displacing liquid was always water. 
All of the waterflood displacement tests were performed under a constant 
pressure drop of 20 psia. The data points reported in Table (F-10) are 


taken directly from the Table II of their paper. 
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BREAKTHROUGH RECOVERY vs. VISCOSITY . 
RATIO OBTAINED FROM TABLE II OF 


EVERETT, et.al. (69) 


Hol Hy Ep () 
1.04 81.6 
2.41 EG 
Rao A128 

15.49 41.8 

48.53 San 

98.6 ZIZ0 


2. Engelberts and Klinkenberg (2) studied the displacement of 
oil by water from a porous formation consisting of grains of sand between 
20 and 30 mesh with a permeability of 200 darcies at a porosity of 38%. 
Kerosene and various white oils were emploved together with sec.butyl alcohol 
or, for the more viscous oils, methyl ethyl ketone, to simulate the oil phase. 
Distilled water saturated with alcohol or ketone represented the water phase. 
In Figure 5 of their paper, they plotted the breakthrough recovery versus a 
dimensionless group dp gk/Un for three different viscosity ratios. isis 
dimensionless group which was obtained from dimensional analysis is the ratio 
of the hydrostatic pressure gradient, dog, to the flow potential gradient 
Un, /K in the water phase. The breakthrough recovery efficiency at 
“Apgh/Un, = 0 was read off from their Figure 5 which corresponds to the 
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absence of any gravity effect in the present study. 
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ABLES per 


BREAKTHROUGH RECOVERY vs. VISCOSITY 
RATIO OBTAINED FROM FIGURE 5 OF 


ENGELBERTS, et al. (2) 


Baal. Ep (3) 
1 83 
4 49 

24 28 


3. Jones-Parra, et al. (70) studied constant rate displacements in 
water wet systems. The porous matrices used were constructed from 2.5 inches 
inside diameter lucite tubes packed with unconsolidated glass beads. The 
packed tubes were of three different lengths, 15.0, 254) aide Scecnaee 1 ne 
beads, of 140 - 170 mesh size were almost perfectly spherical. The permea- 


bility of the cores were reported to be approximately 7 darcies and porosity 


35.5%. The displacing liquids were sodium chloride solutions blended with 
varying amounts of glycerol to produce different WiscoSiticess. | Inesdrspiaced= 


phase liquids were made up of inert water-white petroleum derivatives together 
-with sufficient carbon tetra-chloride to adjust the liquid specific gravities 
to unity. Their results for floods performed on the core of 25.4ycm, Jength 


with unfavourable viscosity ratio are reproduced in Table F-12. 


a) 
> Lal 
: } 
A . 
ray) st ait di y=) a Lvste ae 
ie by ae * 
estarrin $4 “tae biaptt mat 3 


Lad 


Plgged oa tte gato it ne Beno ay 28h 3 70%. poy eae 


ar = 





i? wey UR eae 
Tet, GOA OETA pio 
2 , da : 
wh 
a be 
I . 
p 
1 - 
- =F" erick mA 

















= 


: * hot 
4f scot Sem eo sees > lente aie 


eu a 
. . « ms 
Sa U xy nee” -«z . aq i es 
é ; oy yx. 2 ig rosin 2B 
: ; =i ph : 
= » & . o's » ys : ae 
; pat. 4a reSrib te86) Io ae richagt 


e | o ° rT af 
aaal se a cot OL — Gat Ss 

3 <9 0 ora 
rrp wi Lint Grew STC sai A 
tia mothe exow abi eu suigedh nif 


“ be < o 


ree 


2 \* , ~s a * -i 
Va}'32 Bb taebctia of In ante te asm ES es 

- * - ] , 7 ie ‘ir 7 a 
ri) aarti Jrodi io wm anne Le ee 


rentin ae, ob beaido- sist Geode drat 
lett hee is ; 
jaar de ee 
ee ee 7 Peer ms. 


Mo Bel” S ‘ 
ne > § oi 7 i ryt : 










- - = 








206 
TABLE F-12 
BREAKTHROUGH RECOVERY EFFICIENCY 


OBTAINED FROM TABLE III OF 


JONES-PARRA, et. al. (70) 





0.412 OE 0.92 Bas 59.5 
0.271 Cyd 0.92 5053 58.4 
0.098 ae95 0.92 Ba55 63.9 
0.030 Sey 0.92 aS 3 is 
0.545 ae 0.92 BSG 60.5 
0.598 1.85 0.92 Det 75.8 
0.294 1.85 0.92 Paoli Parsi 
0.151 1.85 0.92 201 64.6 
0.071 185 0.92 Pel 68.7 
0.075 135 0.92 2.0% 66.0 
0.043 1.85 0.92 oO) 61.2 
0.573 2.4 0.92 2°61 7g 
0.280 2.4 0.92 2.61 6255 
0.196 2.4 0.92 2261 74.8 
0.129 a24 0.92 2.61 68.7 
0.063 2.4 0.92 Va ae0 
0.007 eA 0.92 2.61 550 
0.069 2.4 0.92 2.61 68.2 
0.161 2.4 0.92 2.61 7728 
0.153 Dea 0.92 2.61 Haak 
0.473 ee 0.92 5.65 49.2 
0.259 cee 0.92 5.65 Soul 
0.188 Ge2 0.92 5.65 50.3 
0.126 cae 0.92 5.65 62.6 
0.059 me 0.92 5.65 64.0 
0.02 52 0.92 5.65 63.0 
0.129 5.2 0.92 5.65 56.5 
0.003 c2 0.92 5.65 59.5 
0.745 2 0.92 5.65 Ava 







Ct) at 
ac Tet TO ae aaa 
a I Go Ys 7 “Eter Me OE 
(3 cote Snape, 


207 


4. Croes and Schwarz (51) reported results in graphical form of 
a series of experiments carried out for measuring the efficiency of the 
water-drive process affected by the oil/water viscosity ratio. The 
model used in the experiments was a tube of 1m length and 6.2 ™m 
diameter, filled with umconsolidated sand. The sand was 20-30 mesh 
having a permeability of 2 x 10° md and a porosity of 33 percent. dhe 
displacing phase was water and the displaced phase an SAE 30 lubricating 
oil blended with Kerosene to obtain the different viscosity ratios. 
The residual oil saturation was reported to be 15%. In this study they 
plotted the viscosity ratio vs. cumulative oil production as AC ia OL Bey 
place Bag for various values of gross production including breakthrough 
recovery. The data points presented in Table F-13 were read off this 
- graph and corrected for the residual oil saturation by dividing x by 


0.85 - Refer to Appendix D. 
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TABLE F-13 
BREAKTHROUGH RECOVERY vs. VISCOSITY 


RATIO OBTAINED FROM FIGURE 2 OF 
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F.6 Saturation Profiles 

Saturation profiles were determined by the method of light. scanning 
which was discussed in Section (IV-5). The difficulties involved and the 
error introduced during the course of converting the scanning readings into 
the saturation data are thoroughly analyzed in Appendix G.3. The data 
points of Figures (V-21) to (V-37) are presented in Tables (F-14, 15). 

The net number of counts obtained from scanning of line i is 
presented as Oa and the water saturation of the plane corresponding to that 
line is taken as She a Assuming linear dependency of Cy and Sa? one can 


} 
Zs 


write, 
Swyi = Neg! Ni et) 


where ee is the conversion factor. Since the displacement tests were net - 
performed originally for the purpose of the saturation measurement, (the 
amount of the light and colour was not controlled) precalibration of the 
scanning instrument for water saturation was not possible. This means that 
the vaiue of ele was not known in advance and it was necessary to use 
material balance information to determine the ea indirectly. 

The area under the curve of Sw v.s. x at or before breakthrough is 


related to the cumulative volume of water injected by, 


re 
qt. = i S_ dx (F-2) 


O 


or writing in difference fom, 


t nN 
Se i Sy RS ya (F-3) 
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where AX, means 5 447%4)> i = 1 corresponds to the beginning of the bed 
and i = n the position of zero saturation plane. Substituting (F=-1) into 


(F-3) and assuming Ka ; constant for all saturations and/or positions, the 
> 


following expression is obtained for the conversion factor Ka i me Ko 
> 
t nN 
G = Ad TJ Se 8 we es, 


oe ate eee NL (BONS) 


ihe influence, on themverticale bar on ve for the case of more 
uniformly swept frames was significant and the procedure to correct this 
influence is given in Appendix G.3. These data points belong to Figures 
(V-31) to (V-37) which are presented in Table (F-15). Table (F-14) includes 
the datalpoints of Figures (V-2i to V-30) in whichéthe effect of the vertical 


bar was insignificant so that no correction need to be considered. 
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TABLE F-14 


WATER SATURATION DISTRIBUTION FROM 
FILM SCANNING - - EFFECT OF VERTICAL 


BAR INSIGNIFICANT 








ran ia 





x, Cn AX; ,cm ote ae (x) 
Roll No. 8, Run No. 1, qt/¢A = 10.71 cm 
22.05 AUS 3680 Si. 
Sev. Lee Zel, 3740 505 10 
4.60 IRS 2920 Aer ell 
Ba Gs LeeoD 2620 AQ, ho 
oF 09 Be 55 2610 AQ..59 
10. OY, igeraes 2880 44,79 
12.24 Vaya 2420 37.64 
E35 Su 2 Lei 35.0) 
14.79 L528 2150 33.44 
16.06 127 TAL) Sao 
18-61 Zs SD 1720 hoe Ne 
ZA ELG DoS 1820 28.50 
25. bO 2. 4 1280 INNS Sit 
26.25 5.035 1180 UiteWeeie) 
5135 Brey” 1040 someday 
5a. 90) UNS 5/0) Sees 
36.44 Js ill 550, rel 8S 
55.99 Pas 170 2.64 
41.54 2.02 240 Srl 
45.36 3.82 180 2.80 
PCy 4 ax; = 68867.20 
Roli Noz. 8, RuieNo. 3, ot/ok = 158635ecr 

tape) 2.94 4080 54.62 
A ehid 1s) 3960 pe 0 
5.49 UZ 4060 OS: 
Gen 18S 0.64 3920 52,48 
aS L 5.28 3840 San 
10.58 RZ S200) 47.49 
11.86 Te Zs 5550 44.85 
$355.23 Vangel 3840 eee a 
14,41 Les 4100 54.89 
15.68 Lei 3880 51.94 
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TABLE F-14 (Continued) 


oNee 

il 3890 
55 3630 
55 3440 
oe 2810 
Zh 2360 
Dal 227.0 
55 2160 
28 1930 
2 1690 
55 1090 
54 800 
55 520 
28 630 
27 660 
55 160 
Zu S50 
27 230 
55 360 

Ax, = 16747250 

ROLI No. 3) RungNoaes, dt/oA = 13.11 
91 4910 
29 4710 
28 4400 
56 4340 
28 4800 
28 3890 
29 4320 
28 3880 
28 4050 
28 4230 
56 2900 
owe 5110 
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APPENDIX G 


ERROR ANALYSIS 


In this Appendix the major sources of error will be pointed out 
and discussed briefly. Since the scope of the present study is not to 
correlate the experimental data, it is not intended to perform an analytical 


study of error propagation to achieve the variance and distribution of 


random error. 


G.1_ Pressure Drop Data 


The error in the pressure drop data of the single phase flow is 
dependent on the transducer and flow rate measurements. The transducers 
and the pumps were calibrated in advance. To eliminate the possible error 
due to the head of fluid the transducers were always placed at a level below 
the lowest pressure tap so that an equal head was always applied on both 
sides of the transducer. 

Some approximations are involved in measurement and calculation of 
the pressure drop and pressure gradient of the mixed zone for the displace- 
ment tests. There exists a small difference in specific gravity of the 
displacing and displaced fluids. For example in the case of water flood the 
difference in specific gravity was 0.04. Therefore, a change in fluid head 
is being created by advancing the front plane. This has been calculated and 
subtracted from the total pressure drop reading to obtain the pressure drop otf 
the flooded zone reported in Table (F-7). 

The pressure gradient of the mixed zone may be in error because of 


uncertainty in determining the exact length of the mixed zone. The length of 


the mixed zone is approximated as the length of the flooded zone, Xo: 
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The actual length of the mixed zone will be somewhat less than Xo 1tuene 
plane of Sy = 1.0 progress significantly past the inlet screen. Figures 
(V-21) to (V-37), however, show that the latter is not the case even for the 
moderate viscosity ratio; therefore, the error due to this problem cannot be 
Significant. 

The location of the first pressure tap introduces a different kind 
Greerror. ihe first pressure tap is located approximately two centimeters 
above the inlet screen; therefore, the first two centimeters of the flooded 
zone are being excluded from the pressure drop measurement. The: error 
contribution of this segment is larger at the beginning of the displacement, 
and as the frcntal plane proceeds its effect diminishes. To minimize this 
effect, the initial pressure history of the displacement is not being 
considered in calculating the pressure gradient of the flooded zone. In fact 
the pressure drop reading is taken from the time that the frontal plane has 


proceeded at least 10 am. into the bed. 


G.2 Frontal Plane Advance and Breakthrough Recovery 

The data points in Figures (V-5) to (V-18) show some scattering. 
From the procedure described in Appendix F.2, it is clear that any measurement 
errors are small in comparison with the scatter. The scatter is attributable 
to the local non-homogeneity of the bed and also the random appearance of the 
fingers. For example a finger may be created due to the local non-homogeneity 
but does not necessarily grow and become prominent finger, or a prominent 
finger may stop growing and some other new finger may take, over. This creates 


a statistical kind of scatter in instantaneous measurements of frontal advance 


rate. 
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Breakthrough recovery was determined by the slope of the eye-fit of 
Figures (V-3) to (V-18) as well as from the breakthrough time measured by the 


stop watch. Close agreement was always obtained: Refer to Table (F-9). 


Gao acawuratvon Profile 

It is useful to recall that saturation profiles presented in Figures 
(V-21) to (V-37) were determined by scanning individual frames of the film. 
The procedure of constructing saturation profiles from the scanning data is 
given in Appendix (F.6). As described in section (IV.5) the measurement of 
diffused light along a scanning line was performed on a comparative basis with 
two frames of the movie film from each run: one before the displacement, when 
the bed was full of the oil, and the other at breakthrough. The concentration 
profile could be influenced by the steel bars appearing on the photographs and 
by uncertainty in matching position on the frames. 

In Section (IV. 1.2), it was mentioned that there existed four steel 
bars on each side of the bed, three positioned laterally and one vertically 
along the bed. The bars appeared as bright lines on the photographs. These 
bright areas on the photographs gave minimum microdensitometer readings and so 
the lateral lines could be easily located by doing a rapid vertical scan and 
noting where minimum counts occurred. Position of the bar was then deter- 
mined from reading of the position micrometer on the film carriage. Alter 
the positions of these lines were determined, the equivalent positions on the 
two frames were chosen. 

Scanning was generally performed at 0.5 to 1mm intervals. At the 
positions of the lateral bars the intervals were increased to approximately 


3mm so that any effect of the lateral bars could be eliminated. 
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Some error is introduced due to the fact that near to the bars the 
reference frame does not diffuse the same amount of light along different 
scanning lines. This situation arises because of the shadow and reflection 
of the steel bars. This error could be avoided if exact match of equivalent 
positions of the two frames were possible. An estimate based on an uncertainty 
of 0.5 mm in position showed that the error could be as large as 10% at 
positions very close to the tae Actual error from this source is probably 
very much less than this. 

Although the laterally positioned steel bars were excluded in 
scanning, the vertical bar could not be eliminated from the scanning. Since 
the bar always appears the same on both reference and frames after flooding, 
the instrument would always indicate a zero saturation at the position of the 
vertical bar. This effect is insignificant for severe fingering because the 
chance of hiding a finger under the bar in this case is small, and in fact, 
more or less random. 

For moderate viscosity ratios, for which the bed is being swept more 
uniformly, some portion of the displacing fluid is always hidden under the bar. 
For a given scanning line the error in coumts attributable to the vertical bar 
could be as large as 25% of the true counts, depending on the concentration 
beneath the bar. Note from Equation (F-5) that an increase of counts by 25% 
does not necessarily increase saturation by 25%, but always less. The most 
influence of the bar is at the beginning of the bed and the least at the end. 
That is at the position of ns 1, hundred percent of the influence exists 
while at the position of Soe 0 there is none. The existing value of Sy at 
each position gives the percentage of the total influence of the bar at that 


position. To avoid a trial-and-error procedure, the theoretical prediction 
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of Sy was chosen as the weight factor of this error. Therefore, Sw was 
read off from the Figures (V-31) to (V-37) at each position and multiplied 
by the maximum influence of the bar (which was 25%) to obtain the error 
contribution of the bar at that position. The corrected values as well as 
the original data are presented in Table (F-15) of Appendix F. liters 
evident from Table (F-15) that 25% increase in the number of counts obtained 
from scanning changes the saturation data points approximately 10%. There- 
fore, in Figures (V-31) to (V-37), even the effect of vertical bars in its 
extreme limit does not change the shape of the experimental saturation 


profile singificantly. 
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APPENDIX H 


EEPEGEIVE VISCOSITY OF THE MIXED’ ZONE 


Let us assume an element dx of the mixed zone shown in Figure 
(H-1). Assuming that Darcy's law holds, the pressure drop, dp of this 


element is given by, 


qu, dx (H-1) 


where Le is given by Equation (11-34). 





FIGURE H-L § ScHEMATIC REPRESENTATION OF THE 
lxeED ZONE. 
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Total pressure drop of the mixed zone is simply obtained by integrating 


Equation (H-1) between xy and Xo? 


x 
AP ‘aL 0 es dx (H-2) 
ih 


Defining an effective viscosity Ue for the mixed zone, Darcy's law can be 


written as, 


ADA 
AP = ran &, = x4) (H-3) 





From Equations (H-1) to (H-3) therefore, 


ts Xo 
He = es aes dx (H-4) 
O i 


The expressions of Xo and X, are obtained from Equation (11-50) by 





substituting Sy = 0 and unity respectively. The expression of Xy has 


already been given by Equation (11-54) and that of XG is given below, 


_ qt Mtn M : 
X) © (Mei) ¢A ct 


Therefore, 


x - Xy = qt %n M (H-6) 


ower hoeke mele Bp (Ih sweety . ’ Ye 
(cm falta temaded (140 mere: a 


| 
c 
i 
~— 
‘ 
~ a 


rr Gina tity O-= 2 galaenitedaees 
’ 


» 
9 








Al s) 
oe. 


Be: é 


wiv gabzoslte og gabe 


en nstsie 


ar a) 


we o jhe] aorqevedia 


otteams i erie aod : 


nN 






230 
Equation (H-4) can be integrated either directly by substituting He from 


Equation (11-34) and Sw from Equation (11-51) or, as a second approach by 
transforming the variable x into Sy 20Cc the second approach dx is obtained 
from differentiating Equation II-50. The advantage of this appears when 
the procedure is being considered for polymer flood in which case an explicit 
relation for Sw in terms of x is not possible. 

After straightforward algebraic manipulation, the final result for 


Newtonian fluid is obtained as, 


no “> (11-33) 


which means that the effective viscosity of the mixed zone is constant, 
independent of 2S and is simply equal to the arithmetic average of 
the displacing and displaced fluid. 

Due to the complexity of Sw as a function of x shown by Equation 
(11-66), integration for the non-Newtonian fluid is not as simple as 
Newtonian. Equation (H-4) remains the same in which Xo and X, must be 


obtained from Equation (11-66) and Ue is substituted from Equation (II-62). 
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Equation (II-66) is differentiated with respect to Sy to obtain dx, 


da a dn 
qs (qt 2n Mp oA) Le eae ") Bee fou (1M ‘ a) &, | 


Pes) a] boy “ 


(H-10) 





Cie nee M 
x p (H=1 1) 


ds, ae ) +S ( 
WwW W 


By substituting Equations (II-61, 62) and (H-7 to 11) into (H-4), on 


simplification one obtains, 


yl. : 
roe (4 eae 's) (H-12) 





byt 


fouesy aw beast eb 30 ei (ode a ti 
a! \ oe 2 é : -_—<, 


Mg 
— 


fd-tl) corteopl yt mevrg ain ered 


. r : a. 
6 gh0-LT) enotj pupa qnisutisadue: 









(> ' 
m/f ‘. 9 il 

A a .. al ——. 
( wi rh 3) 5 “ 
rT. 


ma me (nb) a 

: Un’ &/> 

Ve AE Sat ph 

Tr: . i) we . 
¥ (we A nee 


le 


h > F _ 


aoe 





= 5 eae 1 (2 aa ae Ce a eee 
(vT H) Mop W ( Ww _(u-1) (AZ) 2 (u-T) (Kuz) 2 ( 


eee RAC aes) eee 


(Se 4 Su - U) 


d 
d (T - | Wu 
(€T-H) a i es ae 





u 


sutsn pure 


:0190UM 





233 


Oe 





























ose W- T) mew { - 
“sp }(, 9 W) (L-DAz -eC Wer) Cs 
K/(1+ S)Z- A/(1+ $)Z- 
(ST-H) 
A 
‘ C/T W- 1 moe . 
- Sa oO 
sp } ( wa) | WO-D+C ed] + Fo Ge) 9+ 4g Cm) Cut) oJ 
K/(t+"sz) - K/(t-"S) Ges) - A/ (I~ 'Sz-) 





234 

As a check it is possible to show that at n= 1 Equation (H-12) 
eventually produces Equation (11-33). It is clear that the integration 
of Equation (H-12) is extremely tedious if possible at all. Therefore 
an alternative method is proposed to approximate the viscosity of the 
mixed zone for the polymer solution. 

The apparent viscosity of the polymer solution is estimated 
from the shear rate encountered in the bed. This estimated apparent 
viscosity is used in Equation (11-33) in place of He to approximate the 
effective viscosity of the mixed zone. 

The wall shear rate for the packed bed based on:the capillary 
tube model can be shown (40) to be equal to, 


onal 1s 8) 


os oat) (H-17) 
W 4n x AUS0KS 





Therefore, apparent viscosity based on wall shear rate, n_, is equal to, 


= a eae eee 
ae a lw) 


In Table (H-1) wall shear rate and the apparent viscosity, ny» are 
reported for different flow rates and polymer solutions applied in this 
work. The effective viscosity and the predicted pressure gradient of the 


mixed zone are given in Table (F-8). 
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TABLE H-1 


APPARENT VISCOSITY OF THE POLYMER 


SOLUTIONS BASED ON THE WALL 


SHEAR RATE 
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TABLE H-1 (Continued) 





ae EE 24.12 1.04 


0.2% Pusher 500, Small Beads 


0.15885 


ac 
0.6354 4.1 


oo) 
ta CA 
fH pA 
HN 


0.2% Separan AP-273, Small Beads | 
| 




















